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Zur plastischen Knickung gerader Stiibe. 
Von A. Pfliiger. 


1. Einleitung. Von Shanley wurde erstmalig 1946 darauf hingewiesen 1, daB die bekannte 
Engesser-Kdrmdnsche Theorie der plastischen Knickung, die bereits um die Jahrhundertwende 
entwickelt und seitdem stets als einwandfrei angesehen wurde, doch einer wesentlichen Ergan- 
zung bedarf. Es ist das Verdienst von Schleicher, in Deutschland die Arbeiten von Shanley be- 
kannt gemacht zu haben ?®. 

Shanley gewinnt seine Erkenntnisse vorwiegend durch anschauliche Betrachtungen, die er 
an Hand von Versuchsergebnissen anstellt. Rechnerisch werden die Verhiltnisse lediglich an 
einem sehr einfachen Modell verfolgt, das aus einem beiderseits gelenkig gelagerten Knickstab 
besteht, der in Stabmitte ein elastisches Gelenk hat und im iibrigen vollkommen biegestarr ist. 
Es bleibt danach die Frage offen, wie die Lésung des Problems in allgemeineren Fallen im einzelnen 
aussieht. Diese Frage erscheint um so wichtiger, als die Ausfiihrungen von Shanley bzw. 
Schleicher bereits angefochten worden sind *. Im folgenden soll die vollstandige exakte Lisung 
fiir den Knickstab konstanten Querschnitts entwickelt und fiir zwei Sonderfalle, namlich fiir den 
Rechteckquerschnitt und den ,,Zweipunktquerschnitt“*, zahlenmaBig ausgewertet werden. Dabei 
wird sich ergeben — um ein wichtiges Resultat der folgenden Betrachtungen vorweg zu nehmen — 
da8B es sich nicht darum handelt, irgendwelche neuen Hypothesen mechanischer Art gegeniiber 
der alten Theorie von Engesser-Kdarmédn einzufiihren, sondern einfach darum, das Problem mathe- 
matisch vollstandig zu beschreiben und bisher vergessene Lésungen zu der seit langem bekannten 
hinzu zu nehmen. 


2. Einzelheiten der Aufgabenstellung. Voraussetzungen. Nach Abb.1 sei ein beiderseits ge- 
lenkig gelagerter Knickstab mit der Koordinate x langs der Stabachse, der Durchbiegung w, 
der Stablange / und der Belastung P betrachtet. 

Die Querschnittsform sei fiir alle x konstant. Es sei nur das ,,Eulerknicken“ in der Zeichen- 
ebene, aber kein Biegedrehknicken untersucht. Fiir den Querschnitt mége gelten (vgl. Abb. 2): 


z Koordinate vom Schwerpunkt S aus gerechnet, pes p 
c z-Koordinate der untersten Faser, i ae SE 
F Flache und ml. i wil 


J Tragheitsmoment in Bezug auf die beim Knicken in 


. Abb. 1. Skizze des Knickstabes. 
Betracht kommende Querschnittshauptachse. 


Zur Erlauterung der beim Knicken im plastischen Bereich fiir das Spannungs-Dehnungs- 
diagramm zu treffenden Annahmen diene zunachst Abb. 3a, in der fiir einen beliebigen Ab- 
stand z die dort herrschende Druckspannung —g, in Abhangigkeit von der Stauchung —e, dar- 
gestellt ist. Bei einer Laststeigerung in irgendeinem Punkte A des Diagramms kommt die stark 
ausgezogene Kurve mit der Neigung E’ in Betracht, bei einer Entlastung die gestrichelt ein- 
gezeichnete Gerade, die der Tangente im Nullpunkt parallel ist und die Neigung E hat. Wird 
weiter die bei der Lisung von Verzweigungsproblemen iibliche und im allgemeinen berechtigte 
Annahme hinzugenommen, dafs die Verformung des Grundzustandes, d.h. die Zusammendriickung 
des geraden Stabes, vernachlassigt werden und das genaue Spannungs-Dehnungsdiagramm im Ver- 
zweigungspunkt durch seine Tangente ersetzt werden darf, so ergibt sich der in Abb. 3b dar- 
gestellte Zusammenhang. 0, ist dabei die Spannung des Grundzustandes, wahrend e, nunmehr 
die beim Knicken zusatzlich auftretende Dehnung bedeutet. Die durch Abb. 3b gekennzeichnete 
. Abhangigkeit gibt die gewohnlichen Annahmen der Engesser-Kdrmdnschen Theorie wieder und 


1 F. R. Shanley, J. aeronaut. Sci. 13 (1946), S. 678 und 14 (1947), S. 261; Proc. Amer. Soc. Civil Engr. 


75 (1949), S. 759. 
‘ Hs Sokivicher, Bau-Ing. 26 (1951), S. 139 und S. 197; Abhandlungen aus dem Stahlbau, Heft 10 (1951), 


SLB: 
3 W. Cornelius, Bau-Ing. 26 (1951), S. 255. 
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wird fiir das folgende als exakt giiltig vorausgesetzt. Unbefriedigend bleibt bei diesem Ansatz, 
daB er weder die Kriimmung des Spannungs-Dehnungsdiagramms, noch die sich bei wiederholter 
Be- und Entlastung ergebenden Verhaltnisse zu erfassen gestattet. Dem wirklichen Werkstoff- 
verhalten wird also die Rechnung nur fiir hinreichend kleine Verformungen gerecht. Bei plastischen 
Verformungen entsteht im Querschnitt ein Belastungs- und ein Entlastungsbereich. Diese Be- 
reiche mégen nach Abb. 2 die Flichen F, und F, haben und seien durch die Grenzlinie z = c — s 


von einander getrennt. 


Abb. 2. Querschnitt mit Spannungsverteilung. Abb, 3a, b. Spannungs-Dehnungsdiagramme. 


3. Ableitung der Differentialgleichungen. Wird zur Abkiirzung 
E 
aan (1) 
gesetzt, so lauten die Voraussetzungen von Abb. 3b 
Oz =O) + 1D. fiir <n 
G) 64 Pn Ee, > ‘a 


Durch Integration der Spannungen bzw. Spannungsmomente iiber den Querschnitt ergeben 
sich in bekannter Weise Langskraft und Biegemomente. Dabei ist es zweckmabig — abweichend 
von Shanley — so vorzugehen, wie es in der Statik des Stahlbetonbaus iiblich ist. Hierdurch 
wird der mechanische Sinn der Ansatze besonders deutlich. Dementsprechend sei einheitlich 
mit dem Modul E’ gerechnet und der Entlastungsbereich fiir den beim Knicken zusatzlich auf- 
tretenden Spannungs- und Verformungszustand als n-fach wirkende Flache angesetzt. Dazu seien 


folgende Griffen eingefiihrt: 


(2 a, b) 


i} 


Ex 
Ex 


I ee ——id ee oe Om shan | 
1 4 fu * 
ea zdF +n | zdF\ , 
A( [xar-+n fear), 
(F,) (F,) ( (3a, b, ec) 
J. =f (e—e?dF +n | (s—e)dF. 
(Fi) (Fs) 


F, ist die .,ideelle** Flache, J, das ideelle Tragheitsmoment und e die Koordinate des Schwer- 
punkts S, der Flache F, (vgl. Abb.2). Ferner sei 

é Dehnung bei z = e, gemessen vom Grundzustand aus, 

N_ Langskraft, 

M, Biegemoment in Bezug auf die Hauptachse z = e des ideellen: Querschnitts. 

Die Annahme vom Ebenbleiben der Querschnitte ]aBt sich in der Form 


dat (4) 


schreiben. Mit (2), (3) und (4) und unter Beachtung der Tatsache, da das statische Moment 
der ideellen Querschnittsflache F, fiir die Achse z = e und das der wirklichen Flache F fir die 
Achse z = 0 verschwindet, ergibt sich dann fiir die Langskraft 


f rt / a ) 
as [oar = | 6) +E’ e,—E' (z aa |F | | 0) +n E's, —nE' (eo) ol aR 
7 dx 
(F) (F,) (F.) 
=o,F+E' F.e, 
und entsprechend fiir das Biegemoment 
WU pasemmey aes fh h walt ates A 
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Mit o) F = N, erhalt man also fiir die SchnittgréBen das Elastizitatsgesetz 


M, = N,e —E’ J, 


APs | (5 a, b) 


das in seiner Anschaulichkeit keiner weiteren Erlauterung bedarf. Befinden sich alle Fasern eines 
Stabes ausschlieBlich im Belastungsbereich, so ist e = Cl, =k and So J. 
Die Differentialgleichungen des Problems ergeben sich aus den Gleichgewichtsbedingungen, 
a - sie hr auf lineare Glieder der beim Knicken auftretenden Zusatzverformungen 
ie Gestalt 


NP 


M,= P(e+w) | Were 
haben. Im Grundzustand sei P — Po; es wird also 
| N, = —P,. (7) 
Mit (6) und (7) und der Abkiirzung 
AP = Pe Pp. (8) 
folgen aus (5) die gesuchten Differentialgleichungen zu 
E’F,2, + AP =0, | 
2. 
eg EP AAP SHO.) ee 
dx? 


Fiir ein Stabelement, dessen simtliche Fasern sich ausschlieBlich im Belastungsbereich befinden, 
wird aus (9) 
EC Ee —_ AP 0 
d?w* 


PAF ere es 


dx2 


(9* a, b) 


wenn zum Unterschied gegeniiber dem durch (9) gekennzeichneten allgemeinen Fall die Dehnung 
der Stabachse mit e* und die Durchbiegung mit w* bezeichnet wird. Die Differentialgleichungen 
(9) und (9*) gelten auch noch fir beliebig veranderlichen Querschnitt. 


4, Lésung der Differentialgleichungen. a) Die Engesser-Karmansche Lésung. Bei Be- 
trachtung der Differentialgleichungen (9) fallt sofort auf, daB es sich — wenigstens zunachst — gar- 
nicht um ein Eigenwertproblem handelt, sondern wegen der Exzentrizitat e, die das Stérungs- 
glied APe verursacht, um ein inhomogenes Knickbiegungsproblem. Es sei jedoch betont, daB 
diese Tatsache zwar durch die gewahlte Darstellung, bei der das Biegemoment auf den Schwer- 
punkt S, der ideellen Flache bezogen ist, besonders deutlich wird, aber nicht etwa durch irgend- 
welche neuartigen Annahmen zustande kommt; die bisherige Rechnung ist vielmehr nichts 
anderes als eine Wiederholung der iiblichen Ansatze von Engesser-Karmén in einer fiir das 
Weitere geeigneten Form. 

Die altbekannte Lésung des Knickproblems ergibt sich nun, wie folgt. Die Gleichungen (9*) 
zeigen zundchst, daB keine Lésung existiert, bei der sich alle Querschnitte im Belastungsbereich 
befinden. Aus der homogenen Gleichung (9* b) entnimmt man namlich, daB es sich um ein 
Eigenwertproblem mit bestimmten Werten P = P, handelt, daf also eine Laststeigerung tiber 
diese Werte hinaus bei ausgebogenem Stab unméglich ist und infolgedessen in (9* a) die Last- 
steigerung AP = 0 zu setzen ist. Fir einen Stab, der zusatzlich zum Grundzustand nur Biege- 
momente bekommt, ist aber die Bedingung ¢, < 0 fiir simtliche Punkte aller Querschnitte nicht 
erfiillbar. Es muB also das Vorhandensein von Entlastungsbereichen angenommen werden. Zu 
ihrer Berechnung seien vorerst die Randbedingungen naher betrachtet. 

Die Voraussetzung eines festen Gelenklagers verlangt 


fir «= 0 ond x=1: w=0, M25 =0. (10 a, b) 
Die Randbedingung (10 b) bedeutet, daB die Differentialgleichung (9 b), in der das Momenten- 


gleichgewicht zum Ausdruck kommt, an den Stabenden mit w = 0 erfiillt sein soll. Es mu8 also 


, 7 dw 
fir x =0 und x =I: E’ J +APe=0 (11) 


° dx 
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2 


gelten. Da an den Stabenden die Belastung nur aus einer Druckkraft im Querschnittsschwer- 


F Pa, 
punkt S besteht, kann sich dort lediglich eine konstante Spannungsverteilung 0 = # einstellen, 


; : 7 dew 
wihrend die beim Ausbiegen des Stabes zusatzlich auftretenden Biegemomente —E’ J, de 
verschwinden miissen. (11) zerfallt also in die beiden Forderungen 

17 @& 

fir x=0 und x=/: E'J,73=0, APe=0. (12 a, b) 


Fiir das Folgende ist vor allem (12 b) von Bedeutung. 

Da eine Liésung gesucht wird, bei der die Stabquerschnitte Entlastungsbereiche haben, liegt 
es nahe, an den Auflagern e 0 zu verlangen, so das dann AP = 0 folgt. Hiermit wird aber 
aus (9b) das Engesser-Kdrmdénsche Eigenwertproblem mit den Kigenwerten P,, deren niedrigster 
mit 


peas (13) 


bezeichnet sei. Das zugehorige Tragheitsmoment J, folgt dabei wegen JP = 0 in bekannter 
Weise aus der Bedingung, daB die Verteilung der zusatzlich beim Knicken auftretenden Span- 
nungen auch wirklich nur ein reines Biegemoment liefert. 

b) Vollstandige Lésung. Wesentlich ist nun, daB die Forderung e # 0 fiir die Stabenden 
nicht die einzige ist, die zu einer Lésung des Problems fiihrt. Die Forderung (12b) kann auch 
mit e = 0, AP £0 befriedigt werden, ohne daB damit Verformungen des Stabes im Entlastungs- 
bereich vollig ausgeschlossen sind. Wie solche Lésungen zustande kommen kénnen, mége zu- 
nichst folgende anschauliche Uberlegung zeigen. 

Bei hinreichend kleinen Lasten kommt nur die Gleichgewichtslage des gestreckten Stabes 
in Betracht, wobei sich alle Querschnitte vollstandig im Belastungsbereich befinden. Mit be- 
ginnender Auslenkung stellt sich cin Biegemoment ein, das in der Stabmitte seinen GréBtwert 
hat und nach den Stabenden hin auf Null 
abnimmt. Es ist danach zu erwarten, daB 
sich der durch das Biegemoment erzwungene 
Entlastungsbereich auch zunachst nur in der 


| Fntlastungsbereich 


l-a a , a . 5 . . . 

8 $44 Stabmitte ausbildet und dann erst mit steigen 
Belastungs = Entlastungs- der Belastung nach den Seiten hin ausbreitet. 
abschutt abschnitt 


Die ganze Stablange wird sich also im all- 
Abb. 4. Knickstab mit Entlastungsbereich. gemeinen nach Abb. 4 in einen ,-ntlastungs- 
abschnitt, der die Linge a haben mige, 
und zwei ,,Belastungsabschnitte‘‘ unterteilen. Die Dicke s des Entlastungsbereiches wird in 
Stabmitte ihren GréStwert haben und nach irgend einem zundchst unbekannten Gesetz bis 
zum Wert Null an der Grenze zwischen dem Entlastungs- und Belastungsabschnitt abnehmen. 
An den Auflagern mu e = 0 sein. Da jedoch im Entlastungsabschnitt e 0 ist, wird sich der 
Stab unter Wirkung des Momentes APe sofort durchbiegen, so daB ein Knickbiegungsproblem 
und kein EKigenwertproblem vorliegt. 

Das geschilderte Verhalten des Stabes wird durch die Rechnung voll bestatigt. Es seien zu- 
nichst die Belastungsabschnitte untersucht, wobei aus Symmetriegriinden die Betrachtung 
eines Abschnitts geniigt. Hierfiir gelten die Differentialgleichungen (9*). Setzen wir konstanten 
Querschnitt voraus, so ist die allgemeine Lésung von (9*) mit den Integrationskonstanten A 


und B 


w* = Asin py st Beos|/ ary e. 


Das Verschwinden der Durchbiegung bei x = 0 verlangt B = 0, so daB 


: it 4 
w* — A sin lea a (14) 
iibrig bleibt. 

Zur Ermittlung der Durchbiegung im Entlastungsabschnitt wird noch folgende geometrische 
Beziehung gebraucht, die sich aus Gleichung (4) herleiten 1aBt. Nach Definition wird fiir z — c —s 
die Dehnung e¢, = 0; es ist also 

d?w 
s) de’ 
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Setzt man fiir «, den Wert von (9a) ein, so wird 


ar = Ae I 
dx? EF’ Fec—e—s~ (15) 


Durch (15) und (9b) wird das im folgenden zu lésende Problem heschrieben. Differenziert man 
(9b) zweimal nach x und eliminiert mit Hilfe von (15) d?w/dx?, so erhalt man eine Differential- 
gleichung vierter Ordnung fiir s. Sie ist nichtlinear, da s in (15) im Nenner vorkommt und e, 
Fund J,irgendwelche durch die Querschnittsform gegebenen Funktionen von s sind. Trotzdem 
ist es gliicklicherweise méglich, auch fiir beliebige Querschnittsgestalt die Lésung auf folgende 
Weise zu gewinnen. 


Aus (9b) folgt unter Benutzung von (15) 


oie de ed 
AP Ee toe Ss ees 


Werden die dimensionslosen GréBen 


s e . oDe P w 
ate = Vw APG a 
eingefiihrt, so wird 
12 
jas Se = Le ears Es 17 
i l—y—yp i es) 
und aus (15) 
ft APE 1 
de Rc il ayy 
Multipliziert man die linke Seite dieser Gleichung mit dw/dx, die rechte nach (16d) mit 
dw AP p dp 
di Pits 
und integriert, so folgt 
Muga to CeO err 1 . 
aa) per | mizgay® rast 


Zur Festlegung der Integrationskonstanten C dient die Bedingung, daB in Stabmitte dw/dx = 0 
sein mu}. Der an dieser Stelle vorhandene Wert von @p sei Yn; bei y = 0 sei mp = @y. Es wird 


dann 
Pm P 
yi sc ec 1 vale 
2 tae) scent et Shag eee ¥ Ke =o oe 


Po Po 


Setzt man zur Abkirzung 


y Pm 
F 1 F I 
ge ca == eee ee ir, 18 a,b 
2 |/: | ees P» Pr \/:f: 1—7n—y ? (18 a, b) 
Po Po 
so folgt 
dw AP? ne abe 19 
eV oe | on We (19) 
Differenziert man (16d) nach x, so ergibt sich 
ARP. 6 
dx 


und mit (19) nach Integration 
~ 
l—a E’F il 
= ——— -dQ. 20 
CET +055 = Y ey 
Po 


Die Integrationskonstante ist dabei gleich 5 (l — a) gesetzt, da x diesen Wert bei p = Y, d-h. 
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an der Ubergangsstelle zwischen Belastungs- und Entlastungsbereich annehmen mu. Das in 
(20) vorkommende Integral sei zur Abkitrzung mit @,, O bezeichnet und die neue Veranderliche 
j= ate sin. oe (21) 
eingefithrt, wobei } als Hauptwert des arc sin definiert sei. Beachtet man, daB®B nach (18a) und 
(21) & bei p= verschwindet, so wird 
? o 


aes Ate sl elle 
Dn OF fone | os cos re 
ig m . 
ae 1 
dp 1 d® Fei—n—y 
dp Pm cos 3 dp 2, sin } cos > 
d | 
= | FU —n—y)sin Odd (22) 
0 
und damit aus (20) 
l—a i Eiko E 
—— —— : 23 
x 9 + c | P P,,, 0 ( ) 
Aus der Bedingung, daB bei x = if , d.h. in Stabmitte, py =, ® = ®, und nach (21) = a 
ist, folgt mit 
n/2 
On = | zl 7 20) cin) aD (24) 
0 
aus (23) 
a (yu te 
Ga yee DAO 23 (25) 


Durch die Einfithrung von @ ergibt sich zwanglos, daB das Integral ®,, O,,, das in der Dar- 
stellungsform der Gleichung (20) wegen der Integration bis ® = @,, uneigentlich ist, einen end- 
lichen Wert hat. Die Benutzung von # erweist sich auBerdem fiir numerische Rechnungen als sehr 
zweckmafig. Durch Aufstellung der Gleichungen (19) und (20) ist die Integration der Diffe- 
rentialgleichungen (9b) und (15) auf Quadraturen zuriickgefiihrt und damit im Prinzip erledigt. 

Als nichstes sind die Ubergangshedingungen zwischen Entlastungs- und Belastungsabschnitt 
zu betrachten. An dieser Stelle miissen die Durchbiegungen und ihre Ableitungen iibereinstimmen, 
d.h. es mu gelten 

l—a dw* dw 


it? Ge a OP ee pss es (26 a, b) 


dw* d , 
w* und — folgen aus (14); w und 7 ergeben sich aus (16d), (17) und (19). Dabei ist zu setzen 


y=0, 2=0, F=Fogl su) ee eo] ot eee 


Man erhalt dann 


A sin Ves 3° = eit, AV ry 08 ee aed Bea jpelee ee 


Nach Elimination von A folgt 


Mit der abkiirzenden Bezeichnung 
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ergibt sich unter Benutzung von (25) 


1 T7eP - 1 
“8 ®,, tg G | re Sara a ®@,, On) == i . (29) 


5. Diskussion der Lésung am Beispiel des Rechteckquerschnittes. Eine Diskussion der in den 
obigen Formeln enthaltenen Lésung erfolgt am zweckmafigsten an Hand eines Beispiels. Hierzu 
sei der Stab mit Rechteckquerschnitt gewahlt. Fiir diese Querschnittsform ergeben sich nach 
elementarer Rechnung die dimensionslosen Ausgangsgréfen zu 


n= n— 3 (2 re! 
2 2+(m—1)y’ 


Fe n—] 
oe 
LF (30 a—d) 
== eh — ly —3y +3) —37 [2 -(n — 1) y] + 2}, 
= 
os 


Im folgenden ist c die halbe Hihe des Rechtecks. Es ist bemerkenswert, da® diese Ausgangs- 
gréBen und damit auch, wie sich zeigen wird, die Ergebnisse der weiteren Rechnung unabhangig 
vom Seitenverhaltnis des Rechtecks sind. 


Die Lésung der Aufgabe, den Verformungszustand des Stabes in Abhangigkeit von der Be- 
lastung darzustellen, wird durch die Gleichungen (16d), (23) und (29) geliefert. Gibt man y als 
unabhangige Verdnderliche vor, so erhalt man aus (16d) nach Berechnung von @ die Durch- 
biegung w und aus (23) nach Ermittlung von @,, und @ die zugehérige Koordinate x, wenn die 
Parameter y,,, P und AP = P — P, ge- Ke 
geben sind. Zwischen diesenletztgenannten 74\4 
GréBen besteht die transzendente Uber- 
gangsbedingung (29), so dafi damit noch 
eine Gréfe, etwa y,, festgelegt ist. Zwei 
Parameter, z.B. P und P,, sind aber 
noch frei wahlbar. Stellt man die Be- 28 
lastung P in Abhangigkeit von der Ver- “s stabiles Gleiohgewioht tir n=4 
formung in Stabmitte dar, so erhalt man 0 /abiles Gleichgemht fir r= 
also eine ganze Kurvenschar mit P, als 
Parameter. Diese Rechnung laBt sichleider 


j > $ : 0 2 ay 06 08 0 72 Ww 
fiir den Rechteckquerschnitt nicht in ge- 4), 5 P/Px in Abhingigkeit von w/e mit P, /Px als Parameter 


by 


oss 


schlossener Form durchfiihren. Da es sich fas ak A eokteckquorohaite. 
jedoch um einfache Quadraturen handelt, 
kénnen diese ohne Schwierigkeiten — und dasselbe gilt auch fiir andere Querschnittsformen — 


numerisch durchgefiihrt werden. Das Ergebnis ist in den Abb.5 bis 9 dargestellt. 

Abb.5 zeigt zunachst fiir den willkirlich gewahlten Wert n = 4 die Kurvenschar P/Px in 
Abhangigkeit von w,,/c mit P,/Px als Parameter, wobei w,, die Durchbiegung in Stabmitte ist. 
Die Kurven beginnen bei w = 0 mit den verschiedenen Werten P = P) und erfiillen den Bereich 
- zwischen der stark ausgezogenen untersten Kurve und der horizontalen Geraden. Diese Gerade 
ist die Engesser-Kdrmdnsche Lisung des Problems, da hierbei 4A P = 0 ist, so dafs Gleichung (9b) 
mit der Lésung (13) gilt. Beim Rechteckquerschnitt ist der bekannte Engesser-Kdrmdnsche 
Higenwert 

ip go AT 

Pe (1+Jn)? 
Die horizontale Gerade ist zugleich fiir alle Kurven Asymptote, die bei unendlich groBen Durch- 
biegungen erreicht wird. Im Grenzfall ist die Lange a des Entlastungsabschnittes gleich 1. Mit 
w = co muB namlich nach (16d) auch y = oo sein, was fiir 1 — 1 —p = 0 eintritt, so daB nach 
(18) ®,, = co und nach (27) a = 1 folgt. Wahrend sich bei der Engesser-Kdrmdnschen Lisung 
der Entlastungsbereich sofort beim Beginn des Knickens ither die ganze Lange erstreckt, fangen 
die tibrigen Kurven bei w,, = 0 mit einem Wert a an, der kleiner als / und um so kiirzer ist, je 
kleiner Py ist. 
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Besonders wichtig ist die stark ausgezogene untere Begrenzungskurve des Bereiches. Sie 
liefert nimlich die stabilen Gleichgewichtslagen des Stabes, da sie die niedrigsten Lasten P kenn- 
zeichnet, die eine bestimmte Durchbiegung zu erzeugen vermégen. Alle itbrigen Kurven, also 
auch die Engesser-Kdrmdnsche Gerade, gehéren zu labilen Gleichgewichtszustaénden und sind 
damit praktisch von geringem Interesse. Die unterste Kurve beginnt mit 4 = Oibert PE=eP x: 
Fir a = 0 wird namlich —~,, = Y und nach (18) @,, = 0, so da in der Tat aus (27) 


0. tei ie =1,, Pa Px 
DE Neve 

folgt. Diese wichtige Aussage, die bereits von Shanley auf Grund anschaulicher Betrachtungen 
als allgemein giiltige Tatsache ausgesprochen wurde, bedeutet also, dafB die Ausbiegungen des 
Stabes gerade bei der Last anfangen, bei der auch »gewohnliches Eulerknicken“ eintreten wirde, 
wenn fiir den ganzen Stab, also auch im Entlastungsabschnitt, nur der Elastizitatsmodul E’ 
gelten wiirde. 

Eine sehr einfache Beziehung gilt fiir die Tangentenneigung der niedrigsten Kraft-Ver- 
formungskurve von Abb. 5. Nach (16d) und (17) wird namlich wegen a = 0 


Pet 
rs PK il 
igen ewe (De = sae a (31) 


Diese Beziehung ! gilt ebenso wie alle bisherigen Aussagen iiber das grundsatzliche Verhalten 
der Kurven zu Abb.5 fiir beliebige Querschnittsform. Der Kurvenverlauf fiir gréBere AP und 
auch die Tangenten der Kurven, die zu den labilen Gleichgewichtszustanden gehéren, sind je- 
doch von der Querschnittsform abhangig. 

Bei gewobnlichen Verzweigungsproblemen hat die Kraft-Verformungskurve im Abzweig- 
punkt eine horizontale Tangente, d.h. bei ein und derselben Last sind verschiedene benach- 


barte Gleichgewichtslagen méglich. Man 


14, — spricht in diesem Fall bekanntlich vom 
je L | fe ea “<i indifferenten Gleichgewicht, das in der 
Regel beim Ubergang von stabilen zu 

j) |i PY ot I a rn labilen Gleichgewichtszustanden durch- 
oe laufen wird. Das vorliegende Problem 

16 Re ist durch die nach (31) von Null ver- 
schiedene Tangente als Ausnahme ge- 

Welne a ee SS Se kennzeichnet. Die Gleichgewichtslage 
_; nee | des geraden Stabes wird bei P = Px 
unvermittelt labil, indifferentes Gleich- 

10 aie if | gewicht tritt nicht auf. In ,,energetischer 
ie Ausdrucksweise“* bedeutet dies, daB die 

185 aE tL 7 ie ss 7 7 ice zweite Variation der potentiellen Ener- 


gie des geraden Stabes unterhalb des 
Abzweigpunktes positiv definit ist, im 
Abzweigpunkt jedoch sprunghaft nega- 
tive Werte annehmen kann. Dieser Sprung wird durch die Unstetigkeit im Elastizitatsgesetz 
verursacht und ist fiir derartige Falle typisch ?. 

Es muf} noch erwahnt werden, da au er der in Abb.5 dargestellten Kurvenschar auch noch 
weitere solcher Scharen existieren, die zu den Knickmiglichkeiten des Stabes mit mehreren Halb- 
wellen gehéren und sich ergeben, wenn man in (29) statt Px die entsprechenden héheren Eigen- 
werte m7? E’ J/P? (m = 2, 3,...) einsetzt. Diese Kurven liefern jedoch samtlich labile Gleich- 
gewichtszustande und sind daher praktisch ohne Belang. 

Zur weiteren Kennzeichnung der Zusammenhinge sind noch fiir verschiedene Werte des Ver- 
haltnisses n in Abb.6 die niedrigsten Kraft-Verformungskurven, in Abb.7 die Abhangigkeit 
zwischen Belastung und Lange des Entlastungsabschnitts und in Abb.8 die Abhangigkeit zwi- 
schen Belastung und gréBter Dicke des Entlastungsbereiches dargestellt. An der letztgenannten 


Abb, 6. P/ Pr in Abhangigkeit von Wm/¢ mit n als Parameter fir stabile 
Gleichgewichtszustinde. Rechteckquerschnitt. 


__* Die Tangentenneigung nach (31) ist bereits in der von F. Schleicher, Bau-Ing. 26 (1951), S. 140, ent- 
wickelten Naiherungslésung des Problems enthalten. 


2 Vel. A. Pfliiger, Stabilitatsprobleme der Elastostatik, S. 66. Berlin 1951. 
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Kurvenschar ist vor allem bemerkenswert, da der Entlastungsbereich um so grofer ist, je kleiner 
nist. Eine anschauliche Zusammenfassung der wichtigsten Aussagen der Diagramme 5 bis 8 
zeigt Abb.9, aus der fiir n = 4 die Zunahme und Ausbreitung des Entlastungsbereiches mit fort- 
schreitender Verformung hervorgeht. Dazu sei noch einmal betont, daB® die Rechnung wegen der 
grundlegenden Annahmen iiber das Spannungs-Dehnungsdiagramm nur fiir hinreichend kleine 
Verformungen das wirkliche Verhalten des Stabes wiedergibt. 


a 
l 
= 
: TE er 
Abb. 7. P/PK in Abhangigkeit von a/l mit n als Abb. 8. P/Px in Abhangigkeit von y,,, mit n als 


Parameter. Rechteckquerschnitt. Parameter. Rechteckquerschnitt. 


Nach der alten Engesser-Kdérmdnschen Theorie galt die Tragfahigkeit des Stabes erst mit 


Erreichung des Wertes Px als erschépft. Nach der vollstandigen Lésung wird dieser Zustand 
schon friither erreicht, wenn auch nicht sofort bei Uberschreitung der Last Px, da ja der Stab 
auch fiir Werte P > Px noch Lasten tragen kann. Fir eine Uberschreitung von Px um 10% 
ergibt sich namlich bei n = 4 eine Durchbiegung 
von etwa 2,5°% der Querschnittshéhe fiir den 
Rechteckquerschnitt (nach Abb. 6) und von 
etwa 7% fiir den Zweipunktquerschnitt (nach 
Abb.10), was in der Regel noch zulassig sein 
wird. Auf die Erérterung weiterer Einzelheiten 
dieser und ahnlicher Fragen praktischer Bedeu- 
tung sei jedoch hier verzichtet, da dies bereits 
in den Arbeiten der in diesem Aufsatz genannten 
Autoren geschehen ist. Es sei nur die Ansicht 
des Verfassers betont, daB keine Veranlassung 
zu einer Anderung der fiir die Bemessung von 
Knickstaben festgelegten w- Werte besteht. Denn 
diese sind ja in Ubereinstimmung mit Versuchs- 


ETT 


PsP, P=10rn, Pa1BR Paisp Pate P=16R 


ergebnissen gewahlt, in denen —— WV Da: auch Abb. 9. Verformung und Entlastungsbereich in Abhangigkeit 
: von der Belastung fiir stabile Gleichgewichtszustande. 
bisher unerkannt — der Shanley-Effekt ent- Hechreck qusrschnren = 4! 


halten ist. Der wesentliche Wert der neuen 
Theorie diirfte vielmehr in der Erkenntnis bestehen, da® das in Versuchen immer wieder fest- 


gestellte von einem Eigenwertproblem verschiedene Verhalten des Knickstabes nicht nur durch 
Vorverformungen und unbeabsichtigte Anfangsexzentrizitaten verursacht sein muB, sondern 
auch durch plastische Verformungen zustande kommen kann. 


6. Zweipunktquerschnitt. Ein I-Querschnit kann bei kraftigen Flanschen und schwachem 
Steg naherungsweise durch einen Zweipunktquerschnitt ersetzt werden, bei dem die Querschnitts- 
flache in zwei Punkten konzentriert gedacht ist. Dieser Querschnitt sei hier noch kurz behandelt , 
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da sich in diesem Fall die erforderlichen Integrationen geschlossen durchfithren lassen '. Wegen 
der z.T. zu einfachen Zusammenhinge ist allerdings der Zweipunktquerschnitt nicht so gut als 
Musterbeispiel gecignet wie der im vorigen Abschnitt behandelte Rechteckquerschnitt. 


Fiir die AusgangsgréBen gilt 


i! Le pee A ehee ay 
Va Fone a Cae 
c ist auch hier die halbe Querschnittshéhe. 


Fir die durchzufithrenden Integrationen wird am zweckmaBigsten y als unabhangige Ver- 


anderliche benutzt. Aus (18) folgt mit dp/dp nach (17) 


a 3 Sa 
ia dy Bim 1 iu 
=o 2 —"l, eas — ¢ 
? Fe a (l—7n—y)?’ 2 VE . a =i =k ce 
0 


Fiir das in (20) vorkommende Integral @,, O ergibt sich 


y | 

| 
W5@ = || = z Me dy | 
| 03, — a ap 
0 


und mit (33) und (17) 


yp | 


l 1 pe 
mi — — —— -— SS eee bre ; ] 
P, re {EF oe | f 1 \Q il \2 (1 — 7 — wy)? dy 
Bcd ier erie a a 


0 


Unter Benutzung von (32) wird 
y 
“ (2 = — Ym\? 


17a dy 


oD) (6) = ea = = —= _2n = | : < —- 
beat (ee Bese) rat, oneal nee 1—n—y 
| eae 
( 


y/ 2n A oun hy — ep, l—7y—yp Ea 
®,,, oe = | Sie ( arc § == Migs ay ° era VER == yt t= 
ee Ags are sin ———— “aes eas: cos cs D.,.0),).  Aaag 


Aus (33) folgt mit (32) 


2 be Dion hee ae / ms 2 
2 = Vorrerp | pee) (ral Veer eee ie 


und mit (34) 


ey Tut Lge XG 
Pn = n+1 '8 (| coe Pr On) ; 
Aus der Ubergangsbedingung (29) wird dann 
of) mae a 1/P n+ 1 
ts(|/ Oy @,, O,.)-t6 ( Q | P Pte P,, On| = os q (35) 


Fiir die gréfte Durchbiegung w,, in Stabmitte ergibt sich nach (16d) und (17) 


Wm A 12 i? 
= n\. 
¢c P \l— aT 


2 n 1 nw — I 

fei ae 3 ib 
ae (| tala On) real 

ae nTD 


' Eine Naherungsrechnung fiir den Zweipunktquerschnitt wurde von U, Miillersdorf, Bau-Ing. 27 (1952) 


5.57 veréffentlicht. Diese Nah h i ii i 
pakiairigcn Beekeencn aherung hat den Vorteil, daB genaue gekriimmte Spannungs-Dehnungsdiagramm 


und mit (34) und (32) 
ier ty 


reine D 


(36) 
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a Die Formeln (35) und (36) stellen die Lésung des Problems dar. Bei gegebenem P, und AP 
olgt aus (35) zunichst @,, @,, und dann aus (36) w,,. Nach (25) und (28) ist mit i? = 1 


ase gery PF 
P,,, @.. idee’ 1 9 | Px 2 


Damit nimmt (35) die Form 


Fo ee ea ‘ Picay a a a ong 

s(| 2 ee | Peepers | jobs a he Vee 

an, aus der sich besonders bequem die Grenzfille von P ablesen lassen, die sich fiir a — 0 und 
a =I ergeben. Im ersten Fall wird 


az 1/P 
tes \po = ParPz 
und im zweiten 
a w= 1 ve i 2n = 
t {= — —s - —— 
I s(| ee Al Pel tne? Dante < jane eae 


also gleich dem Engesser-Kdrmdnschen Wert, da in der Tat beim Zweipunktquerschnitt 


aaa ay 

; n+] 

ist. 

Zur Weiteren Veranschaulichung der Zusammenhange sind fiir verschiedene Parameter n 
in Abb. 10 die zu den stabilen Gleichgewichtslagen gehérenden Kraft-Verformungskurven und 
in Abb.11 die Abhangigkeit zwischen P/Px und a/l aufgetragen. Grundsitzlich Neues gegen- 
tiber den Verhiltnissen beim Rechteckquerschnitt AP 


ergibt sich dabei nicht. py te 
© 1765 
5 
78\~ 
eee ee oe ee = 
oat 2S SS ee See | 
| n=6 
14 + 
i) Es ee eae ae — —+ ai 
| | n=2 
12) 
10 
| Yn a 
0B L | set ke ae L 
0 02 GY 06 48 UL lé bie ad. G2 G4 46 G8 10 
Abb. 10. P/Px in Abhangigkeit von m/e mit n als Parameter. Abb. 11. P/Px in Abhingigkeit von a/l mit n als 
Zweipunktquerschnitt. Parameter. Zweipunktquerschnitt. 


7. Zusammenfassung. Fiir einen beiderseits gelenkig gelagerten Knickstab mit beliebigem, 
aber iiber die Stablange konstantem Querschnitt wird das ,,Eulerknicken“ (nicht das Biegedreh- 
knicken) im plastischem Bereich den Gedankengangen von Shanley folgend untersucht. Zur Er- 
fassung des plastischen Verhaltens des Stabes werden die Annahmen der Engesser-Kdrmdnschen | 
Theorie benutzt. Die Rechnung fithrt mathematisch auf eine nichtlineare Differentialgleichung 
fiir die Breite des sich beim Knicken ausbildenden ,,Entlastungsbereiches*. Mechanisch gesehen 
ergibt sich ein Knickbiegungsproblem, da im Entlastungsbereich eine Schwerpunktverschiebung 
des ,,ideellen‘‘ wirksamen Querschnitts gegeniiber der Wirkungslinie der Last und damit zu- 
satzliche Biegemomente entstehen. Bedingt durch die Unstetigkeit im Elastizitatsgesetz geht die 
urspriinglich stabile Gleichgewichtsform des geraden Stabes in die Gleichgewichtslage des aus- 
gebogenen Stabes iiber, ohne dafs dabei ein Punkt indifferenten Gleichgewichts auftritt. Die 
Engesser-Kdrmdnsche Lésung erscheint als eine Sonderlésung des Problems, die jedoch zu labilen 
Gleichgewichtszustainden gehért. Die Zusammenhange werden zahlenmafig am Rechteckquer- 
schnitt und am Zweipunktquerschnitt untersucht. 


(Eingegangen am 20. Marz 1952.) 
Anschrift des Verfassers: Prof. Dr.-Ing. habil. A. Pfliiger, Hannover, DieterichsstraBe 36. 
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Eigenschwingungen idealer Fliissigkeiten in Rohrleitungen 
mit verschiedenen Endquerschnitten. 


Von Hsien-Chih Liu. 


1. Einleitung. Die analytische Behandlung von Schwingungen einer Fliissigkeitssdule in Rohr- 
leitungen mit veranderlichem Querschnitt fiihrt auf die Zuhilfenahme einer vereinfachenden 
RechengréBe, der sogenannten reduzierten Linge. Dies ist seit Jahizehnten bekannt. Unter der 
reduzierten Lange versteht man die Linge derjenigen zylindrischen Fliissigkeitssaule, die der 
betrachteten Fliissigkeitssdule mit verdnderlichem Querschnitt in ihrem physikalischen Ver- 
halten aquivalent ist. Wahrend der Schwingung einer beliebig geformten Fliissigkeitssaule andert 
sich mit ihren Ortskoordinaten im allgemeinen auch die jeweilige GréBe der entsprechenden redu- 
zierten Lange. Bei sehr groBen reduzierten Laingen und sehr kleiner Amplitude kann diese 
Anderung vernachlassigt werden. Trifft ersteres nicht zu, so darf die Anderung der reduzierten 
Lange unter Umstinden nicht unberiicksichtigt bleiben; d. h. die Gri Be der Anderung der redu- 
zierten Lange kann nicht mehr als vernachlassigbar klein gegen die urspriinglicne reduzierte Lange 
betrachtet werden. Das Ziel der vorliegenden Arbeit! ist, den KinfluB® zu untersuchen, den diese 
Anderung der reduzierten Linge auf die Eigenschwingung hat. Der Einfachheit halber wird hier 
nur derjenige Fall behandelt, bei dem die Enden der Fliissigkeitssaule wahrend der Bewegung 
zwar beliebig groBe aber jeweils konstante Rohrabschnitte bestreichen. 


2. Bezeichnungen und Voraussetzungen. Das Schwingungssystem miége durch die in Abb. 1 
gegebene Anordnung schematisch dargestellt sein. Die Rohrleitung ist im Ruhezustand von 
Stelle 1 bis 2 mit Fliissigkeit gefii'lt. Die Feder c, der Kolben F’, und die dariiber stehende Fliissig- 
keitssaéule bilden also das Schwingungs- 
system. 

Es sei s,(t) der Abstand der an die 
Fliissigkeit grenzenden Kolbenflache von 
einer festen Bezugsebene 0—O, ferner 5, 
die Kolbengeschwindigkeit, m die Masse 
des Kolbens, ¢ die Federkonstante, F, der 
Kolbenquerschnitt bzw. Querschnitt des 
vom Kolben bestrichenen zylindrischen 
Rohrleitungsteiles, s(#) der Abstand des 
beliebigen Rohrquerschnittes F von der im 
Raum festen Bezugsebene 0—0, F, der 
Endquerschnitt der vom Kolben hin- und 
hergeschobenen Fliissigkeitssaule, s.(t) ihr 
Abstand von der festen Bezugsebene, 
P2 == Po der Atmospharendruck am freien Ende der Fliissigkeitssaule, p, der vom Kolben auf 
die Fliissigkeitssaule ausgettbte Druck und p, F, die entsprechende Kolbenkraft. Es soll keine 
. Kavitation auftreten. F = F(s) ist eine von der Zeit unabhangige, konstante Funktion. 

AuBer dem atmospharischen Druck py wirkt auf die untere Kolbenseite die Federkraft cs,; 
letztere verschwindet an der Stelle s; = 0. 

Die Fliissigkeit sei reibungsfrei und inkompressibel. Die Reibung zwischen Kolben und Rohr- 
wand wird vernachlassigt. 


i] \ C ‘% 


Abb. 1. Festlegung der Koordinaten fiir das Schwingungssystem. 


Es wird weiterhin \F <K $2. — Ss; vorausgesetzt und die seitlichen Geschwindigkeitskomponenten 
kénnen vernachlassigt werden. 

Endlich sei v die Geschwindigkeit der Flissigkeit an der beliebigen Stelle s, ferner « der Winkel 
zwischen der Vertikalen und der Tangente an die Mittellinie der Fliissigkeitssaule an einer belie- 
higen Stelle s und h der Héhenunterschied zwischen den Enden der Fliissigkeitssaule. 


' Diese Arbeit ist von mir gr6Btenteils fertiggestellt, wahrend ich (1943) als standiger Assistent am Lehr- 
stuhl und Institut fiir Mechanik der Technischen Hochschule zu Berlin tatig war. An dieser Stelle méchte 
ich nicht versiumen, Herrn Prof. Kucharski meinen Dank auszusprechen fir einige Anregungen, die zur 
Vollendung dieser Arbeit beitrugen. 
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3. Ableitung der Bewegungsgleichung aus der Bernoullischen Gleichung fiir instationire Stré- 


mung. Die Bewegungsgleichung fiir den Kolben ergibt sich unmittelbar aus der Newtonschen 
Grundgleichung 


Ss 
m 7 (Pi — Po) Fi, + ¢5, == 0. . (1) 


Die Druckdifferenz p, — py = py — Pp» lat sich ebenfalls sofort nach der bekannten verall- 
gemeinerten Bernoullischen Gleichung hinschreiben: 


Pi—Po=yh + © vs) + oh SE mit La | pds. (2) 


Hierin bedeuten @ die spezifische Masse und y das spezifische Gewicht der Flissigkeit. 


Setzt man 


z fede" gst 
_F, = 3 Aig 
ee ye? ey 1 (3) 
dann wird ut aft 
f= ras 
2 A sp wl a | ds, \? ea abe 
OB Ol (73 == ) ‘i) ; S2min—-1*-S2ma 
a b€/ 84=Symin\ O6/ S4=Symax 
So Ss | 
mit | ee if nds , (4) ~ om 
a 59=Spq b/ 5720 
53=0 
und hiermit lautet die Bewegungs- ; 
: onderfall des Rohres, 1=konst. bei Stelle =F,=kon 
gleichung (1) Abb. 2. Sonderfall des Roh Be ai sts t. bei Stelle 1, F st. 
(m +o | nds) wet +S hGi—y (Gy testy hh=o, (5) 


sy 
Hierin sind 4, 75, h urd s, bekannte Funktionen von s,, das sich mit der Zeit indert. Bei der redu- 
2 2 1 


zierten Lange L andern sich mit s, die Grenzen des Integrales; also ist auch L eine bekannte Funk- 
tion von s,. 


4, Ein Sonderfall der Bewegungsgleichung. Bei den Stellen 1 und 2 in Abb. 2 sollen geniigend 
lange Sticke mit F', = konst. und F,, = konst. vorhanden sein. Diese beiden Stiicke mégen ver- 
schiedene Neigungen besitzen, also a, ~ dy. Jetzt gilt 


Go L= al Pr Se honey (6) 


2 


Die reduzierte Lange schreibt man zweckmabig 


L=f'n ds = fy ds — fy ds +I'y ds ; (7) 
hierin ist i ; E 
= ic ds = konstz, fin doz si fin As ay 
also : 


L=1L,+ (fz—1)s,, h = ho + (Hp cos & — COs a) 8; 


Damit lautet die Bewegungsgleichung 


@ aks 
(m + @ F, Ly + 0 Fy (43 —1) 5) My | 2 Fi (te wa) 4 | 


Wenn man 


2 Fyh 
y F, (nq cos & — cos a) + €= 4, Oe Orgs A ar 


oF, (n3 — 1) = 2a, $— 31 
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in (8) eingesetzt, so erhalt man schlieBlich die Bewegungsgleichung in der Form 


(m, + 206,) 3 4 “(ae) Led = 0. (9) 


5. Integration der Differentialgleichung fiir den Fall F', = konst., F, = konst. a) Die not- 
wendige Einschrankung der Schwingungsamplitude. Der Wert 77, in (9) ist je nach der 
Rohrgestaltung sehr mannigfaltig, das Gleiche gilt dementsprechend auch fiir die Lésung dieser 
Differentialgleichung. Wegen der groBen mathematischen Schwierigkeit beschranken wir uns 
hier nur auf die Untersuchung fiir den Sonderfall F, = konst. und F’, = konst. Nach dieser Ejin- 
schrinkung miissen die beiden Enden der Fliissigkeitssdule sich stets in Rohrstiicken von kon- 
stantem Querschnitte bewegen. Das Rohr wird durch die Ziffern J, If, II], und IV in drei Strecken 
unterteilt. Es wird also 

IV UI 


Te : F. 
Lye= | 7 dst See | ra ds + Wo So max + (10) 
t II 
Nach der oben erwahnten Einschrankung sollen die Enden der Fliissigkeitssaule sich wabrend der 
Bewegung niemals innerhalb II—III befinden. Deshalb muf die Anfangsauslenkung (¢ = 0) sein: 
$19 See (11) 


b) Integrationder Bewegungsgleichung. Man kann feststellen, da die Multiplikation 
von (9) mit é, 


d /m, ; : c 
a pelt Cott FAl= 0 
und nach der ersten Integration 
™m, ¢ e s 
H+ og s+ Pas A (12) 


liefern, wobei A = c, &7, durch die Anfangsbedingung t = 0, &; = &,) und é, = 0 bestimmt ist. 
Gleichung (12) laBt sich schreiben 


ate fp Bete (12a) 
| 1 


Gleichung (12a) besagt, da®B die Maximalausschlage auf beiden Seiten der Gleichgewichtslage 
gleich sind. Auffallig ist das Auftreten von €;im Nenner unter der Wurzel. Daraus kann man 
schlieBen, daB die absoluten Werte der Geschwindigkeit bezogen auf gleiche Abstande zu beiden 
Seiten der Gleichgewichtslage nicht gleich sind. Wenn man é, gegen é, auftragt, so erhalt man 
ellipsenartige Kurven (Abb. 4), deren Maximalwerte jedoch nicht auf der Achse &, = 0 liegen. 
Nach welcher Seite und wie weit sich diese Maximalwerte gegen die Mittellage verschieben, wird 
durch den Wert 7, bestimmt. Der Wert 


fee ee les (13) 


liefert den Maxiamalwert von &,. Das positive Zeichen vor der Wurzel gilt fiir den Fall a> 0, das 
negative Zeichen fiir den Fall a<0. Man kann noch zeigen, da 


Ghar a4 


ist, wobei das negative Zeichen fiir a<0 und das positive Zeichen fiir a<0 zu setzen sind. Ferner 
kann man schreiben 


Paes time t aft 
oT (aera eye tek ee 


Die Beschleunigung ist also nicht symmetrisch zur Gleichgewichtslage, (Abb. 3 und Abb. 5), d. h. 
sie hat bei gleichen Ausschlagen zu beiden Seiten der Nullage verschiedene Werte. Im Gegensatz 
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zu den einfachharmonischen Schwingungen ist hier die Beschleunigung bei &, = 0 nicht mehr 
Null, sondern 


al Gisth 
dt eg rer Meera 16 
dt” ]s 9 pina (16) 
Fiir die maximale Be schleunigung hat man 

ae Cy Emax 

( ae), = m, + 2a Gnas ‘ (17) 


Da die rechte Seite von Gleichung (17) keine symmetrische Funktion von Eq, in bezug auf die 
Gleichgewichtslage ist, liefert sie bei gleichen positiven und negativen Ausschlagen zu beiden 
Seiten der Gleichgewichtslage ungleiche absolute Werte. Auf der Seite, auf der die reduzierte 
Lange mit zunehmendem Ausschlag anwachst, ist der absolute Betrag der Beschleunigung gréBen- 
mafig tiberwiegend. Dies gilt auch umgekehrt. 


ys * : 
7000 ade lS hg 16 
0 ; Tae ) 7000\- 100 = 4 
| : 4 
800;- 40} 30 = ‘re 800|- 80), — fs =i 
| é, 
600|- 30} 45 = 
oe I a ie 
400 20; 40 cn ai bf : 
f } i Y00\- 40 ‘ina ; 4 
F é 
E 73 
200; 10\- 5 ay r] zob 20 = 
Diao 0) 0 1 £ bs br bn y 
r oY Noosa Mae paz oe o- Sa \=ag Lo aS 7 eT a 
-200\- -10'--0.5 ee Sa -200+- ~20 
bo 
-400\}- -20|--40 = t -400\- -40 
E 
-600"— -30'—-4.5 -600'— -60 
Abb. 3. Wege, Geschwindigkeiten und Beschleunigungen Abb. 4. Geschwindigkeiten und Beschleunigungen 
in Abhangigkeit von der Zeit fiir 7,—6. in Abhingigkeit vom Wege fiir 7) = 6. 


Nach Trennung der Veranderlichen und Ausfiihrung der Integration nach der Zeit in (12a) 
ergibt sich 


Er» — 
1 ( y/m, + 2¢8,-,, 
(ig ey ee Go eS 18 
ye | aoa ~ 
Man setzt &, = &) cos yin (18) ein und erhalt 
ie Se lala 
dike = 1/222 | ya ide 1 
ue [EE [ath Fora " 
Pr 
Nach weiteren Umfo: mungen ]a8t sich schreiben 
eee ew, 
ota) SSE | JIB ainty dp ; (20) 
7 Yo 
darin sind 
4aé 
ae [cn ea Sel Oa 21 
GS Og Zp 0< ory Ge (21) 


Der Ausschlag sei willkiirlich als positiv zu betrachten, wenn die reduzierte Lange L mit wachsen- 
dem Ausschlag zunimmt und umgekehrt. Dann wird die Schwingungszeit fiir die Schwingungs- 
abschnitt von Null bis &,) und wieder nach Null zuriick mit der iiblichen Darstellung fiir ein ellip- 
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tisches Integral zweiter Art A 

Im, +2 0&1 ( Aad S 29 
7, = 4/™ teste Bl iB ( ) 
Entsprechend wird die Zeit fiir die Halbschwingung der entgegengesetzten R'chtung von Null 

bis — ,, und wieder nach Null zuriick 
g/m + 24 S10 ( 4a "| 23 
fe p V Cy Rites Ie +2a&,° 4 ai 
Die Addition von (22) und (23) liefert die Schwingungsdauer einer vollstandigen Schwingung 


i aa d Mee a ler <= oo E , (24) 
1 


worin E ein vollstandiges, elliptisches Integral zweiter Gattung ist. T* ist absichtlich eingefiihrt 
um von dem Fall unterschieden werden zu kénnen, 


E ee . bei dem die reduzierte Lange wahrend der Schwin- 
gung konstant bleibt. 
Tene 
4 600 
6 { AZ, 
16 \-a 
+ 400 
44 — ny | 
x oe 
a 13 = | 
SII he: 45 9% 
ha 
sl! = 
0,52 | 
——— -400 an 
~600 ve 
Te 
* sky sab ER 3 ¥ 5 6 7 
Abb. 5. Die maximale Beschleunigung in Abhangigkeit Abb. 6. T,,,/T,p, in Abhangigkeit von 7, mit Anfangs- 


von den gréBten Auslenkungen mit 7, als Parameter. auslenkung €,, als Parameter. 


Durch Reihenentwicklung des Integranden der Gleichung (24) und mit M = 


erhalt man 


1-3-5-+-(2n—1) 
2-4-6225 Ore 


vere m,+ 2a Ey) 7 M2 Aa Ex y 
2 =2n)/ , r pig i, te Gere ailiaa fe 
Aus (22) und (23) lassen sich noch die folgenden Verhaltnisse aufstellen: 
E—E(,,™ : 
Fe eee ( 4 ) ste eee (26) 
Tu ae a) Tu E(k ae Tab [Di 15 k, Oe i 
4 “4 4 


Da nach der Eigenschaft der elliptischen Integrale stets 
Tt m 
E(k. 7) < E<2 E(k 7) 
ist, so ergibt sich noch die folgende Tatsache: 


are Ih, (giiltig fiir fo a 1), (27) 


d.h. die Schwingungszeit in der Richtung, in der die reduzierte Lange der Rohrleitung mit Zu- 
nahme der Auslenkung wachst, ist gréBer (Abb. 6), als die in der entgegengesetzten Richtung. 
Je mehr sich die beiden Endquerschnitte voneinander unterscheiden, desto gréBer ist die Ab- 
weichung der Schwingungszeiten in beiden Richtungen voneinander. AuBerdem wird auch diese 
Abweichung von der VergréBerung der Anfangsauslenkung verstarkt. 


XX. Band 1952. Hsien-Chih Liu: Eigenschwingungen idealer Flissigkeitenin Rohrleitungen usw. 307 


c) Vergleich zwischen T* und T. T sei die Schwingungsdauer eines ahnlichen Systems, 
das sich aus einem Kolben mit der Masse m, einer Feder mit der Federkonstanten c und einer 
Fliissigkeitssaule von konstantem Querschnitt mit einer Lange L, zusammensetzt. Die Schwin- 
gungsdauer fiir ein solches System ist offenbar unabhangig von der GréBe der Amplitude. Die 
Schwingungsdauer T' dieses vereinfachten Systems soll als Vergleichswert fiir den Abschnitt b) 
behandelten Fall herangezogen werden. Die Neigungen der Rohrenden im behandelten und im 
vereinfachten System sollen an den Stellen 1 bzw. 2 jeweils einander gleich sein. Fiir den verein- 
fachten Fall ist 


; = 
T=2n)/ ; m + oF, Ly . 
c+ y F, (cos x, — cos «,) 


und fiir den behendelten Fall ist 


T* — 97 | mtehi—D So | 1 — S I m2 (_2eFi G2 — 1) S10 ‘4 ; 
cy ee m, + @F, (nj — 1) &s 


Dann erhalt man 
ee = c + y F, (cos a, — cos &,) y oF, (73 — 1) y Fh, (28) 
/E,,>0 c +» F, (7, cos «, — cos &,) c,(m+oF,Iy) ~ 


Mit den Bedingungen 7, > 1, ¢ +y F; (7? cos x, — cos a1) > 0 und c +» F, (cos a, — cos 0%) > 0 
kann man aus (28) die Giiltigkeit folgender Beziehung feststellen: 


Tee eal (29) 


: ‘ Ons ; 
Die Bedingung aE, = 0 hefert 
1 


m, = n 4aé = 1 ; 4aé ic 

it 4 a M2 $10 2 10 a 
aA Sa tees eas jie ese ae) As 
mit der Lésung &,) = 0. Mit &,) = 01aBt sich bestatigen, daB 


(az _ a FH(ni—1P 
05To Ex5=0 4c, mi 


n—1 


Hiermit ist bestatigt, da T* bei &,, = 0 ein Maximum besitzt. Infolgedessen ist damit auch 
bewiesen, da prinzipiell stets 

TX 7, (30) 
ist. Die (T* &,))-Kennlinie besitzt bei &;, = 0 eine waagerechte Tangente und fallt mit wachsen- 
dem £,, monoton gegen die &,)-Achse. 

d) Verlauf von T,,und T,,inAbhangigkeit von &,. Aus der Beziehung 
ae _ @F, (4; —1) 
O10 F,o—0 im, ey 

ist ersichtlich, daB die (Ty, , &,9)-Kurve bei €;) = 0 nur horizontal verlaufen kann, wenn 7, = I, 


d.h. F, = F, ist. Ist 7, >1, so steigt die Kurve mit wachsendem &;, , sie fallt, wenn 7.< 1. Der 
Verlauf von T,,, kann von der Differenz T* — T,,, leicht abgeleitet werden. 


e) Vergleichzwischen T/2 und (T,,)z,,—9. Aus Gleichung (22) erhalt man 


a m + QF, Ly tel? es 3] 
(Teu)e—0 = 2 c+yF,(y, cosa, —coso,) [e+ y Fy (7, cos x, — cos o,)]? ep 


Daraus 1aBt sich schlieBen, dab 


(Tees, (32) 


wird. 

6. Zusammenfassung. Die Darstellung der Schwingung einer idealen Flissigkeit in Rohr- 
leitungen beliebiger Gestalt fihrt zu nichtlinearen Differentialgleichungen zweiter Ordnung; bei 
diesen ist der Koeffizient des zweiten Differential quotienten eine lineare Funktion des Ausschlages. 
Die Differentialgleichung ist wegen des Vorhandenseins eines quadratischen Gliedes des ersten 

22 
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Differentialquotienten nicht linear. Das Verhaltnis der beiden Leitungsquerschnitte an den 
Enden der Fliissigkeitssaule tritt als ein wichtiger Parameter in saimtlichen Koeffizienten der 
Bewegungsgleichung auf. 

Die oben erwahnte Differentialgleichung gewinnt wieder ihre lineare Eigenschaft, wenn und 
nur wenn die beiden Leitungsquerschnitte an den Enden der Fliissigkeitssaéule gleich macht. 
Wenn die Leitungsquerschnitte an den Enden der Fliissigkeitssiule verschieden groB aber kon- 
stant sind, laBt sich die Differentialgleichung mathematisch streng lésen. Die Schwingungszeiten 
in beiden Richtungen weichen voneinander ab. Diese Abweichung wird von der VergréBerung 
des Ausschlags monoton verstarkt. Die Schwingungszeit fiir die Richtung, in der die reduzierte 
Lange mit wachsendem Ausschlag zunimmt, iiberwiegt unter Umstanden um ein betrachtliches 


MaB. 


(Eingegangen am 25, Januar 1952.) 
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Uber Schwingungen einer idealen Fliissigkeit in ellipsen- 
und kreishogenférmigen, kommunizierenden Réhren. 


Von Hsien Chih Liu. 


1. Einleitung. Das Ziel der vorliegenden Arbeit ist die Untersuchung von Schwingungen einer 
idealen Fliissigkeit in ellipsen- und kreisfirmig gebogenen Roéhren von konstantem Querschnitt. 
Der erste Fall wird behandelt mit der Einschriinkung, da®B die beiden Enden der Flissigkeitssaule 
in der Gleichgewichtslage mit einem Hauptdurchmesser der ellipsenférmig gebogenen Réhre zu- 
sammenfallen. Fiir beide Falle gilt die Annahme, da®B die Endflachen der Fliissigkeitssdule stets 
senkrecht zur Mittellinie der Réhre sind. 


2. Schwingungen einer idealen Fliissigkeit in ellipsenformiger, kommunizierender Réhre von 
konstantem Querschnitt. Es sei angenommen, daf die ellipsenformig gebogene Riohre konstanten 
Querschnitt besitzt. Man wird aus Nachfolgendem ersehen, 
daf diese Kinschrankung gar nicht notwendig ist. Es geniigt, 
daB die beiden Enden der Fliissigkeitssaule sich wihrend der 
Bewegung in Rohrabschnitten mit konstanten Querschnitten 
bewegen. Die Gestaltung des dazwischen liegenden Teiles 
der Roéhre ist vollig belanglos, wenn man die Gestaltung des 
Rohres bei Errechnung der reduzierten Lange der Fliissigkeits- 
sdule entsprechend beriicksichtigt. Die beiden Enden der 
Fliissigkeitssaule sind auch wahrend der Bewegung gleichem 
Druck unterworfen. Die Fliissigkeit sei auBerdem inkom- 
pressibel; Kavitation soll nicht auftreten. . ee ne en trae pases 

Nach der verallgemeinerten Bernoullischen Gleichung _ jn chiipsenférmiger, kommunizierender Réhre 
erhalt man mit den Bezeichnungen von Abb. 1 be eee eae 


Con = 1 
[5a t2er=9, (1) 


$1 


wobei y der Abstand des Fliissigkeitsspiegels von der Gleichgewichtslage und v die Geschwindig- 


keit der Strémung an einer beliebigen Stelle s bedeutenten. Ferner gilt 


$2 P2 F. 
te | = | je sin? p + b? cos? p dp (2) 
51 1 
und 
dv ds as (dp\?_, Os dy 3) 
dt dt® dg*\dt Og dt® ’ 


hierin bezeichnen L die Linge der Fliissigkeitssiule zwischen irgend zwei Grenzstellen s; und sy, 


g den Winkel entsprechend Abb. 1 und t die Zeit. Das Einsetzen von (2), (3) und in (1) liefert die 
Bewegungsgleichung 


2 Ip\2 , 2¢b.;>-—G— ute: $e 
(a? sin? p + 6? bert ge) ee (a? — 6?) sin p cos 7 (a + T. \@ sin?y+ bcos*ysing=—0, (A) 


dt? 
wobei L, die Lange der Fliissigkeitssdule bedeutet, die gemaB unserer Voraussetzung die Halfte 


des Ellipsenumfangs betragt. Es laBt sich zeigen, daB die Beschleunigung nicht linear von p 
abhangt, sie aber zu beiden Seiten der Ruhelage bei gleichem Abstand, gleichen absoluten Wert 


hat. Nach Umformung 1aBt sich Gleichung (4) mit 
_1-3-5++-(2n—1) 


ye 2-4-6---2n 
schreiben " 
' \2 Sig oo 2__ p2\" cos2"t! 
d \a?sin® py + b? cos? y (dp 286 Diy oe 2gabwy 1 m(* le 0X6 
emcee ft orpcegta i oe) aa as 
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oder nach Integration und Bestimmung der Integrationskonstanten durch die Randbedingung 
Y = Go - dp/dt = 0 fiir eine volle Schwingung 


Po 
° 


(a? sin?» + b? cos? p dy . (7) 


= 00 l (PIP ” cos? "+1 y — cos?" t1 gy, 
COS P — COS Po Deen | = re 


Die Integration der rechten Seite von (7) ist als allgemeiner Ausdruck schwierig auszufiihren; sie 
kann aber im Prinzip als gelést betrachtet werden. 


0 


3. Schwingungen idealer Fliissigkeiten in kreisférmiger Réhre mit zwei verschieden groBen 
konstanten Querschnitten. Nach der erweiterten Bernoullischen Gleichung gilt fiir die beiden 
Enden der Fliissigkeitssaule (Abb. 2), wegen 
PP; = Pp = Atmospharendruck 


2 yf Lig B 
vmod yy, m=z | ee: (8) 


28 


In Abb. 2 ist y, der Abstand der Flissigkeits- 
spiegel bei Gleichgewicht von der Bezugsebene 
0 — 0, und gy, der dazu gehérige Winkel. In 
Gleichung (8) bezeichnen v, und v, die Spiegel- 
geschwindigkeiten in der tangentialen Rich- 
tung des Kreises, v die Geschwindigkeit an 
einer beliebigen Stelle s zwischen den beiden 
Spiegeln, y, und y, die lotrechten Abstande 
der Spiegel von der Bezugsebene und ds ein 


Abb, 2. Zur Schwingung idealer Flissigkeit in kreisformiger Réhre - 2 hs 
mit zwei ungleich groBen konstanten Querschnitten. In der Strémungsrichtung gemessenes Langen- 


element. Ferner folgt wegen 


F : 
is UF, = v2 F;, y1 = Rsing, 


Po = 17 (91 —Po)> = 2 =E— Pot NoV1—%o)> Yo = Rsin(Y — 9) 
aus (8) 


oe (1 — 73) + R[sin y, + cos Go sin Ho(~ — Yo) — sin gy cos No (Y, — %)] = = i dv ds. (9) 


Unter Beachtung von 


ds, 


Y=: — 
a dt’ 


oF =»v,F,, i 


wird 


2 
zs) 


1 ds,\? ; 1a ise 
L Aas (1 — 95) (zi) + Risin y+ cos Ho sin H2(P1 — Yo) — sin Yo COS 1) (Y1 — Yo) ] = 0. (10) 


Bei der Voraussetzung F, == konst. und F’, = konst. gilt 


hierin ist 
Lo = / 4 ds = konst., Jn Mea oO In Osi sa 
0 0 
Demnach erhalt man 
LS ii 4 ds = L, — (v2? — 1) s,, 


sy 


wobei Ly die reduzierte Linge der Fliissigkeitssaule in der Gleichgewichtslage bedeutet. Wenn man 
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die Beziehung s, = R(p, — Yo) beachtet, so erhalt man mit 2 = nz—1 


d* 1 dp 
(Ly +4 Ry —A Rq,) ie . AR( z) | g {sin My ++ sin [2 P1 — (1+ 2) Gol} =(. (11) 


Dies ist die gesuchte Bewegungsgleichung. Es laBt sich zeigen, das die Beschleunigungen anj e 
zwei zur Ruhelage symmetrischen Stellen verschiedene Werte haben. 


Man kann (11) in die folgende Form umschreiben: 
d [Lo tARG ry 3 ae a 2 


dt | 2 21 Vi — & [2 cos py + COS (1) Py s—= (1 75) #)]) =0 (12) 


oder 
f dy,\? , , ; 
No(LotARGy—ARQ) (ze) + 2 2 {N2 (COs Pye — cos Hy) + cos [479 Yq — (1+ 72) Vol | 


— cos[1»~, — (1+ me) Pol} =0, 
worin die Integrationskonstante durch die Randbedingung dy,/dt = 0 bei v, = 9 bestimmt 
wurde. Nach Umformung erhalt man 


(13) 


is a, 


& ee: | be Eg FAR G—ARG, 
| 


—dg,. (14) 


Der Integrand der rechten Seite von (14) ist eine ungerade Funktion von ¢, und das bedeutet, daB 
die Schwingungszeit fiir entgegengesetzte Richtungen verschieden sind. Die Integration fiir 
praktische Falle ist stets numerisch oder numerisch-graphisch auszufiihren. 


4. Schwingung idealer Fliissigkeit in kreisbogenformigen Réhren von konstantem Querschnitt. 
Wenn man 77. = | setzt und den dadurch belanglos gewordenen Index @ fallen 1aBt, erhalt man 


aus Gleichung (11) die Bewegungsgleichung fiir diesen Fall 
Ly 5! + glsing + sin(y —294)] =0- (15) 
Nach der ersten Integration ergibt “ek 
16(ir) + 2p[cos g, — cos p+ cos (Pa — 2 Fo) — cos (p — 2 )] = 0, (16) 


wobei gy, der Auslenkungswinkel am Anfang mit 0—0, als Bezugsebene bedeutet. Mit den Sub- 
stitutionen 


Pa—Yo A) 


if == 590 eae sin a =ksin yp 


14350. Cnet) 
Redo one 


erhalt man die Zeit einer Vollschwingung mit M — 


4 


Pe 
tas Vi g : 


ara maa 


0 | as Aa =22)/, 1+ >) M?sinte OP) (17) 
cos Py J V1—K* sin? yp 2 g COS Yo — : 


Fiir den Spezialfall vy, = 0, d. h. wenn die Enden der Fliissigkeitssaule in der Gleichgewichts- 
lage gerade mit dem waagerechten Durchmesser des Kreises zusammenfallen, erhalt man aus (17) 


4 [Ly “ 2 otn2n Qa 18 
eee at oY sin?” | . (18) 
LaBt man den Durchmesser des Kreises unbeschrankt wachsen, dann wird gy, = Sy) OTund 


damit aus (18) 


Das entspricht auch dem Newtonschen Resultat fitr Flissigkeitsschwingung in U-férmig gebo- 
gener Réhre mit konstantem Querschnitt. 
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Fir kreisbogenférmige Réhren kann man noch Ly = R(a-+ 2%) schreiben, und somit wird | 


rman) ee 3 Mes), (19) 


& COS Po i 


oder fiir den Sonderfall gy, = 0 


T=2n e() ae > Me? ants) (20) | 


ios l 


5. Vergleich zwischen Schwingung eines mathematischen Pendels und einer idealen Fliissigkeit 


in kreisbogenférmigen Réhren von konstantem Querschnitt. Die bisherigen Ergebnisse behalten — 


ihre Giiltigkeit ohne weiteres auch fiir solche Roéhren bei, die mehrere Windungen besitzen. Man 
hat nur dann L, = R(2v2-+-% +2) zu setzen mit y= 0, 1, 2,..., und erhalt damit die 
Schwingungsdauer fiir den allgemeinsten Fall 


R(2vx+2+2 Mp) Slag t en Veo 
va Selbstverstandlich miissen die vollen 
fs aes Windungen der Réohre vollstandig 
[| | mit Fliissigkeit gefillt sein. 
a 1S ‘ Die in Abb. 3 aufgetragenen Kur- 
Ee | t Ni | ven sind errechnet aus (21) mit y= 0. 
/ NN | Da ein AbreiBen der Strémung nicht 
/ y | auftreten soll, galt die EKinschran- 
/ Ni | kung q@ <2/2. Mit wachsendem 
‘es if Ni | Filllungswinkel wachst die Schwin- 
Si; gungszeit und sinkt die maximal er- 


reichbare Amplitude. Die Schwin- 
gungsdauer wachst auch mit der Zu- 
nahme von dem Ausschlag @, . 


Wenn nur der Fall »y = 0 betrach- 


‘nach der prositiven Richting 


/ 


| S 
Pa f2 pon P42 Po 
Re aCe : (22) 


-240° -210 -780 ~10 -10 -90 -60 -3 0 JO WO 9° 


Abb. 3. Schwingungsdauer einer idealen Fliissigkeit in kreisférmiger Rébhre 
von konstantem Querschnitt in Abhingigkeit von dem Winkel gp. 


ic 


und fiir das mathematische Pendel 


Le = Ao bi ee 
: (2 ! Den 5): (23) 


i y, 


Was aa 2 It 


wobei 2. der vollstandige Schwingungswinkel und | die Pendellange bedeuten. Der Winkel 
Pa —Yo in (22) entspricht dem Winkel a =p, —p in (23). Fiir den Vergleich soll deshalb 
X = Py — Yo gesetzt werden, d. h. die beiden Schwingungssysteme machen den gleichen Winkel- 
ausschlag. Es ist hiernach 


F nt l 
(ies GE fiir (= - v6) sun COS Do. (24) 
Der interessanteste Fall von (24) ist 
: R 
TT st TS +9) Sno. (25) 
Die zweite Gleichung (25) ist in Abb. 4 mit R/l als Parameter dargestellt. Die Schnittpunkte der 


Geraden mit der Kosinuslinie sind diejenigen q-Werte, bei denen beide Schwingungssysteme 
gleiche Schwingungsdauer haben. Aus der oberen Grenze fiir E folgt die Kinschrankung 


SQ : ; ae oan | 
S tet wird, so ist fiir die Fliissigkeits- 
& saule 
Ss = nee 
Ss 

: | 2 g COS Po 
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fiir das Verhaltnis 1/R 
1 l 


a Ut 1 ey < i 
oe ee ee a Os “0 
0& 1: 
Da (sea) fs | ist, besteht nur fiir R</ die Méglichkeit, daB die beiden Schwingungssysteme 
Po=—z/ 


gleiche Schwingungsdauer besitzen. 
Fir @) = 0, y.—> 0 und a— 0 gelten die folgenden Benepe: 


Es ergibt sich die itherraschende Liésung 
Ta=Thro (giiltigfirg,—->0, o—+O0undzR=21). (28) 
Mit anderen Worten: bei sehr kleinen Schwingungsausschlagen besitzen beide Schwingungs- 


systeme gleiche Schwingungsdauer, wenn die Pendellange gleich dem vierten Teil des zugehérigen 
Kreisrohrumfangs ist. Die Pendellange ist dabei gleich der Linge der halben Fliissigkeitssaule. 


6. Vergleich der Schwingungsdauer einer idealen Fliissigkeit in U-férmig gebogener und kreis- 
bogenformiger Réhre. Aus T = 27 VL/2 g und (19) erhalt man 


-R@ +2 Po) 
Ess Ao) =e ecg wi 
Mithin laBt sich schlieBen, dab, wenn die exakten Langen der Fliissigkeitssdule beider Schwin- 
gungssysteme einander gleich sind, die Schwingung in U-formig gebogener Réhre wenigere Zeit 
benétigt. Fir ~ = 0 und g,— 0 ent- 
E steht der Grenzfall, daB die Schwingung 
der Fliissigkeit in kreisbogenférmiger 


a 
\/i & WY ss : : areas 
INS aA Rohre der Schwingung in U-férmiger 
fare Réhre gleich ist. 
G75, 
02 
O5- 
G25 
= — %o 
Te _% Ea Ed o 
2 a 0 a D- 
Abb. 4. Die g,-Werte, bei denen beide Schwingungssysteme gleiche Abb. 5. Schwerpunkt eines 
Schwingungsdauer besitzen. Kreisringsegmentes. 


7. Vergleich zwischen der Schwingung eines physischen Pendels und einer idealen Fliissigkeit in 
kreisbogenférmigen Réhren von konstantem Querschnitt. Die Schwingungsdauer eines physischen 
Pendels ist 

jee ae Ja_ (1 +S Misinte 9), (29) 
worin J, das Massentragheitsmoment des Koérpers um cine horizontale Drehachse durch einen 
Aufhangepunkt, m die Masse des Korpers, e der Abstand des Schwerpunktes vom Aufhangepunkt 
und ¢, der maximale Ausschlagswinkel bedeuten. 

Wenn der Kérper die Form eines Ringsegmentes (Abb. 5) hat, dann wird e = R sin y/y und 
J, = mR? sein, worin y die Halfte des Zentriwinkels und R der mittlere Radius des Ringsegmentes 
sind. Fiir diesen Fall kann man fiir (29) schreiben 


= Ee — 2 ann £9 30 
PAV eare (+ >oM sin }: (30) 

Wenn dasselbe Ringsegment aus Fliissigkeit statt aus festem Stoff besteht, dann erhalt man 
mit den Resultaten der Rise SSE SEE Abschnitte 


= (R(x + 2%) ee) 31 
Aa | “2g C08 Yo (+37 sa 3): GD 


roel 
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Wenn man in Gleichung (31) 2 +2q) =2y und cos Y = cos (y —2/2) = sin y setzt, dann wird 


A Ry < 2 otn2an £0 39 
Die Gleichungen (30) und (32) sind identisch, d.h. eine kreisférmige Flissigkeitssdule und ein 
gleich groBes Ringsegment aus festem Stoff schwingen mit gleicher Schwingungsdauer. Die 
Schwingungsdauer ist nicht von der Dicke sondern nur von dem Radius und der Bogenlange des 
Ringsegmentes abhangig. 


8. Betrachtung der Grenzfille. Wenn man in Gleichung (31) @y)> + 2/2 gehen 1aBt, so wird 
Ty, 9, —nj2 = © - (33) 


Dieses Ergebnis ist plausibel, weil @ = +-7/2 eine Vollringréhre darstellt, die mit Flissigkeit 
vollstandig gefwllt ist. In diesem Grenzfall kann die Fliissigkeitssdule iberhaupt nicht mehr 
schwingen, oder héchstens kann sie mit einer unendlich groBen Schwingungszeit schwingen. 
Diese Betrachtung scheint natiirlich auch richtig aber man darf nicht vergessen, daB man die 
Erscheinung stets im Sinne der Schwingung zu deuten hat. Gleichung (29) zeigt, daB bei einem 
Vollring der Schwerpunktsabstand eim Nenner verschwindet und den Ausdruck zu Unendlichem 
fihrt. Vom Standpunkt der Schwingung gesehen ist dieses System gegen kleine Stérungen in- 
different, d. h. die fiir eine Schwingung notwendige Riickstellkraft ist in diesem Fall nicht mehr 
vorhanden, weil jetzt der Kérper in seinem Schwerpunkt gestiitzt ist. 
La8®t man anderseits in Gleichung (31) yp — 2/2 gehen, so wird 


af IR. Sere + 2/2 
(Pa a ey | _ ( oF > M2? sin?” are (34) 
Ww j 


Dieses Ergebnis stimmt mit der Schwingungsdauer eines mathematischen Pendels iiberein. Die 
Masse konzentriert sich bei g—> — 2/2 in einem Punkt. 


9. Zusammenfassung. Die Schwingungszeit idealer Fliissigkeiten in elliptischen Réhren mit 
konstantem Querschnitt und in kreisbogenférmigen Réhren von ungleichen, konstanten Quer- 
schnitten 1aBt sich numerisch errechnen. 

Die Schwingung idealer Fliissigkeit in kreisformigen, kommunizierenden Réhren von kon- 
stantem Querschnitt wurde mathematisch streng gelést. 


(Eingegangen am 25, Januar 1952.) 
Anschrift des Verfassers: Hsien-Chih Liu, Tsinan (China), Technische Hochschule. 
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Nichtlineare Biegeschwingungen des axial-pulsierend 
belasteten Stabes. 


Von F. Weidenhammer. 


1. Zielsetzung. Zahlreiche neuere Arbeiten untersuchen die Stabilitat einfacher Grund- 
bewegungen elastischer Kérper unter pulsierenden duferen Lasten, und inzwischen ist von 
E. Mettler * eine sehr allgemeine Theorie dieser kinetischen Stabilitatsprobleme auf der Grund- 
lage der Methode der kleinen Schwingungen geschaffen worden. Man kann bei diesen Problemen 
stets gewisse Frequenzen der pulsierenden Lasten abhangig von deren Amplitude angeben, die 
zu kinetischen Instabilitéten einer Grundbewegung fiihren. Die Theorie liefert fiir diese aus- 
gezeichneten Werte von Frequenz und Amplitude streng nur die Aussage, daB die Grund- 
bewegung instabil ist und sagt iiber die dann eintretende Bewegungsform selbst nichts aus. Als 
bevorzugtes Modell fiir diese ganze Problemklasse ist mehrfach der axial-pulsierend belastete 
Stab behandelt worden. In diesem Beispiel ist die Grundbewegung die reine Langsschwingung 
des Stabes in seiner Ausgangslage, und die Theorie sagt nicht befriedigend aus, in welche stabile 
Bewegungsform die genannte instabile Grundbewegung infolge einer kleinen Stérung iibergehen 
kann. Da die in der Technik vorkommenden pulsierenden Lasten klein gegen die statische Knick- 
last sind, wird man keine unbegrenzt wachsende Resonanzlésung erwarten, sondern eine Be- 
wegung mit endlicher Amplitude vermuten. Es miissen daher gewisse stabile, noch unbekannte 
Nachbarbewegungen vorhanden sein. Zweck der vorliegenden Arbeit ist es, einen Beitrag zur 
Auffindung dieser stabilen Nachbarbewegungen zu liefern. Dabei beschrankt sich die Unter- 
suchung auf den technisch wichtigen, axial pulsierend belasteten Stab mit freiaufgestiitzt ge- 
lagerten Enden. An den Stabenden soll also der Stab in solcher Weise gelenkig gelagert sein, daB 
dort keine Querverschiebung und kein Biegemoment auftritt. Hingegen mu ein Lager in Langs- 
richtung verschieblich sein (vg]l. Abb.1). 

Da die mathematische Beschreibung der Nachbar- 
bewegungen auf nichtlineare partielle Differential- 
gleichungen fiihrt, darf man zunachst nur bescheidene 
Naherungsaussagen erwarten. Immerhin lassen sich 
gewisse Aussagen iiber endlich bleibende Stabbewe- 


gungen bei instabiler Grundbewegung machen. 


70(2,t) 


Plt) | wm Plt) @& 


Abb. 1. Der Stab im Koordinatensystem. 


2. Allgemeine Stabilitatsaussagen. Die allgemeine Stabilitatstheorie gestattet fiir axial har- 
monisch pulsierend belastete Stabe folgende Aussagen. Wenn die Frequenz w der pulsierenden 
Axiallast weit unterhalb der tiefsten Eigenfrequenz der Langsschwingungen liegt, ist die Grund- 
bewegung, das ist hier die reine Langsschwingung in der gestreckten Ausgangslage des Stabes, 
im allgemeinen stabil. Instabil kann dieser Grundzustand nur in der unmittelbaren Nahe dis- 
kreter Frequenzen j, werden. Und zwar findet man fiir den freiaufgestiitzten Stab nur Fre- 
quenzen in der Nahe von 


_ 203 


(p> 0 ganz), (1) 


worin «; die i-te Biegeeigenfrequenz des Stabes hezeichnet. Diejenigen Anregungsfrequenzen We, 
welche gerade noch zu Instabilitaten fithren, hangen von der Amplitude der pulsierenden Last 
ab. Sie schlieBen einen Frequenzbereich ein, den man Bereich erster Art der Ordnung p nennt, 
wenn er von Frequenzen (1) ausgeht. pari teeens 

Der komplizierter gelagerte, z. B. eingespannte Stab hat iibrigens weitere zusatzliche instabile 
Anregungsfrequenzen in der Nahe von 

wi + Oj : ae 
Wy = or (p>0 ganz, ij). (2) 

Diese sogenannten Instabilitatshereiche zweiter Art in der Nahe der Frequenzen (2) kénnen bei 
dem freiaufgestiitzt gelagerten Stab nicht auftreten und bieten daher hier kein Interesse. 


1 FE. Mettler, Ing.-Arch. 17 (1949), $.418 und die dort angegebene Literatur. 
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Weiter entnimmt man der allgemeinen Theorie und zahlreichen Zahlenbeispielen, daB die Be- 
reiche erster Art erster Ordnung die bei weitem gréBte praktische Bedeutung haben. Das folgt 
cinmal daraus, daB diese Anregungsfrequenzen von jeweils etwa der doppelten Biegeeigentrequenz 
eine besonders starke Instabilitat ergeben, und daher die Dampfung bei diesen Frequenzen nur 
geringere Wirkung hat als bei den anderen Bereichen. Zum anderen ist die Frequenzbreite dieser 
Bereiche die bei weitem gréBte. Man wird daher sein Hauptinteresse auf die Instabilitaten bei 
Langsschwingungen des Stabes mit Anregungsfrequenzen in der Nahe von Ww.) = 2; richten. 


3. Ausgangsgleichungen. Die systematische und auch die naherungsweise Behandlung wird 
zweckmabig tiber das Hamiltonsche Prinzip gefithrt. Demnach verlauft die Stabbewegung so, daB 


ty ty 
Saf CAG Veer) dee [oAdt=0 (3) 


wird. Hierin bedeuten A; die Formanderungsarbeit, V das Potential der duBeren Krafte, T die 
Bewegungsenergie und 6A die virtuelle Arbeit der nichtkonservativen Dampfungskrafte. Mit den 
in der technischen Biegelehre bewahrten, vereinfachenden Annahmen ist 1 


1/2 1/2 


A; = se | [u, + a dx + 2e i Wi, 0% , Viz Pine & , 


2 
Beye iP . 
1/2 1/2 (4) 
T=% | (uitui)de, v4 | Pujdwde. | 
~—i/2 —1/2 
Hierin bezeichnen (vgl. Abb.1) 
% = die Ortskoordinate, 
u(x, t), w(x, t) = die gesuchten Verschicbungskomponenten in Stablangsrichtung und quer zur 
Stabachse,. 
t = die Zeit, 
J = konst. — das fiir die Biegeschwingung maBgebende Flachentragheitsmoment, 
F = konst. = die Querschnittsflache des Stabes, 
js = konst. = die Stabmasse je Langeneinheit, 
l = die Stablange, 
E = den Elastizitatsmodul, 
P(t) = P, + P,cos wt die rein harmonisch pulsierende Axiallast, 
i) = die Dampfungskonstante der Biegeschwingungen. 


Der erste Summand in der Formanderungsenergie (4) riihrt von der Dehnungsarbeit her. Nur 
durch diesen Energieanteil entsteht in der obigen Naherung fiir A; eine Kopplung der Langs- 


und Querbewegung des Stabes. Man iiberzeugt sich leicht aus der Anschauung, da = w nicht 


gegen u, vernachlassigt werden darf, wenn man die vorliegenden Instabilitatsfragen untersuchen 
will. Ein Element der Stabachse habe namlich vor der Ver- 

ad p formung die Linge dx, nachher dx. .Dann ist die Dehnung: 

S a(vru) va Ey, = dx/dx —1*. Aus der Abb.2 liest man _ ferner 


Abb. 2. Der Stab mit einem Langenelement dt = d(x w)* dw? ab und findet mithin 
der verformten Stabachse und seinen ee 


Komponenten. = VQ +4u,) | we le ux +5 w, . 


Die niedrigste Naherung verlangt also bereits die Mitfiihrung des quadratischen Gliedes a Ware 


Unter Voraussetzung des Hookeschen Gesetzes ergibt sich dann aus 4 Je dV (dV Volumelement) 


und damit der erste Summand von (4). Der zweite Summand von (4) rithrt von der Formande- 
rungsarbeit der Biegemomente her und kann bei kleinen Langskraften noch unabhangig von den 
Langsverschiebungen berechnet werden. 

Die Voraussetzung der strengen Giiltigkeit des Hookeschen (linearen) Elastizitatsgesetzes be- 
darf im vorliegenden Falle einer Rechtfertigung. Diese 1aBt sich bis zu einem gewissen Grade 


1 E. Mettler, Ing.-Arch. 13 (1942), S. 97. 
2 Vel. z. B. A. Pfliiger, Stabilitatsprobleme der Elastostatik, S. 10. Berlin 1950. 
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aus der Durchrechnung eines benachbarten Problemes gewinnen 1. Es zeigt sich namlich bei der 
Berechnung der freien nichtlinearen Biegeschwingungen eines beidseitig unverschieblich ge- 
lagerten Stabes, da infolge eines nichtlinearen Werkstoffgesetzes im allgemeinen zwar die 
Koeffizienten der nichtlinearen Schwingungsgleichung etwas andere Werte annehmen, daB je- 
doch bei normalen technischen Werkstoffen und Abmessungen der Charakter der Schwingungs- 
gleichung nicht gedndert wird, so da keine andersartigen Phinomene auftreten, wie sie der 
vereinfachte Ansatz (4) liefert. Dieser Ansatz soll daher auch fiir den pulsierend belasteten Stab 
beibehalten werden. 

Aus dem Hamiltonschen Prinzip (3) lassen sich nun die gekoppelten Bewegungsgleichungen 
fiir die Langs- und Querbewegungen des Stabes herleiten. Man erhalt gekoppelte nichtlineare 
partielle Differentialgleichungen mit nichtlinearen und wegen der pulsierenden Liangskraft zeit- 
abhangigen Randbedingungen. 


Im einzelnen findet man fiir die Langsbewegung 


—EF (u, - = of fy = 0 (5) 
mit 
Tye ty (DT Rai g x= —1/2, (5’) 
a 
EF(u, oe ee ee 0) Stir > are /2, (5’’) 
Fiir die Querverschiebung gilt 
“2h FR he, (u, ae | + E J Wyxxy + Bw, Wy = 0 (6) 
mit ‘4 
1 — WO LED a. (6’) 


Die strenge Liésung der vorstehenden Differentialgleichungen scheint kaum méglich zu sein. 
Man wird daher zum naherungsweisen Studium des zeitlichen Verlaufes der Bewegungen zweck- 
maBig die schon vielfach bewahrten ,,gemischten Ritz-Ansatze‘‘ * 


Hes ee demu Oe (7) 


i=1 


wv 
w(x.) = 3 X,(x) Mi (3) 
ii 
direkt in (3) einfiithren. In (7) und (8) sind die (dimensionslosen) Ortsfunktionen Y; und X; so 
zu wahlen, daB alle wesentlichen Randbedingungen fiir alle Zeit erfiilit sind. Zur Erleichterung 
der Zahlenrechnung werden Y; und X; als Funktionen der dimensionslosen Ortskoordinate 


ee 1 (9) 
durch die Bedingungen 


[ Gdiaw, f Yidg=1 (10) 
—1/2 1 


normiert. 

Man hofft, in den Ansitzen (7) und (8) mit einer méglichst kleinen Zahl N der Ansatzglieder 
fiir eine brauchbare Naherung auszukommen. Zur Erfiillung von (3) hat man eine Variation 
auszufiihren, und diese fiihrt auf nichtlineare Eulersche Gleichungen fiir die Zeitfunktionen U; 
und W;. Durch die Ansitze (7) und (8) kommt man so von partiellen zu nichtlinearen gewohn- 
lichen Differentialgleichungssystemen, iiber deren Lésungen man gewisse Aussagen machen kann. 

Um in der weiteren Rechnung weitgehend unabhangig von den speziellen Stababmessungen 
zu werden, sollen zunichst geeignete Parameter und Veranderliche cingefiihrt werden. Auber 
der dimensionslosen Ortskoordinate £ und einer dimensionslosen Zeitzablung 


s=wt (11) 


1 F, Weidenhammer, Z. angew. Math. Mech., 32 (1952) Heft 8/9. 
2 L. Collatz, Numerische Behandlung von Differentialgleichungen, 5. 334. Berlin 1951. 
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werden auch die Langs- und Querverschiebungszeitfunktionen in (7) und (8) durch 
USSU, WWI (12) 


dimensionslos eingefiihrt. Durch 
Dee 
WO =| egal (13) 


wird weiter eine dimensionslose ,, Langsanregungsfrequenz‘‘ x in Einheiten f = jE J/ul4 der Biege- 
schwingungsfrequenzen gemessen. Die i-te Biegeeigenfrequenz @,; wird entsprechend durch 


oO; = Hit (14) 


bezeichnet. In der pulsierenden Axialkraft P(t) = P) + P,cos s werden die ruhende Vorlast 
P, und die Amplitude der pulsierenden Last P, in bezug auf die (kleinste) Eulersche Knicklast 
P, des Stabes fiir gleiche Lagerung gemessen: 
bale suede 
ope ae (15) 
Bei technischen Aufgabenstellungen ist ¢ und ¢ iiblicherweise klein gegen 1. Fiir die statische 
Vorlast P, wird zweckmaBig auch noch der dimensionslose Parameter 

ene? 


Dh A Ey == 0. (16) 


€o 


verwendet. Auch die Dampfungskonstante / kann durch das mit der jeweiligen t-ten Biegeeigen- 
schwingung gemessene logarithmische Dekrement dimensionslos ausgedriickt werden: 


asia 
Der Schlankheitsgrad geht ebenfalls in die Rechnung ein und wird wie iiblich mit 
A= 
VJ/F 


bezeichnet. 


Den Ansatz (7) fiir die Langsschwingung u (&, s) kann man eingliedrig und besonders einfach 
wahlen, wenn man voraussetzt, daB die Lingsanregungsfrequenz » weit unterhalb der ersten 
Langseigenfrequenz bleibt. Die Langsresonanz wird dann durch eine solche Rechnung nicht mit 
erfaBt werden kénnen. Unter dieser Einschrankung wird sich im Stab ein nennenswert 6rtlich 
verschiedener Langsspannungszustand nicht ausbilden, und man kann dann (7) hinsichtlich der 
Ortsvariablen einfach linear ansetzen. In dieser Naherung hangt die Langsdehnung u;/I und da- 
mit die Langskraft allein von der Zeit und nicht vom Ort ab. Auf die strenge Erfiillung der 
komplizierten nichtlinearen und zeitabhangigen Randbedingung (5’’) wird man verzichten und 
deren naherungsweise Erfiillung von der Lisung des Variationsproblemes (3) fordern. Der Ansatz 


w(E,s) = 73 (: _ a U(s) (18) 


erfiillt also nur die Randbedingung (5’) und nicht (5’’). Der Faktor V3 ist gemaB (9) als Nor- 
mierungskonstante zugefiigt worden. Die Zeitfunktion U(s) wird spater als Lésung einer 
Variationsaufgabe zu berechnen sein. 

Als Ortsfunktionen X; fiir den Ansatz (8) bieten sich die Eigenfunktionen (Eigenschwingungs- 
formen) des statisch gedriickten und ungedampft schwingenden Stabes an. Sie sind also aus den 
Lésungen der Schwingungsgleichung des gedriickten Stabes 


Weeee +2 0? Wee +x? w,, = 0 (19) 
mit den Randbedingungen (6’) zu gewinnen. Der Produktansatz w = X (£)T(s) fithrt auf das 
Kigenwertproblem 

AVE 2 a (20) 
mit den Randbedingungen 


X(t | = X} (4 a ==) (20) 


und eine hier nicht interessierende Schwingungsgleichung. In (19) und (20) sind die Abkiirzungen 
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(10) bis (16) verwandt und ein Akzent bezeichnet Differentiation nach &. Die orthogonalen und 
durch (10) normierten Eigenfunktionen lassen sich leicht angeben. Mit ihnen sind im weiteren 
einige Zahlenwerte zu berechnen: 


1/2 


= | ®) 
1p 
ijnr= | XUX/X, Xp de, (22) 
==} 
1/2 " 
rinse | XE Xp y 


Zwischen (21) und (23) 1aBt sich aus (20) durch geeignete Integrationen iiber die Stablange die 
Beziehung 
a Fi ; (Litare. a7. 
Wij — 20° W;; = €; 5%) mit j= \0 fir i Aj (24) 
ableiten. Damit ist auch der Ansatz (8) festgelegt. Er erfiillt die Randbedingungen fiir alle Zeit. 
Nach diesen Vorbereitungen kénnen (7) und (8) in (3) eingesetzt werden. Nach Ausfithrung 


der Integrationen iiber € bleibt dann eine Variationsaufgabe fiir die Zeitfunktionen U und W;, 
allein ibrig: 


ss N N N 
leit) BaaF an ee ae oe 
a | (eam tyae > YW; +S 2s ee i he > “ij WW; + 


ij=l i,j, ky l= i,j=l 


= N 89 N 
P3207 = 1 age 7a “ —.9 | | 2% a AY 
+ P(s) cae Ae Da W: i a mA ELA deta 0 
Hierin bezeichnen Punkte die Ableitung nach der dimensionslosen Zeit s. Die hierzu gehérenden 
Eulerschen Gleichungen sind ein inhomogenes System nichtlinearer gewéhnlicher Differential- 


gleichungen: 


a. ai V8e SN = 13202 
PU" + 32 U4 Pe Adee tegieaar aes (25) 
a 
Sg 2 Wadena vel ee eee 
Wi + Dinix Wi+ > 18; W; +732 U > eG rata > tiiniVjiW,W,=0 (26) 
" viaes Poa jr ky b=1 
GSlk Ql Nx 


Die Gleichungen (25), (26) enthalten von den speziellen Stababmessungen nur noch den Schlank- 
heitsgrad A, die Eulersche Knicklast Py, die Dampfungskonstante in den #2; und die Zahlen- 
faktoren (21), (22), (23), die durch Integrationen aus den EKigenfunktionen X; gebildet sind. Um 
die Lésungen der Gleichungen (25), (26) abhangig von Amplitude und Frequenz der harmonisch 
pulsierenden Langslast P(s) untersuchen zu kénnen, miissen die Zahlenwerte (21), (22), (23) jetzt 
angegeben werden. 

Zuvor soll jedoch darauf hingewiesen werden, da die Bewegungsgleichungen (25) und (26) 
nicht nur fiir den speziell freiaufgestiitzt gelagerten Stab gelten. Vermeidet man in (16) die Be- 
zugnahme auf die Eulersche Knicklast und rechnet dementsprechend mit dem dimensionslosen 
Parameter 2 a? gemaB (16), so gilt das Eigenwertproblem (20) mit entsprechend abgeanderten 
Randbedingungen auch fiir andere Lagerungsarten. Z. B. ist fiir den beidseitig eingespannten 


Stab nur X;' | -- e durch Xj {-+ ) zu ersetzen. Hingegen bleibt die EKigenwertgleichung (24) 


bestehen. Allerdings haben die Koeffizienten (21), (22) und (23) nicht mehr so einfache Eigen- 
schaften wie im vorliegenden Falle. Insbesondere reduziert sich die Matrix der Koeffizienten y; ; 
nicht mehr auf eine Diagonalmatrix (vgl. (29)), was weitreichende Konsequenzen fiir die Behand- 
lung dieses allgemeineren Stabproblemes hat. ; . 

Fir den freiaufgestiitzten Stab dagegen wird das EKigenwertproblem (20), (20’) zur Bestim- 
mung der Koordinatenfunktionen besonders einfach. Man findet 


X, = Zein in( £ + eae ate (27) 
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und die Eigenwerte (Biegeeigenfrequenzen) 
n= Pat Yl—el?, (28) 
wenn 2? mit seinem Wert (16) eingesetzt wird. Auch die Zahlenwerte (21) und (23) lassen sich 
leicht ganz allgemein angeben: 
ypy=O GA). Pura, (29) 
2,=0 (#j), xt =itat. (30) 
Zur Berechnung der 7;;;; hat man nach (22) ebenfalls nur bestimmte Integrale auszuwerten. 


Fiir einen zweigliedrigen Ansatz (8) mit den ersten heiden Eigenfunktionen (27) braucht man 
z. B. nur die Zahlenwerte 


3 meal 3 
Manis eS, Vaan ge ak oe None = At 28, (31) 
da die Yo. = Yoo: = 0 sind. Allgemein gilt 
4 
Hiv = FDO, (32) 


was spater gebraucht werden wird. Damit sind alle Gleichungskonstanten in (25) und (26) all- 
gemein gefunden. 


4, Beschreibung der Biegeschwingungen durch erginzte Mathieu-Gleichung. a) Auf stellung 
der Bewegungsgleichung. Einen ersten Einblick in die durch die Gleichungen (25), (26) 
beschriebenen Stabbewegungen kann man gewinnen, wenn man sich auf einen eingliedrigen Ansatz 
(8) fiir die Querverschiebung beschrankt. Aus den dann verbleibenden beiden Gleichungen kann 
man in mechanisch sinnvoller Naherung eine einzige Differentialgleichung fiir die insbesondere 


interessierende Verschiebung W; herleiten, wenn das Tragheitsglied x2 U" gegen das Riickstell- 


glied 322U in (25) vernachlassigt werden darf. Dann wird also die Langsbewegung etwa statisch 
betrachtet, und es mufB auf dieser Stufe der Naherung die Langsanregungsfrequenz x sehr klein 


gegeniiber der in (25) auftretenden ersten Langseigenfrequenz x; ~ y3a sein. Dies ist offenbar 
der Fall, wie man aus dem Vergleich der Vorzahlen x” und 3/? fiir in Betracht kommende Werte 
der Anregungsfrequenz % ersehen kann. Von den Bereichen erster Art sind ja diejenigen der 
ersten Ordnung (p = 1) besonders interessant und fiir technische Stabe ist die verlangte Un- 
gleichung 
<— 2%; a Mes [9 
xo $ = <<y3a 

auch schon fiir p = 1 sicher erfiillt. Es ist also eine sinnvolle Naherung, U aus (25) rein statisch 
zu berechnen und in (26) einzutragen. Daher findet man unter Verwendung von (24) fiir einen 
eingliedrigen Ansatz mit der i-ten Eigenfunktion: 

P 
Oe, 


In die nunmehr homogene Gleichung (33) ist wieder die urspriingliche Verschiebungsfunktion W; 
mit der Dimension einer Lange eingefiihrt worden, um fiir die Zahlenrechnung unbequeme Fak- 
toren /? zu vermeiden. An diese erginzte Mathieu-Gleichung kniipfen alle weiteren Rechnungen 
an. Sie kann als Bewegungsgleichung eines gedimpften Einmassenschwingers mit periodisch 
veranderlicher Federsteife aufgefaBt werden. Dabei ist jedoch die Federkennlinie iiberlinear, und 
zwar hat das kubische Glied gema8 (32) einen stets positiven Vorfaktor und entspricht so formal 
einer mit wachsendem Ausschlag harter werdenden Feder. Man wird daher endlich bleibende 
Schwingungen erwarten diirfen. 

Die nichtlineare Differentialgleichung (33) ist in der L.teratur schon gelegentlich aufgetreten. 
Thre allgemeine Lésung ist nicht bekannt. Hingegen sind spezielle periodische Lésungen von 
N. W. McLachlan! mit einem zweigliedrigen Fourieransatz berechnet worden. Experimente 
scheinen ebenfalls in Einzelfallen fiir die Existenz dieser Lésungen gesprochen zu haben. So 
fanden I. Utida und K. Sezawa? bei experimentellen Untersuchungen eines axial pulsierend 


* IN. W. McLachlan, Ordinary non-linear Differential Equations in engineering and physical Sciences, 
S. 122. Oxford 1950. 

* I. Utida u. K, Sezawa, Rep. aeronaut. Res. Inst., Tokyo Imp. Univ. 15 (1940), S. 193. Leider ist das 
dort beschriebene Experiment nur sehr schwer theoretisch auszuwerten. Fiir die entsprechenden Saiten- 


“2? We + S O:ninWi + xiW, — en ®y;; cos sW, + (yiiti pi) Wi= 0. (33) 


XX. Band 1952. Weidenhammer: Nichtlineare Biegeschwingungen usw. 321 


belasteten Stabes mit Endmasse ebenfalls periodische Querschwingungen mit Frequenz-Ampli- 
tuden-Abhangigkeiten, wie sie nur bei nichtlinearen Schwingungen mdglich sind. Da auch 
N. Minorsky! eine stationire Lésung fand, liegt es nahe, nach stabilen periodischen Quer- 
schwingungen zu suchen, wenn die gestreckte Gleichgewichtslage des Stabes durch Anregungs- 
frequenzen in der Nahe von x, = 2x;/p instabil geworden ist. Im weiteren wird zunichst das 
Dampfungsglied in (33) vernachlassigt. Die Berechnung der Wirkung dampfender Krafte wird 


in Abschn. 3e in gewissem Umfange nachgeholt. 

b) D as Verzweigungsproblem. Es soll jetzt durch eine einfache Anwendung der Inte- 
gralgleichungstheorie von E. Schmidt? gezeigt werden, daB es zur gestreckten Stablage fiir Lings- 
schwingungen kleiner Amplitude benachbarte kleine Querschwingungen gibt, die halb- oder ganz- 
periodisch in bezug auf die Lingsanregung sind. Es ist zweckmaBig, hier die Bewegungsgleichung 


(33) mit den Abkiirzungen 


) axa y | tiii Wii) (34) 
und 
v(s) =€ ae COS S = € COS S (35) 
in der Form : 
Ww 4 et W—wW+yW=0 (36) 
“0 


zu schreiben, wobei von der Wirkung daémpfender Krafte ganz abgesehen wird. Ferner wird die 
Bewegungszeitfunktion W ohne den Index i geschrieben und der Anregungsfrequenz % werden 
spezielle Werte x, beigelegt. In (36) rithrt v(s) von der erzwungenen Liangsschwingung her und 
ist klein mit ¢ = P,/Pg. Die von der Theorie geforderte Beschrankung auf kleine Amplituden 
von v und W ist bei in der Technik vorkommenden Schwingungen daher von selbst gegeben. 
Unter diesen Voraussetzungen wird sich zeigen, daB Stabquerschwingungen nur bei den Fre- 
quenzen (1) méglich sind. Dies ist jedoch gerade der besonders interessante Fall, in dem die 
Grundbewegung instabil ist und andere stabile Bewegungsformen nicht bekannt sind. 
Zunachst werden beispielsweise nur spezielle periodische Lésungen von (36) betrachtet, welche 


durch die ,,Randbedingungen“ 
W (—x) = W(+2) =0 (37) 


festgelegt werden. Dann ist die Integralgleichung 
eae S z 
W(s) —(2) K(s, t) W(s) dt = — K(s, t) v(t) W(t) dt +y | K (s, t) W(t) dt (38) 


0 
— —7 1 


mit (36) gleichwertig, wenn 


ns n 23 <7, mMS sinm if d 
a aa a Boa (39) 


x n\2 : Be m m\2 m gerade 
aye a 


ist. Von den Lésungen der Integralgleichung (38), die den Randbedingungen (37) geniigen, sollen 
zunachst nur die in bezug auf die Zeitvariable geraden Lésungen berechnet werden. Dann ist 
es bei dem vorliegenden Problem zulassig, die zweite Summe in (39) mit den ungeraden Eigen- 
funktionen des Kernes fortzulassen, was zur Vereinfachung der Formeln schon hier geschehen 
soll. Von der Berechtigung dieser Auslassung kann man sich durch Mitfiithrung der ungeraden 
Kernfunktionen nach einiger Rechnung iiberzeugen. Man ersieht aber auch aus der Aquivalenten 
Differentialgleichung (36) unmittelbar, da sowohl gerade periodische Lésungen als auch un- 
gerade Lésungen je fiir sich allein bestehen kénnen. Die antimetrischen Eigenfunktionen des 
Kernes (39) kénnen also nicht in die Lésung eintreten und fallen auch tatsachlich in einer aus- 
fiihrlichen Rechnung infolge der Orthogonalitat der EKigenfunktionen wieder heraus. Gesucht 
sind nun Nachbarbewegungen W zur Grundlésung W = 0 und v = 0, wie sie durch (37) und (38) 


schwingungen vgl. C. V. Raman, Physic. Rev. 35 (1912), S. 449 (zitiert nach Beibl. Ann. Physik 37 (1913), 


Salone 
1 N, Minorsky, C. R. Acad. Sc. Paris 231 (1950), S. 1417. 
2 E. Schmidt, Math. Ann. 65 (1908), S. 370; vgl. auch R. Iglisch, Mh. Math. Physik 27 (1930), S. 325, 


iiber das Duffingsche Problem. 
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beschrieben werden. Nach E. Schmidt sind diejenigen Werte (x;/x9)*, die Eigenwerte der linken 


Seite von (38) und (37) sind, Verzweigungsstellen, da in ihrer Nahe mehrere Lésungen existieren - 


kénnen. Diese Verzweigungsstellen 
gung 


Hi = | Ho (n ganz) 
sind also nach (1) gerade die kritischen Frequenzen der Bereiche erster Art. Dies ist hier zwar 
nur fiir ungerade n hergeleitet worden, kann aber durch andere Lisungstypen mit anderen Rand- 
bedingungen als (37) allgemein bestatigt werden. 

Die Verzweigungslésungen selbst sollen hier nur fiir die Randbedingungen (37) und den wich- 
tigsten Bereich erster Ordnung fiir %, = 2%; aufgesucht werden. Nach E. Schmidt hat man zu 
dem Zweck K(s, t) von der zugehérigen ersten Kigenfunktion zu befreien. Nimmt man als Trans- 
formationsfunktionen p,, g, im wesentlichen die normierten Eigenfunktionen selber, so wird 


co 


E(s, 1) = K(s,1) —B@a _ $2 coi (wm ungerade) (40) 


(alin emia 


und damit 


Mit 


mu man nun die Verzweigungsgleichung als algebraische Gleichung fiir z bestimmen. Da die 
linke Seite in (41) keine Null-Lésung mehr hat, 1aBt sich zur Bestimmung von W(s) eine (Pseudo-) 
Resolvente R(s, t) aus } 


7% 


Ris) Els.) —1 | Els) Re Hde=0 (43) 


—Tt 


als 


(n ungerade) 


und durch 


auflésen. Im einzelnen liefert diese Scheinauflésung 
Ss 
W(s) = zcos — — | G(s, t) v(t) W(t) dt +» | G(s, t) W3(t) dt, (44) 
wenn mit dem Kroneckersymbol 


sft fir n=1 
“ol \0) fittest” 


* Vel. z. B. W. Schmeidler, Integralgleichungen, Band I, S. 286. Leipzig 1950. 
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G(s, t) zur Abkiirzung fiir 


——s 1 1 = . ns ent/s 
G(s, t) = K(s, t) + 7 i R(s, t) K (ct, t) dr =- ay Cae cde 27! (i ungerade) 


It WAN] 2 
e) = oennres 


yesetzt wird. Von W, = z cos, ausgehend, laBt sich W(s) fiir hinreichend kleine Maximalbetrage 


von v und W durch schrittweise Annaherung finden, wofiir sich ein Konvergenzbeweis bei 

. . 1 . . . . . : 
L. Lichtenstein | findet. Die Verzweigungsgleichungen findet man nach E. Schmidt am einfachsten 
durch folgendes Iterationsverfahren. Die erste Naherung 


s 
W, = 2008 — 


werde mit noch unbestimmtem z in (44) eingefiihrt und so eine zweite Naherung bestimmt: 


BLA Er A - 


Ws) ~ cos | — | G(s, t) Ecos tzcos — dt +y G(s, t) 2 cos 5 Ae 


€ 


ae —n 


Hiervon werden nur die quadratischen Glieder €z, ¢2 und @ beibehalten und damit W, gebildet: 


V2 a co agele Sal 3 
3 as 7 008 GS: 
Mit der Summe von W, und W, bildet man ebenso eine nachste Naherung, von der man nur 
kubische Glieder beibehalt. Dann ist nach diesem Schritt genahert: 


W(s) ~ cos 5 E 2ez- = é7z +3 y2z3! + cos sdlimareacse 


worin alle weiteren Summanden fiir die Aufstellung der Verzweigungsgleichung nicht interessieren. 


Diese ergibt sich aus (42) bis auf héhere Glieder zu 
= Oe 
hie er LE char eos ieee © 2s (45) 


Fiir geniigend kleine ¢ und z gibt es genau drei reelle Wurzeln und damit auch ebenso viele 
Lésungen W(s), die die Periodizitatsbedingungen (37) erfiillen. Sie gehéren zu den Werten 


oe a. 
Ui, me +0(éye) , Za = —%- (46) 


AuBer der trivialen Grundbewegung W = 0 sind demnach noch die beiden periodischen Nachbar- 
bewegungen 

W(s) = + | 5,008 5 Lo 
méglich, worin das Doppelvorzeichen selbstverstandlich ist, da die Bewegungsgleichung (36) nur 
angerade Funktionen von W enthalt und homogen ist. 

Auf ganz entsprechende Weise kénnen die Verzweigungsgleichungen fiir Bereiche hoherer 
Ordnungszahlen und fiir ungerade Lésungsfunktionen aufgestellt werden. Wenn jedoch die be- 
rachteten Querschwingungen um eine verschobene Mittellage erfolgen und demnach das Zeit- 
ntegral { W(s) ds iiber eine Periode nicht verschwindet, ist der Lésungsweg tiber die Integral- 
sleichung (38) nicht ohne weiteres gangbar ?. Da man aber schon bei Bereichen zweiter Ordnung 
wegen der entsprechenden Eigenschaften der Mathieuschen Funktion ce, derartige verschobene 
Schwingungen zu erwarten hat, ist die Integralgleichungstheorie in dieser einfachen Form fiir die 
weitere Untersuchung nicht geeignet. Da sich jedoch durch vorstehende allgemeine Erérterungen 
las Problem als Verzweigungsproblem ergeben hat, und insbesondere auch die Natur der Ver- 
“weigungslésungen bekannt ist, werden im weiteren die periodischen Lésungen einfacher mit der 
Stérungsrechnung bestimmt, wobei man leider auf die von der Integralgleichungstheorie ge- 
potenen Existenz- und Konvergenzsicherheiten verzichten mu8. 


1 [., Lichtenstein, Vorlesungen iiber einige Klassen nichtlinearer Integralgleichungen usw. S. 6. Berlin 


1936. 
2 Dies hangt damit zusammen, dafB nur eine Greensche Funktion im weiteren Sinne existiert. Vgl. 


.B. Ph. Frank u. R. v. Mises, Differential- und Integralgleichungen, 8.570. 2. Aufl. Braunschweig 1930. 
23 
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c) Periodische Lésungen und ihre Frequenz-Amplituden- Gleichungen. 
Zur einfachen Durchfithrung einer Stérungsrechnung werden in die Bewegungsgleichung (33) des 


Stabes ohne dimpfende Krafte folgende Abkiirzungen eingefithrt : 


a (47) 
ue 
2)  %€) 48 
MOSS pet (48) 
F (5:33 —Vii) 
NS) ; 49 
i se (49) 


Dann lautet die Differentialgleichung (33) allgemein mit den hier zweckmaBigen Parametern: 


g@W; +(1—écoss)W, +Tew7=0. (50) 


Durch (49) ist also in die Differentialgleichung (50) der kleine Parameter ¢ zusatzlich auch als 
Faktor bei dem kubischen Glied eingefithrt worden '. Fiir die erste Eigenfrequenz wiirde man 
zufolge (29), (30) und (32) mit 1 = 1 die speziellen Werte 


Pg eee AT’ | 
ere, (47’) 
Pee 48’ 
tS Geet (48) 
F ; 
= 49 


finden. In der Technik sind die statische Vorlast und die pulsierende Last im allgemeinen sehr 
klein im Vergleich mit der Eulerschen Knicklast. Es ist daher ¢ ~ ¢ <1 und im weiteren wird 


daher stets kiirzer ¢ statt @ geschrieben. In der Gleichung (50) wird der Index i zur Bezeichnung — 


der i-ten Bewegungszeitfunktion (des Ansatzes (8)) entbehrlich, da im weiteren nur noch mit 
einem jeweils eingliedrigen Ansatz gearbeitet wird. 

Die Lisung W und die dividierte Anregungsfrequenz (48) werden nun nach dem kleinen Para- 
meter € entwickelt: 


RY 2s Ss Ee" Wan > (51) 
m=0 
iS ae (52) 


worin nun der Index m bei den Funktionen W,,, zur Kennzeichnung des Koeffizienten der Potenz 
e™ benutzt wird. Fiihrt man nun die Entwicklungen (51) und (52) ein, so lassen sich die Funk- 
tionen rekursibel als periodische Funktionen aus den linearen Differentialgleichungen 


gwrtihk=o, (50°) 


gg Wi + W,=—2q49, Wo +coss W,—TW3, (50%) 


gg Wz + Wi= —[2 904 +47] Wo —2 qq, Wi + cossW,— 3D W2W,,..., (50%) 


berechnen, wenn man die Entwicklungskoeffizienten q,, jeweils geeignet bestimmt. Um das 
Rekursionsverfahren eindeutig zu machen, werde festgesetzt, da® sich ohne Riicksicht auf An- 


fangsbedingungen fiir [° = 0 die (nichtnormierten) Mathieuschen Funktionen ce,, oder se, er- 
geben sollen. Es wird sich zeigen, daB die Frequenzgleichungen dann direkt in die entsprechen- | 


den Eigenwertgleichungen der Mathieuschen Funktionen iibergehen. 


x) ce;- Typ. Nach den Ergebnissen des Abschn. 4b gibt es halbperiodische Lisungen in der 
Nahe der Anregungsfrequenzen %) = 2%;, die jetzt berechnet werden sollen. Mit q) = 2 kommt 
von den Lésungen der Gleichung (50°) hier nur 


W, = Acos = 


in Betracht. Fir = 0 und ¢ = 0 ist dies wie verlangt die Funktion ce,, wenn speziell A = 1 | 


gewahlt wird. Die ,,Amplitude“‘ A muB jedoch hier allgemein in die Rechnung aufgenommen 


' Vel. z. B. das entsprechende Vorgehen beim Duffingschen Problem in J. J. Stoker, Nonlinear Vi- | 


brations, S. 101. New York 1950. 
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werden, um die bei nichtlinearen Schwingungen vorhandene Amplituden-Frequenz-Abhangig- 
keit aufzudecken. Weiter fithrt nun (501) auf 


f A 3 ¢ 
4W; + W, = q, Acos a +. 3 (cos 35 -- cos =) —]"A8 ia cos = s+ Z cos 7 
und liefert 
1 3 
g=—>4 ra, (54) 
A PAs 3 
We =(— 5 +g 00 (55) 
So fortfahrend findet man nach einiger Rechnung 
s A IA® 3 
W = Acoss +e|— 7 a. 35] 0% 5 5 +0), (56) 
lel aa: A PRESS 116 
q=2 + ( ol pate +( a lcs rag FAS) ot + OC) (57) 


In diesen Entwicklungen ist A? von der Gré®enordnung ¢ angenommen. Diese Annahme diirfte 
wegen der Kleinheit von F/J zufolge (49’) fiir eine technische Aufgabenstellung zutreffen. Der 
Koeffizient des linearen Gliedes in der Frequenzgleichung (57) ist bereits von N. W. McLachlan } 
auf anderem Wege gefunden worden. Insbesondere hat N. W. McLachlan auch bemerkt, daB 
fiir einen Spezialfall die strenge Lésung von (50) leicht anzugeben ist. Diese findet man aber auch 
aus der Reihenentwicklung (56) fiir eine spezielle Amplitude mit: 


W=+ |) pcos = 


und der zugehérigen Frequenzgleichung : 


q=2te— Peto). 


Dieser Sachverhalt rechtfertigt ein gewisses Vertrauen zur Liésung (56) und deren Frequenz- 
gleichung (57). 

Ganz entsprechende Rechnungen wie die vorstehende fiihrex auf die durch die Ausfithrungen 
von Abschn. 4b nahegelegten Lisungstypen. Die Ergebnisse werden hier nur fiir Bereiche erster 
und zweiter Ordnung ausfithrlich angegeben, da Bereiche héherer Ordnung kaum technische Be- 
deutung haben diirften. 


fb) se,-Typ. Hier gilt 


TA? Soa 
W=Asing +e! : 19/8 8 +O). (58) 
1 3 1 i 15 
q=24|e+glAlet |g gl Mig? +000). (59) 
y) ce,-Typ. Hier gilt 
W = Acoss +e[5 — coeds + yp cosds| +000), (60) 
5 15 
q=1+5P4te +| = sep P24] & + 0168). (61) 


Aus Gleichung (60) liest man ab, da die Schwingung um ¢ A/2 aus der Mittellage verschoben ist. 
6) seg-Typ. Hier gilt 


3 
W = Asins +e[—{ sin2s— 59 sin 3s] +(e), (62) 
3 iit la 
Sit PA + [agg 7At er 0), (63) 
q 8 24 256 | 


1 Vel. FuBnote 1 auf S. 320, In den dort angegebenen Rechnungen werden allerdings andere Parameter 


benutzt. 


23% 
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é) Cen- Typ, Sem-Typ allgemein fiir m > 2. Hier gilt: 
ae i els | m?A ele | 1)s mA san 1) | PAscos 3, + 0(e2), 


sin 2  8(m—1)sin\ 2 S732 sin 2 


2 ‘A? m. 15 . 
ae an A ~: oo 8 (m? — 1) 128 Tt oa A 
Hierin gelten die cos-Funktionen und das obere Vorzeichen fiir den ce,,-Schwingungstyp, die 
sin-Funktionen und das untere Vorzeichen fiir den entsprechenden se,,-Typ. Bis auf Potenzen & 
gilt fiir beide Typen die gleiche Frequenz-Amplituden-Beziehung. — 
Im weiteren mu nun noch mit der Methode der kleinen Schwingungen die wichtige Frage 
untersucht werden, unter welchen Bedingungen die Lésungen (56), (58), (60) und (62) stabil sind. 


d) Stabilitat periodischer Lésungen. Zur Untersuchung der Stabilitat der berech- 
neten periodischen Lisungen iiberlagert man diesen kleine Nachbarbewegungen w,) und priift 
deren Verhalten fiir unbegrenzte Zeit. Geht man dementsprechend mit W + wy in (50) ein und 
linearisiert in bezug auf wy, so erhalt man eine Hillsche Gleichung fiir wo: 


gw +(1—ecoss +3I'e W*)w = 0. (64) 

Zur Entscheidung der Stabilitatsfrage hat man nun abhangig von den Parametern q undée< 1 
zu prifen, ob wy fiir alle Zeit endlich bleibt und damit W stabil ist oder nicht. 

Zunichst findet man fiir die gestreckte Stablage die bekannten Instabilitaten wieder. Setzt 


man niamlich die strenge Lisung W = 0 in (64) ein, so erhalt man eine Mathieu-Gleichung, deren 
instabile Lésungen innerhalb der von den Frequenzkurven gebildeten Bereiche 


= Dayz 1 1 2 
te, 25 ee (65) 
— 5/24 2 


foo, = 1 4 yyogen F (08) 
liegen. In diesen Frequenzbereichen ist also die gestreckte Stablage instabil, und man hat nach 
stabilen periodischen Querschwingungen zu suchen. 


x) ce,-Typ. Die Lésung (56) der Bewegungsgleichung des Stabes iiberdeckt zufolge der Fre- 
quenzgleichung (57) bei geeigneter Wahl der Amplitude A den ganzen Frequenzbereich (65). Sie 
ist daher besonders interessant, und es soll jetzt nachgewiesen werden, daB sie in diesem Bereich 
stets stabil ist. Fiihrt man W, aus (56) in (64) ein und beachtet die Frequenz-Amplituden-Ab- 
hangigkeit (57), so wird (64) wieder eine Mathieusche Gleichung, wenn man sich auf lineare 
Potenzen von € beschrankt: 


(2— peg late) ws + (1 —|>— a Fate cos We= 0'(e)x (67) 


In (67) sind alleidings die Parameter bereits so umgerechnet, daB die Mathieu-Gleichung in der 
Normalform erscheint, wodurch der Vorzeichenwechsel des letzten Summanden im Frequenz- 
parameter von (67) bedingt ist. Stellt man nun die zu (67) gehérenden Grenzkurven der In- 
stabilitatsbereiche gemaf (65) mit den Parametern der Gleichung (67) auf, so erkennt man, daB 
die Lésung (56) stets stabil ist, denn die Kurve 


qe) == 2 ae Z 


5 pie 


liegt fiir kleine endliche Amplituden A ganz auBerhalb des Instabilitatsbereiches, der von den 
beiden Grenzkurven 


1 3 
q@)=2—fe+ Pate und g(e)=24+pe—+ Ate 


gebildet wird. Uber die mechanische Bedeutung dieses bisher wichtigsten Ergebnisses wird noch 
zu sprechen sein. 

B) ce,-Typ. Zufolge der Frequenzgleichung (59) gehért die Lésung (58) nicht zu den durch 
(65) gegebenen Frequenzen, fiir welche der gerade Stab instabil wird. Die Bewegungsform (58) 


ist auBerdem nicht stabil, da ihre Frequenzkurve 


q(e) =2+5e—4 Tae 
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innerhalb des von 
vs 1 Bish 1 ZY siete 
G2) = 2 eT AMS und ie aa i 
gebildeten Bereiches liegt, wie man aus den den Gleichungen (64) und (67) entsprechenden 
Gleichungen abliest. 


Y) ce,-T yp. Da sich bei Bereichen héherer Ordnung die Grenzkurvenentwicklungen der In- 
stabilitatsgrenzen erst in den Faktoren der héheren Potenzen von ¢ unterscheiden, hat man hier 
in (64) mindestens ¢ zu beriicksichtigen. Dies wiirde mit (60) zu einer Hillschen Gleichung mit 
einer langeren Koeffizienten- Funktion und einer entsprechend weitlaufigen Stabilitats-Diskussion 
fiihren. Wenn man aber Potenzen ¢? und hihere vernachlassigt und weiter speziell qg = 1, d.h. 
etwa Bereichmitte, annimmt, so erhalt man aus (61) eine quadratische Gleichung fiir [A®. Von 
. 5 : ‘ 5 : 
ihren Lésungen ist nur /'A2 — é +0(e) brauchbar, da die andere Lésung fiir ¢ > 0 un- 
beschrankt und daher hier nicht sinnvoll ist. Es wird daher allgemein [42 von der Gré®en- 
ordnung ¢ mit der freibleibenden Vorzahl a? angesetzt: 


De Noe oat a (68) 
Mit dieser Voraussetzung vereinfacht sich (61) zu 
5 3) eoyil ue 
q=14 ae a e +.0(e) (69) 


und man kann leicht ersehen, daB durch geeignete Wahl von a der ganze Stabilititsbereich 
zweiter Ordnung (66) des Stabes zu itherdecken ist. Fiir instabile Stablagen existiert hier also 
eine ganzperiodische Lésung, deren Stabilitat jetzt leicht nachgewiesen werden kann. Mit (68) 
wird (64) hier 


ge wy + ff +5 a? e* —eécoss +5 a? cos 2 ) Wy = 0(é°) 


und die Stabilitatsgrenzen lassen sich in bekannter Weise! durch die Berechnung periodischer 
Lésungen mit Hilfe der Stérungsrechnung finden: 


is late T 
7a igs es val + 0(e). (70) 
Lea + ale 


Durch Vergleich von (69) und (70) erkennt man, dafs (69) niemals in den durch (70) gebildeten 
Bereich fallen kann. Die Bewegung (60) ist daher stets stabil. 

6) seg-Typ. Macht man wieder die spezielle Voraussetzung (68) hinsichtlich der Amplituden- 
GréBe, so vereinfacht sich die Frequenzgleichung (63) zu: 


Offenbar liegen diese Frequenzwerte niemals in den durch (66) gegebenen Bereich. Diese Be- 
wegungsform kann daher nur eintreten, wenn der Stab in seiner Ausgangslage stabil ist. Da dieser 
Fall hier nicht von Interesse ist, soll diese Schwingung nicht aufihre Stabilitat untersucht werden. 


e) Dampfungseinflu8. Um den Einflu8 der bisher vernachlassigten dampfenden Krafte 
wenigstens genihert erfassen zu kénnen, sollen die Untersuchungen des vorstehenden Abschn. 4 
durch Ejinbeziehung geschwindigkeitsproportionaler Dampfungskrafte erganzt werden. 

Die giundlegende Naherungsgleichung (50) fiir die Querbewegung Jautet dann: 


@W; +D;~ W;+(1—écos s) W, + TEW? =0. (71) 
Hierin sind die Abkiirzungen (47), (48), (49) und auBerdem 
De 
mT 


zweckmafig, worin #; das mit der i-ten Eigenfrequenz gemessene logarithmische Dampfungs- 
dekrement bedeutet. Da #; in der Technik sicher kleiner als 0,1 sein diirfte, ist daher auch D;< 1. 


1 Vel. z. B. N. W. McLachlan, Theory and Application of Mathieu-Functions, S. 140. Oxford 1947. 
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Es soll im weiteren nur noch fiir die speziellen Anregungsfrequenzen x) = 2%; und %) = x; unter- 
sucht werden, wie groB der Schwellwert der Anregungsamplitude mindestens sein mufs, um eine 


periodische Querbewegung anzufachen. 
X) % = 2x;. LaBt man in (71) den Index i und ~ fort, so wird hier speziell 


4W +2 DW + (1—ecoss)W+TeW? =0 (71’) 
und im AnschluB an N. W. McLachlan! kann man durch den zweigliedrigen Fourieransatz 
W = A, cos = + B, sin eine Naherungsaussage fiir die ,, Amplitude“ A = ) 42+ B? finden. Aus 


den beiden Koeffizientengleichungen fiir cos s/2 und sin s/2 


DB, ee +2743 eel Bid Bia 0% 
—DA, +24 42 Pe AiB, +7 Pe 0 
1aBt sich namlich A? berechnen: 
2 et she, Di 
A lee ae 


Aus diesem von N. W. McLachlan angegebenen Resultat erkennt man, da eine Bewegung 
aus der Ruhelage W = 0 heraus nur dann erfolgen kann, wenn der Schwellwert der Anregungs- 
amplitude 

Emin = 2D 
iiberschritten wird. 

B) x) =x; Zur Berechnung des entsprechenden Schwellwertes fiir die kritische Anregungs- 
frequenz des Bereiches zweiter Ordnung hat man in 

W- + DW + (1—ecoss)W+TeW?=0 (71) 
mit dem dreigliedrigen Ansatz 
W = A, + A,coss + B,sins 
einzugehen. Hier muB man bereits in erster Naiherung drei Fourierglieder ansetzen, da zufolge 
(60) die Querschwingung um eine verschobene Mittellage erfolgt. Aus den Koeffizientengleichungen 
fiir 1, cos s und sin s 


=A 3 
Ay eA te Aya) 08 
DB, —2 4, +4 -Te A, (4? + Bi) +31 42 A, = 0, 


—DA, ++ Te B,(A? + BY) + 3T'e APB, =0 


laBt sich zwar nicht mehr so einfach die ,,Amplitude“‘ A, aber doch der ,,Phasentangens‘‘ 
_ B, 

tg @P a= ‘A, 
genahert berechnen. Aus der zweiten und dritten Gleichung findet man namlich ohne Vernach- 
lassigungen eine quadratische Gleichung fiir Dtg @ 


(Dig do —* Digd 4+ D®=0. 
1 


Da der Phasenwinkel reell ist, muB 
Ag Bs. 
Tene Dea 0 
sein. Um diese Bedingung auswerten zu kénnen, braucht man eine Naherung fiir A,/A,. Eine 
solche kann man aus der ersten Gleichung durch Vernachlassigung der dritten Potenzen der 
Fourierkoeffizienten gegen erste Potenzen finden. Dies fiihrt auf: 


A 
A; Be : 


1 Vel. FuBnote 1 yon S. 320. 
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was sich auch in (60) fiir den ungedimpften Stab auf anderem Wege ergeben hatte. Demnach 


ist der Schwellwert anzugeben: ¢,,;, = 2 yD. Zur Anregung einer Querbewegung in dem Bereich 
zweiter Ordnung ist also eine gréBere Amplitude nétig als im Bereich erster Ordnung, was nach 
den bisherigen Aussagen der Stabilititstheorie fiir gedampft schwingende Stabe zu vermuten 
war 1. 

f) Ergebnisse. Die wichtigsten Ergebnisse des Abschn. 4 sollen nun noch in den urspriing- 
lichen mechanischen Groen angegeben werden. 

Fir Anregungsfrequenzen ~ in der Nahe von , = 2@;/p und weit unterhalb der ersten Langs- 
resonanzfrequenz lieBen sich stabile periodische Querschwingungen finden. Da bei den genannten 
Frequenzen die gerade Ausgangslage des Stabes mit reinen Langsschwingungen instabil ist, kann 
der Stab aus der gestreckten Lage in Querschwingungen endlicher Amplitude iiberspringen. Da- 
bei wurde der erste Fourierkoeffizient der Bewegungsgleichung seiner iiberwiegenden GréBe wegen 


als ,,Amplitude‘‘ bezeichnet. Diese Amplituden lassen sich fiir den freiaufgestiitzten Stab ohne 


dampfende Krafte leicht allgemein angeben. Wenn die Anregungsfrequenz genau den Wert 
Wy = 2; annimmt, findet man aus der Frequenzgleichung (56) in erster Naherung [°A? = 2/3 


oder mit den Abkitrzungen (47), (49) und den Zahlenwerten (29), (30), (31) und (32): 


[- 3 efi yh P. . z 
A=LV sale) 3, (i ganz, P,< Pr), 2) 


worin t von der Ordnungszahl der i-ten Biegeeigenfrequenz herrithrt. Ganz entsprechend ergibt 


sich aus (60) und (68) fiir m) = ; 


(73) 


Vergleicht man (72) und (73), so erkennt man, dafs die Amplituden vom Wert Null auf einen Be- 
trag von der Gré®e des Tragheitsradius springen. Diese Spriinge sind in den Bereichen erster 


Oy Se +] a he ad hye P, Pi< a 


972 Py— P,)?) F Pr ( ee Pye Py 


Ordnung proportional }/ P,/P; und damit gréfer als in den Bereichen zweiter Ordnung, wo sie 
nur proportional P,/P, sind. Die Einbeziehung einer geschwindigkeitsproportionalen Dampfung 
in die Rechnung ergibt daher plausiblerweise auch, daB in Bereichen zweiter Ordnung der Schwell- 


wert zur Anregung der Amplitudenspriinge um einen Faktor von der Gré®enordnung VPi/Pr 
gréBer sein mu als in den Bereichen erster Art. — Beiden Fallen (72), (73) ist gemeinsam, daB 
die Amplitudenspriinge mit wachsenden Ordnungszahlen der zugehérigen Biegeeigenfrequenzen 
wie 1/1? abnehmen, und demzufolge die Grundeigenfrequenz die gréBte Bedeutung hat. Alle diese 
Ergebnisse passen sich gut in das Bild ein, das man sich nach den Ergebnissen der bisherigen 
Stabilitatstheorie von dem Mechanismus der Instabilitat machen muBte. 


5. Zusammenfassung. Die vorliegende Arbeit versucht, die Aussagen der ,,Allgemeinen Theorie 
der Stabilitat erzwungener Schwingungen elastischer Kérper“‘ von E. Mettler ? im Hinblick auf 
die Nachbarbewegungen instabiler Grundbewegungen zu erganzen. Als Modell fiir die ganze 
Problemklasse wird wegen seiner Bedeutung in der Technik der axial pulsierend belastete Stab 
benutzt. Die Stabilitatstheorie sagt fiir den freiaufgestiitzt gelagerten Stab insbesondere aus, 
daB in der Nahe kritischer Anregungsfrequenzen Ww, = 2 w;/p (p> 0 ganz, w; t-te Biegeeigen- 
frequenz) der harmonisch pulsierenden Langslast die gestreckte Ausgangslage des Stabes mit 
reinen Langsschwingungen instabil gegen eine Stérung in Querrichtung wird. Hingegen sind 
die dann eintretenden Querbewegungen des Stabes noch nicht zufriedenstellend berechnet wor- 
den, denn die bisher nur gefundenen, unbegrenzt exponentiell anwachsenden Schwingungen sind 
mechanisch wenig plausibel. — 

Ausgehend vom Hamiltonschen Prinzip wurde daher das Problem naherungsweise nichtlinear 
behandelt. Dabei konnte man zuniachst durch den vielfach bewahrten gemischten Ritzansatz 
die gekoppelten nichtlinearen partiellen Differentialgleichungen fiir Quer- und Langsbewegung 
des Stabes auf gleichartige, jedoch gewohnliche Differentialgleichungen fiir die Zeitfunktionen 


- yzuriickfithren. Bei Beschrankung auf den technisch allein bedeutungsvollen Fall einer An- 
 regungsfrequenz weit unterhalb der Langsresonanzfrequenz lief sich die Querbewegung des Sta- 


bes durch eine um ein kubisches Riickstellglied erganzte Mathieu-Gleichung naiherungsweise be- 


rechnen. Mit der Integralgleichungstheorie liefBen sich gerade fiir die kritischen Anregungs- 


1 Vgl. z. B. fiir den Stab: E. Mettler, Mitt. Forsch.-Anst Gutehoffnungshiitte-Konzern 8 (1950), 5. 1. 
2 Vel. FuBnote 1 von 8. 315. 
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frequenzen w, in bezug auf die Langsanregung halb- oder ganzperiodische Querschwingungen | 
auffinden. Diese nichtlinearen Schwingungen sind mit Hilfe eines Stérungsansatzes berechnet || 
worden, und gewisse von ihnen lieBen sich mit der Methode der kleinen Schwingungen als stabil i] 
nachweisen. Demnach gibt es zur instabilen Ausgangslage des Stabes benachbarte stabile Quer- | 
schwingungen. Die ,,Amplitude“‘ springt dann also von Null auf einen endlichen Wert der Gréfen- | 
ordnung des Flachentragheitsradius. Diese Spriinge sind auBerdem in Bereichen erster Ordnung | 


(d. hh. in der Nahe von wy = 2;) proportional |/P,/P,z, in Bereichen zweiter Ordnung (d. h. in 
der Nahe von Wy, = w;) dagegen nur proportional P,/Py < 1. 


Bezieht man eine geschwindigkeitsproportionale Dampfungskraft in die Rechnung ein, so | 


findet man dementsprechend, da infolge der Dampfung der Schwellwert der Lastamplitude um 


etwa den gleichen Proportionalitatsfaktor in Bereichen zweiter Ordnung gréfer sein muB als | 
in Bereichen erster Ordnung. Bereiche erster Ordnung haben also die gréBere technische Bedeu- | 


tung, da hier die Querschwingung durch geringere pulsierende Krafte anzuregen ist. 


(Aus dem Institut fiir Mathematik und Mechanik der Bergakademie Clausthal.) 
(Eingegangen am 8, April 1952.) 
Anschrift des Verfassers: Dr, Fritz Weidenhammer, Clausthal-Zellerfeld, SchulstraBe 2. 
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Ein einfaches Verfahren zur niherungsweisen Berechnung simtlicher 
Torsionseigenfrequenzen eines Stabes veranderlichen Querschnitts. 


Von H. Wittmeyer. 


1, Einleitung. Es wird cin Verfahren entwickelt, das nur eine unbestimmte Integration zur 
naherungsweisen Berechnung samtlicher Torsionseigenfrequenzen eines Stabes erfordert. Der 
Stab habe das Massentragheitsmoment J — J (x) je Lingeneinheit und die Torsionssteifigkeit 
G Jag = K = K (x). Im Intervall 0 < x <I seien K > 0, J > 0 und seien beide GréBen viermal 
stetig differentiierbar. Fiir die Drehamplitude ¢, (x) und die Frequenz f, = W,/2% der n-ten 
Kigenschwingung ergibt sich folgendes Eigenwertproblem fiir die ,,reine“‘ (d. h. nicht mit Biege- 
oder Langsschwingungen gekoppelte) Torsionsschwingung : 


d dp, ; 
are F) + oh Ton = 0. (1) 


Das entwickelte Verfahren gestattet dessen angenaherte Lésung fiir die folgenden drei ver- 
schiedenen Randbedingungen: Fest-frei, fest-fest und frei-frei. In dieser Arbeit werden nur die 
ersten Randbedingungen, d. h. 


G0) = Dig eK * |e bo, (2) 


behandelt. Die beiden anderen lassen sich mit Hilfe der ,,Zerlegungsmethode“ (Ziff. 6) auf die 
erste Aufgabe zuriickfiihren. Eine ausfiihrlichere Darstellung samtlicher hier behandelten Fragen 
findet man in der friiheren Arbeit des Verfassers!, die im folgenden mit W zitiert werden soll. 

Die in dieser Arbeit angegebene Methode erfordert einen geringeren Arbeitsaufwand als das 
Verfahren von Rayleigh- Ritz und das Iterationsverfahren, sowie Kombinationen von beiden, 
und zwar einen relativ umso geringeren, je héher die Ordnung der berechneten Frequenz ist. 
Diese Einfachheit ist dadurch erreicht, daB bei der hier entwickelten Methode auf die Berechnung 
von Schwingungsformen verzichtet ist und ebenso auf die Méglichkeit, die Eigenfrequenzen 
belie big genau zu berechnen. (Die erreichte Genauigkeit ist jedoch gut und fiir die Praxis 
im allgemeinen ausreichend.) — Die neue Methode liefert in den durchgerechneten Beispielen 
mit ,,glatt veranderlichen“‘ Massen- und Steifigkeitsverteilungen genauere Werte als eine Methode 
etwa gleichen Rechenaufwandes von P. B. Walker (die nur die Grundfrequenz liefert). Das 
liegt unter anderem daran, daB P. B. Walker den Stab durch einen masselosen mit einem Einzel- 
tragheitsmoment ersetzt, wahrend in der vorliegenden Arbeit der (nach einem anderen Gesichts- 
punkte gefundene) Ersatzstab glatt veranderliche Querschnittswerte hat. 

Die in der vorliegenden Arbeit angegebenen asymptotischen Formeln zur Berechnung der 
n-ten Eigenfrequenzen unterscheiden sich von den bisher bekannten (vgl. Collatz*, Ince‘, 
Hohenemser- Prager *) dadurch, da® sie nicht nur fiir n > oo asymptotisch genau werden, sondern 
auch fiir n = 1 richtig oder-angenahert richtig sind. 

Die in Ziff. 6 entwickelte ,,Zerlegungsmethode“‘ zur Berechnung verfeinerter Werte fiir die 
Oberfrequenzen benutzt einen Gedanken, der sich schon bei 0. Féppl und R. Grammel findet °. 
Dieser besteht darin, da man bei den héheren Eigenschwingungsformen die Knoten annehmen 
soll und dadurch die Berechnung der héheren Frequenzen auf die Berechnung der Grund- 


‘frequenzen der zwischen je einem Knoten und einem benachbarten Bauch gelegenen Teilsysteme 


zuriickfiihren soll. Eine Frequenzberechnung nach diesem Gedanken war bisher jedoch nur bei 


1H. Wittmeyer, Asimple method for the approximate computation of all the natural frequencies of 
a member with variable cross-section. Technical Notes publ. by the Division of Aeronautics of the R. Insti- 
tute of Technology, Stockholm. Nr. 13. 1951. 

2 P. B. Walker, Simple formulae for the fundamental natural frequencies of cantilevers, R& M No. 1831, 

3 L. Collatz, Eigenwertprobleme und ihre numerische Behandlung, S. 117 ff. Leipzig 1945. 

4 EF. L. Ince, Ordinary Differential Equations, Dover Publications, S. 270 ff. New York ohne Jahres- 
zahl, First American Edition. Vorwort von 1926. 

5 K, Hohenemser, Die Methoden zur angenaherten Lésung von Eigenwertproblemen in der Elasto- 
kinetik. Ergebnisse der Mathematik und ihrer Grenzgebiete, Bd. 1, H. 4. Berlin 1933. — K. Hohenemser 
und W. Prager, Dynamik der Stabwerke. Berlin 1933. 

6 Siehe bei K. Hohenemser, a. a. O. S. 62. 


332 Wittmeyer; Ein einfaches Verfahren zur naherungsweisen Berechnung usw. _ Ingenieur-Archiv 


masselosen Staben mit groBen Einzeltragheitsmomenten oder bei Staben konstanten Quer- 
schnitts praktisch durchfithrbar. Mit Hilfe der in der vorliegenden Arbeit aufgestellten neuen 
Formel zur Berechnung der Grundtorsionsfrequenz ist es nun méglich geworden, diesen Ge- 
danken auch fiir Stabe mit glatt veranderlichen Querschnittswerten zu einer praktisch brauch- 
baren Methode (Ziff. 6) auszubauen. 

Bei der Ableitung des neuen Verfahrens wird in einem entscheidenden Punkte die Stérungs- 
rechnung benutzt. Bemerkenswert ist dabei, daB sie nicht, wie sonst iiblich 1 an den Differential- 
gleichungen durchgefithrt wird, sondern an den zugehérigen Variationsproblemen. Das ist hier 
moglich, weil nicht nur das gestérte Problem selbstadjungiert ist, sondern auch das ungestérte 
Problem ebenso gewahlt ist. Diese Art der Durchfiihrung ist hier deswegen vorteilhaft, weil 
heide Variationsprobleme die gleichen (kiinstlichen) Randbedingungen haben, dagegen die zu- 
gehorigen Differentialgleichungen sich in einem Teil der Randbedingungen unterscheiden (ndm- 
lich in den natiirlichen Randbedingungen, die beim Variationsproblem ja fortgelassen werden 
kénnen). 


2. Das Verfahren. Aus den gegebenen Funktionen J = J(x) und K = K(x) bildet man (wie 
spater zu begriinden sein wird) 


L = log (eJ K) (c > 0, beliebig, konstant) (3) 


und 
x 


i= Ww (a) ‘ ye dx (w (Ne w;) (4) 
6 
und findet daraus durch Elimination des Parameters x (eventuell zeichnerisch) die Beziehung 


L = L (w), aus der man (Abb. 1) 
Ipite: 


L (0,123 wi) —L (0,877 wi) (5 
0,754 ) 
berechnet. Damit berechnet man weiter 


F, = F, (k) = 0,25 + 0,05833 k + 0,003173 k2 — 0,0001598 k? (giiltig fir —1,7 <k <7), (6) 


4 =A (k/ 
42 
L (ro) 
ap 10 = 
\ 
NG 08 
+ 
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‘ ee en . k 
0 w=OT3Ny, 7=087770;, 10=70(U/=70, =) ey 0 2 Y C Sik Bile 
Abb. 1. Graphische Ermittlung von k [GI.(5)]. Abb, 2. 7 Q i go ii i u 
Ee ee ante one pe en f: -2, Zur Ermittlung yon F, fiir die Gl. (7a, b). 
als A hat. 


bzw. liest den Wert von F; aus einer ein fitr allemal gezeichneten Kurve F, = F, (k) ab (Abb. 2). 
Die Torsionseigenfrequenzen f,, (Perioden in der Sekunde) sind dann 


1 
ia: eee 

liek US sae —2in=1 A 

nine gee eke a oe (7a, b) 


Die Werte f, werden bei stark veranderlichem L = L (w) fiir n = 2 und n = 3 zweckmabiger- 
weise noch durch das in Ziff. 6 angegebene Verfahren verbessert. 


3. Beispiel. Mit K(x)/K = J(x)/J, =(1—08 “al wird fiir ¢ = (Kg Jo) zunachst L(w) 
= § Le log (1—0,5 w/w))|. Mit (5) und (6) erhalt man k = 2,367 und F, = 0,4037. Die sich damit 


1 Velo. Collatz, a. a, O, S305 tf 
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aus (7a, b) ergebenden Frequenzwerte f,, f, und f, sind in der folgenden Zahlentafel den exakten 
Werten gegenibergestellt, die durch Umrechnung aus Rechenergebnissen einer Arbeit von 
W. Meyer zur Capellen? gewonnen sind: 


n | I | 2 3 | 4 


EV Taf 0,4018  0,8265 | 1,300 | 1,783 
LY Jol Ko fn, Naherung 0.4037  0,8012 | 1,281 1,772 
Fehler +0,5% | 31% | —1,5% | —0,6% 


Auch weitere durchgerechnete Beispiele ergaben ahnlich gute Resultate. GroBere Abweichungen 
ergaben sich nur bei f,(~ — 10%). Solche Fehler bei den Oberfrequenzen kénnen jedoch mit 
der ,,Zerlegungsmethode“ (Ziff. 6) rasch verringert werden. 
4. Ableitung der Formel fiir die Grundfrequenz. Statt der Differentialgleichung (1) mit den 
Randbedingungen (2) wird der Untersuchung das zugehérige Variationsproblem 
Somete 
| K (as) dx — a? | J ¢? dx = stationar (8) 


0 0 
mit der einzigen Randbedingung 
¢ (0) =90 (9) 


zugrundegelegt. Die zweite in (2) enthaltene Randbedingung stellt sich von selbst als natiirliche 
Randbedingung ? ein. Benutzt man? die neuen Groen 


i = i Vz dx (we = / 2 ax) ; (10) 
10) —1RF 90), ay 
js ith Ue (12) 


setzt abkiirzungsweise 
h(é) = + log (cK J) (c > 0, beliebig, konstant) (13) 


und bezeichnet Differentiationen nach € mit Strichen, so geht das Eigenwertproblem (8), (9) 
tuber in 
1 
| [(4’ —h'n)? —ArP\ dé = stationar, 17(0)=0. (14) 
0 
Es liegt nahe, den Eigenwert /, des Problems (14) dadurch angenahert zu berechnen, daB man die 
Funktion h’(é) durch eine Konstante a ersetzt. (Uber die Wahl von a vgl. weiter unten.) Dann 
hat man das neue Problem 
1 
is [(y’ — an)? —ArP| d& = stationar, (0) =9 (15) 
i 
fiir den genadherten Eigenwert 7, und die geniherte (normierte) Eigenfunktion 7, zu lésen. 
Der Ubergang von (14) zu (15) bedeutet physikalisch, daB der gegebene Stab durch einen 
anderen fest-freien mit einfacherem Steifigkeits- und Tragheitsmomentenverlauf ersetzt wird, 
dessen Eigenfrequenzen fiir n—> oo das gleiche asymptotische Verhalten zeigen wie die des ge- 
gebenen Systems. Ap 
Schon Hohenemser* hat darauf hingewiesen, daB es zur Frequenzberechnung zweckmabig ist, 
den gegebenen Stab durch einen einfacheren mit gleichem asymptotischen Frequenzverhalten 
zu ersetzen. ‘ 
Durch Lésung der zu (15) gehorigen linearen Differentialgleichung mit konstanten Koeffi- 
zienten erhalt man folgende Naherungsformel, deren Fehler im Intervall —4 <a <1 kleiner 


1 W. Mey Crpiicn, Aun, Physik 6), 8 (1931) 57297. ; 

Ott. Cae Hilbert. Methoden der Mathematischen Physik, Bd. 1, S. 179 ff., 2. Aufl. Berlin 1931. 
3 Vel. E. L. Ince, a. a. O. S. 270. 

cand. Wire O yur Sar 


334 Wittmeyer: Ein einfaches Verfahren zur naherungsweisen Berechnung usw.  Ingenieur-Archiv _ 


als 0,5% ist: 
Ve 0,25 _. 0,10132 a + 0,009575 a2 + 0,0008375 a?. (16) 
ELA 


Nun ist noch ein zweckmafiger Wert fiir a zu bestimmen. Dazu denkt man sich die Pro- 
bleme (14) und (15) mit Hilfe eines Parameters ¢ durch folgendes Problem verbunden, das fiir 
¢ = 0 in (15) und fiir ¢ = 1 in (14) ttbergeht. (Damit 7 eine eindeutige Funktion von ¢ wird, fiigt 
man eine Normierungsbedingung fiir 7 hinzu.) 


1 1 1 
i M(y, 7/3 €) dé atl eae == state, 70) 00s [ yr dete |, (17) 
0 0 0 
wobei 
M(n, 7's. 0) = (4 —an), M(m, 931) = (yf — BP (18) 
gelten soll. Der niedrigste Eigenwert von (17) sei gleich “4, = p,(e) und die zugehdérige Eigen- 
funktion 7, = 7, (&: 2). Bezeichnet man die entsprechenden Eigenwerte und normierten Kigen- 
funktionen von (14) mit /,, 7, und von (15) mit A,, 7,, so ist 
(0) =A, 4) =A (19) 


und 

mh (530) =m (8), 7m (E31) =m (E)- (20) 
Da der Naherungswert J, méglichst wenig von dem exakten Wert /, verschieden sein soll, so 
liegt es nahe, als Bestimmungsgleichung fiir a zu fordern 
Aus (21) folgt, da 7, die Eulersche Variationsgleichung von (17) mit « = (1,(0) erfullt, 


1 
TMs = 0% fir oO, erate (22) 
0 


(21) 


Nun setzt man 


M = M, = {n' —[a +e (h’ —a)|n}?. (23) 
(Uber die ZweckmaBigkeit dieser Wahl vel. Ziff. 8.) Durch Einsetzen von (23) in (22) erhalt 
man nach geeigneter Umformung, wenn man zur Abkiirzung noch 
P, =, Gf. —@m) (24) 
setzt, 
1 
[Ww P, dg 
ee (25) 
[ Pydé 
0 


Die Konstante a ist also ein gewogenes Mittel von h’(&) mit der Gewichtsfunktion P,. Ersetzt 
man in (25) h (&) durch ein Naherungspolynom vierten Grades 


h=h +m, f al + m, f =a 1 ms [ — Al +m, [ — ral : 


so ergibt sich, da die rechte Seite von (25) nur sehr schwach von a abhangt (jedenfalls im prak- 
tisch wichtigen Bereich a < 0,7), 


__ (0.877) —h (0,123) 
0,754 (25a) 


Fir a muf8 somit der Differenzenquotient der Funktion h (€) zwischen den Werten & = 0,123 
und € = 0,877 gewahlt werden. Fiihrt man zur Vereinfachung der numerischen Rechnung 
statt der GréBen &, h (€) und a die GréBen 


w=wé, L=A(loge)h, k = — 4.(2° loge) a (26) 


ein, so erhalt man das in Ziff. 2 geschilderte Verfahren zur Berechnung der Grundfrequenz f,. 


a ~m, + 0,142 m, 
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5. Ableitung der asymptotischen Formeln fiir die Oberfrequenzen. Es ist nicht zweckmabig, 
auch die héheren Frequenzen f,, nach der Methode des vorangehenden Abschnittes anzunahern, 
weil die der Gewichtsfunktion P, entsprechenden Gewichtsfunktionen P, dann im Intervall 
0 <§ < 1ihr Vorzeichen wechseln wiirden, was zu recht komplizierten Ausdriicken fiir a fiihren 


witirde. Statt dessen werden zwei andere Methoden benutzt, von denen die erste hier und die 
zweite in Ziff. 6 beschrieben wird. 


Setzt man in Verallgemeinerung der Formel fiir f, 


foe = =F, ( j | £ dx), (27) 
0 


so gilt nach Ince? folgende asymptotische Formel fiir F,,: 


2n—1 I 
Fi, fick ro (3 — i 4 ce) 
Man erkennt, daB 
1 
ors 
2n—] 1 4 
3 a ag ae eves mn cae (29) 


fiir n = | genau ist und fiir groBe n in erster Naherung gemaf (28) asymptotisch genau ist. 

Bei der praktischen Anwendung wird man in (29) fiir F, den nach (6) gefundenen Wert ein- 
setzen. Das ergibt das in Ziff. 2 dargelegte Verfahren zur Berechnung der Oberfrequenzen. 

6. Verbesserung von Niherungswerten fiir die Oberfrequenzen (Zerlegungsmethode). Beziig- 
lich des Rechenschemas sei auf W verwiesen. Hier soll nur die Ableitung des Verfahrens skizziert 
werden. 

Schneidet man einen Stab an den Bauchen der n-ten Eigenschwingungsform durch und hilt 
ihn an den Knoten dieser Schwingungsform fest, so ist die Grundfrequenz f,, ; (i =1, 2, ..., 2n—1l) 
eines jeden der so entstandenen (2n—1) fest-freien Stabe gleich der Frequenz f,, der n-ten Kigen- 
schwingungsform. Die w-Werte (4) der Knoten und Bauche seien mit w; bezeichnet. Ferner 
sei p; die Steigung einer Geraden, die in dem i-ten Intervall w;_, < w < w; in entsprechender 
Weise an die Kurve L = L (w) konstruiert ist, wie die Gerade durch A, B in Abb. 1 in dem 
Intervall 0 Sw < w, (p ist also eine mittlere Steigung von L (w) in dem i-ten Intervall). Dann 
ergibt sich mit k; = (— 1) p; (w; — w;_,) nach Ziff. 2 


F, (k; 
eeeg ee 


W.—W. 
Eine von der Lage der Knoten wenig abhiangige Darstellung von f,, ist offensichtlich 


(30) 


i t—1 


1 2n—l1 2n—1 


Z=h aes > Bit, i (g = a fi). (31) 


i=4 al 


wenn die g; noch so gewahlt werden, dah 
ees, (32) 


ist. Gleichung (32) ergibt fiir die Gewichte g; naherungsweise 


C061) p:\ae 
& = (4 — 1) (1 2 my , 


(33) 


Die Gleichungen (31) und (33) werden zur iterativen Verbesserung von fn benutzt, indem man in 
die rechten Seiten die bei einer irgendwie angenommenen Intervalleinteilung bekannten Nahe- 
1 nee ; 
rungswerte wi), p??, sowie die mit ihnen nach (30) berechneten Werte fr und schlieBlich einen © 
geschatzten Wert f statt der wahren Werte w;, p;, f;,; und f, einsetzt und so auf der linken Seite 
1 F 

von (31) den verbesserten Wert fase erhalt. ; 

Eine neue Intervalleinteilung berechnet man dann in folgender Weise: Aus der Naherungs- 
formel 


F, (ky) ~ 4 + 0,06 k; = 4 + 0,06 (— 1)" ps (10; — 105-1) (34) 


tp AC a. O..S; 270T£. 
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ergibt sich durch Auflésung nach w; — w;_, unter Beriicksichtigung von F, (k;) = fp > (Wi — w-1) 
[siehe (30)] ; 


Aw; = w; — w;_,= [4(fr — 0,06 (— 1)‘ p,)I". (35) 
Diese Gleichung wird in der Weise zur Iteration benutzt, daB man rechts statt f, und p; die 
GréBen fo und p?? einsetzt. Auf der linken Seite erhalt man neue GroBen (Aw;)” , die nach 


einer Normierung, so da ihre Summe gleich w, wird, verbesserte Gréfen (Aw,)” ergeben. Durch 
diese (A w;)” — w\) — wl), ist eine verbesserte Intervalleinteilung mit den Abszissen w = wi?) 
gegeben, mit der der nachste Iterationsschritt beginnen kann. — Da z. B. zur Berechnung von 
f, das Intervall Ow Sw, nur drei Teilintervalle hat, ist der fiir die Zerlegungsmethode er- 
forderliche Rechenaufwand nicht gro8. Fiir groBe n braucht die Methode nicht angewandt zu 
werden, da fiir diese ja die asymptotischen Formeln genau genug sind. 

Als Beis piel fir die Zerlegungsmethode ist in W die Frequenz der ersten Oberschwingung 
eines bei x = 0 einseitig gehaltenen Stabes mit den Querschnittswerten K/K, = J/Jy = 
(1 — 0,75 x/l)* berechnet worden. Als Ausgangslésung wurde eine EKinteilung des Intervalles 
0 <w <w, in drei gleiche Teile benutzt, sowie ein nach (7b) berechneter FS?-Wert mit einem 
Fehler von —10%. Das Verfahren ergab nach einem Iterationsschritt einen verbesserten Wert 


FY) mit einem Fehler von — 0,2%. 

7. Fehlerabschaitzung. Wie in W naher gezeigt ist, ist der Fehler bei der Grundfrequenz f, [bis 
auf den in (6) enthaltenen Fehler von maximal 0,5%] proportional zu Max (h’’ ())?, wobei 
h(€) nach (10) und (13) zu berechnen ist. Er kann positiv oder negativ sein. Die nach (7a, b) be- 
rechneten Werte von f, und f; weisen im allgemeinen gréBere Fehler auf als der von f,. Die nach 
der Zerlegungsmethode verbesserten Werte fiir die Oberfrequenzen sind jedoch durchweg ge- 
nauer als der Wert fiir die Grundfrequenz, und zwar 

a) weil die zu einem Teilstabe gehérige Funktion h(€) (0 <€ < 1) nicht so stark gekriimmt 
ist, wie die zum ganzen Stabe gehérige (andere) Funktion h(€) (0 <& < 1), 

b) weil f, durch die Mittelbildung (31) mit der Nebenbedingung (32) gewonnen wird. 

8. Bemerkung zur Wahl von M = M, in Gleichung (23). Die Wahl von M = M, in (23) er- 
scheint zunachst willkiirlich, da M, nur cine von den unendlich vielen M-Funktionen darstellt, 
die die Bedingung (18) erfiillen. Es laBt sich jedoch zeigen, dafi M, unter allen gemaB (18) zu- 
lassigen M-Werten dadurch ausgezeichnet ist, daB es als einziges M die folgenden Eigenschaften 
hat, auf deren Vorhandensein die Kinfachheit und die Genauigkeit der in dieser Arbeit ent- 
wickelten Methode beruht: 

a) ¢ tritt in M, explizit héchstens quadratisch auf. 

b) Die Wahl M = M, fithrt zu einer Gewichtsfunktion P, [vgl. (24)], die einfach gleich dem 
Produkt der linken Seiten der kiinstlichen Randbedingung 7, = 0 fiir € = 0 und der natiirlichen 
Randbedingung 7; —a7, = 0 fiir € = 1 ist. 

c) Der Fehler der zu M = M, gehérigen Methode verschwindet mit der zweiten oder héheren 


Potenz von Max |h’ (é)|. 
OSes 


9. SchluBbemerkungen, Das in dieser Arbeit beschriebene Verfahren wird man zweckmabiger- 
weise dann anwenden, wenn die Schwingungsformen nicht interessieren und bei den Frequenzen 
ein Fehler von wenigen Prozenten zulassig ist. (Bei den unter Zuhilfenahme der Zerlegungs- 
methode bisher durchgerechneten Beispielen war der absolute Fehler héchstens 1,8%.) Verlangt 
man gréBere Genauigkeiten als etwa 1%, so benutzt man zur Berechnung der Frequenzen 
niederer Ordnung besser eins der bekannten Verfahren, mit denen die Frequenzen mit jeder 
vorgeschriebenen Genauigkeit berechnet werden kénnen. Unter diesen Verfahren wird das von 


R. Grammel angegebene den relativ geringsten Mehraufwand an Rechnung gegeniiber dem 


neuen Verfahren haben. 

Die in der vorliegenden Arbeit entwickelte Methode 1aBt sich auBer auf die Berechnung von 
Torsionseigenfrequenzen auch auf die Berechnung der Eigenfrequenzen der Langsschwingungen 
von Stiben und der Querschwingungen von Saiten, sowie auf die Berechnung der Knicklasten 
von Staben mit verdnderlicher Langslast und veradnderlichem Querschnitt anwenden. 


(Eingegangen am 16. April 1952.) 
Anschrift des Verfassers: Dr.-Ing. H. Wittmeyer, Linképing (Schweden), Valkebogatan 14 A. 


i = R. Grammel, Ein neues Verfahren zur Loésung technischer Eigenwertprobleme. Ing.-Arch. 10 (1939), 
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Uber die Differentialgleichungen fiir die Bewegung eines idealisierten 
Kraftwagens*. 


Von B. Stiickler. 


1, Vorbemerkungen. In dieser Arbeit werden die Bewegungsgleichungen eines Kraftwagens 
aufgestellt, der auf einer Ebene rollt, ohne zu gleiten. Die Ebene kann geneigt sein, und es 
mogen irgendwelche eingepragte Krafte auf das Fahrzeug einwirken. Selbstverstindlich kommt 
man bei der Kompliziertheit eines modernen Kraftwagens nicht ohne mehr oder weniger weit- 
gehende Idealisierungen ans Ziel. Daher miissen auch als erstes diejenigen Vereinfachungen 
angegeben werden, von denen ausgegangen wird, wobei aber hier nicht erértert werden soll, 
inwieweit mit ihnen in der Praxis gerechnet werden kann: 

1. Das Fahrzeug wird wie ein System starrer Kérper behandelt. 

2. Es wird nur das reine Rollen des Fahrzeuges untersucht. Das bedeutet u.a., da&B vom 
Steuermechanismus vorausgesetzt werden muB, da durch ihn das reine Rollen in keiner Weise 
beeintrachtigt wird, was bei den meisten Kraftwagen nicht ganz der Fall ist. Auch ein Abheben 
der Rader wahrend der Bewegung soll nicht stattfinden. 

3. Die Beriihrungsflache zwischen einem Rad und der Fahrebene sei so klein, da®B von einem 
Berithrungspunkt gesprochen werden kann. Dieser soll seine Lage gegeniiber der Radachse 
nicht andern. 

4, Bei Beriicksichtigung der Relativbewegung von Motor- und Kraftiibertragungsteilen 
gegeniiber dem Fahrgestell werden nur die rotierenden Massen, wie Schwungrad, Kurbelwelle 
und die gréBeren Zahnrader, beachtet. 

5. Die Relativbewegung von Insassen gegeniiber dem Fahrgestell wird nicht beriicksichtigt ; 
diese werden gewissermafen in ihrer normalen Stellung auch als starr angesehen. 

Fiir einen Kraftwagen, der diesen Voraussetzungen genau entsprechen wiirde, werden die 
genauen allgemeinen Bewegungsgleichungen aufgestellt, wahrend die bisher in der Kraftwagen- 
Technik benutzten Gleichungen nur fiir die reine Geradeausfahrt und andere Sonderfalle gelten. 
Dabei wird so vorgegangen, daf} sich die Bewegungsgleichungen anderer Fahrzeugarten, wie des 
gewohnlichen vierradrigen Wagens und des Dreirades, als Sonderfalle ergeben. 

Die betrachteten Fahrzeuge gehéren zu den nichtholonomen Systemen mit einer endlichen 
Anzahl von Freiheitsgraden. Daher liegt es nahe, die Bewegungsgleichungen mit Hilfe der sog. 
Lagrange-Eulerschen Gleichungen aufzustellen. Es erweist sich aber, da man auf andere Weise 
erheblich Jeichter und schneller ans Ziel gelangt, naémlich durch Anwendung des Prinzips von 
D’ Alembert und des Energiesatzes fiir skleronome Systeme. Hiermit werden in Ziff. 2 die genauen 
Bewegungsgleichungen aufgestellt. 

Bei Beriicksichtigung von Motor- und Kraftiibertragungsteilen wird an die gebrauchlichste 
Antriebsart gedacht, namlich an den Antrieb durch einen Motor ither Ausgleichsgetriebe. Dabei 
wird gleich vom Vierrad-Antrieb ausgegangen, doch so, daB sich aus den aufgestellten Glei- 
chungen ohne weiteres sowohl der reine Hinterradantrieb als auch der reine Vorderradantrieb 
als Sonderfalle ergeben. Der Antrieb wird aber vom iibrigen Fahrzeug gesondert behandelt, 
so daf leicht zu iibersehen ist, was zu andern ist, wenn es sich um andere Antriebsarten handelt. 

In Ziff. 3 ergeben sich wesentliche Vereinfachungen durch Beschrankung auf den Fall kleiner 
Ausschlagwinkel der Vorderrader, der beim Kraftwagen in der normalen Fahrt meistens vorliegt. 
Hierzu wird ein Zahlenbeispiel durchgerechnet mit Zahlenwerten, wie sie bei W. Kamm und 
C. Schmid! zu finden sind. SchlieBlich wird als erste Anwendung der aufgestellten Bewegungs- 
gleichungen die verhaltnismabig einfache Frage untersucht, wie sich ein Kraftwagen bei einer 
Geradeausfahrt mit konstanter Geschwindigkeit gegenttber kleinen Stérungen verhalt. Es zeigt 
sich, da® unter den getroffenen Voraussetzungen ein Wagen mit Vorder- oder Vierrad-Antrieb 
von selbst wieder zur reinen Geradeausfahrt iibergeht, wahrend dies beim Hinterradantrieb 

* Gekiirzte Fassung der von der Technischen Universitat Berlin 1951 angenommenen Dissertation 


(Berichter Prof. Dr. W. Haack, Mitberichter Prof. W. Kucharskt). ; ; 
1 W. Kamm und C. Schmid, Das Versuchs- und MeBwesen auf dem Gebiet des Kraftfahrzeugs, Berlin 


1938. 
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nicht der Fall ist. Das beruht im wesentlichen auf dem Unterschied in der kinematischen 
Beschaffenheit der Antriebsarten. 

In Ziff. 4 wird die kriftefreie Bewegung des gewohnlichen vierrddrigen Wagens und anderer 
Fahrzeuge untersucht, deren Bewegungsgleichungen als Sonderfalle aus den Gleichungen des 
Kraftwagens hervorgehen. Dabei gelangt man zu dem wohl iiberraschenden Ergebnis, dal 
die genauen Gleichungen dieser Fahrzeuge exakt gelést werden kénnen, und zwar durch ele- 
mentare Funktionen, so da8 zur Ermittlung der Bahnkurven nur noch gewéhnliche Quadraturen 
erforderlich sind. Fir die wichtigsten Falle werden mit Zahlenbeispielen Bahn- und Schlepp- 
kurven errechnet und gezeichnet. 


2. Aufstellung der Bewegungsgleichungen. a) Kennzeichnung der Fahrzeuge. 1. Der 
Kraftwagen. In Abb. 1 ist das Fahrgestell des Kraftwagens durch das Rechteck B B, B, B, 
angedeutet, wihrend die Rader einfach 
durch die Projektion ihrer Berithrungs- 
kreise mit der Fahrebene dargestellt 
sind. Die Vorderréder seien um zur 
Fahrebene senkrechte Achsen durch B 
und B, drehbar gelagert und durch den 
Steuermechanismus miteinander ver- 
bunden (in Abb. 2 schematisch skiz- 
ziert). Die Lenkachsen haben vonein- 
ander den Abstand b’, ihr Abstand von 
der Hinterradachse sei !. Der Halb- 
messer der Vorderrader sei r, derjenige 
der Hinterrader r,;. Die Beritthrungskreise 
der Vorderrader mit der Fahrebene 
haben von den entsprechenden Lenk- 
achsen den Abstand c, derjenige des 
linken Hinterrades habe den Abstand a, 

Abb. 1. Bezeichnungen. derjenige des rechten den Abstand b 

von B B,. Gewohnlich ist beim Kraft- 

wagen b—6'-+a. Wenn aber diese Abhangigkeit zunachst nicht vorausgesetzt wird, kénnen 
gleich andere Fahrzeugarten miterfaBt werden, namlich: 

2. Der vierrddrige Wagen und der lenkbare Schlitten. Bei diesen ist b’=0, b—a. 

3. Das Dreirad. Bei diesem ist 6’ —0, b=a, c=0. 

4. Das Fahrrad oder Motorrad mit Beiwagen. Bei diesen ist b’=0, c=0, a—0, 
bA0. Allerdings sind bei diesen Fahrzeugen die Lenkachsen geneigt; auch beim Dreirad ist 
dies in der Regel der Fall. 

b) Anzahl der Freiheitsgrade und der Bedingungen. Es seien x, y die rechtwinkligen 
Koordinaten von Bin der Fahrebene. Die Langsrichtung des Fahrzeuges sei durch den Winkel 0, 
festgelegt, die Stellung der Vorderrider gegenitber der Langsrichtung durch die Winkel y und 7p. 
Der Rollwinkel des linken Vorderrades, gemessen vom Beriihrungspunkt bis zu einer Marke auf 
dem Rade, sei gy. Entsprechend sei ¢~, der Rollwinkel des linken Hinterrades, y, der Rollwinkel 
des rechten Vorderrades, p; der Rollwinkel des rechten Hinterrades. Sind diese 9 GréBen bekannt, 
so ist die Lage des Fahrgestells und der 4 Rader auf der Ebene vollstandig bestimmt. Die Drehung 
von Motor- und Kraftiibertragungsteilen um deren Achsen wird hier noch nicht beriicksichtigt, 
sondern erst im Abschnitt h). Das bedeutet, daB die in den folgenden Abschnitten aufgestellten 
Bewegungsgleichungen zundchst nur bei nichtlaufendem, ausgekuppeltem Motor gelten. Vor- 
laufig hat man also nur mit der 9 Koordinaten x, y, 01, 7. Yo, P, Py, Pa» Pz Zu tun. Wird fir diese 
die allgemeine Bezeichnung q, (vy = 1,2,..9) verwendet, ist jeder Ortsvektor T des Fahrzeuges 
eine Funktion dieser g, und eines Vektors a, der den einzelnen Fahrzeugpunkt charakterisiert 
und zeitunabhangig ist: 


T= T (G1) qa +> Go3 4). (1) 
Unter der Voraussetzung, dai die Fahrebene im Raum ruht, kommt die Zeit t in TF nicht explizit 
vor; man hat es dann mit einem sogenannten skleronomen System zu tun. Fir die Geschwin- 
digkeit jedes Fahrzeugpunktes gilt dann 
fF wer. 
ah Oq. q - (2) 
vy=1 y 
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Sie ist also eine lineare Funktion der q,- Diese sind aber bei reinem Rollen nicht voneinander 
unabhangig, denn dann haben die Berithrungspunkte der Rider mit der Fahrebene die Geschwin- 
digkeit Null. Fiir diese Punkte ist in (2) dr/dt = 0 zu setzen. Dies ergibt fiir jedes Rad zwei, 
bei 4 Radern also 8 skalare Bedingungsgleichungen. Da aber das eine Hinterrad nicht seitlich 
gleitet, wenn das andere dies nicht tut, gibt es nur 7 voneinander unabhangige Bedingungs- 
gleichungen. Nun folgt aus (2), daB fiir jeden Fahrzeugpunkt auch jede beliebige Komponente 


des Geschwindigkeitsvektors dr/dt linear in den q, sein muB, so daB die Bedingungsgleichungen 
die Form 


9 
bus qi== 9 (Tra ee a (3) 


haben miissen, worin die Koeffizienten b,,, gewisse Funktionen der Kordinaten q, sind. Demnach 
handelt es sich um ein System mit 9 Freiheitsgraden im Endlichen und 7 Bedingungsgleichungen, 
so daf} das Fahrzeug im Unendlichkleinen 2 Freiheitsgrade hat. 

c) Uber die Ableitung der Bewegungsgleichungen. Das Aufstellen der Bewegungs- 
gleichungen fiir derartige Systeme mit linearen Bedingungsgleichungen bereitet keine grund- 
sdtzlichen Schwierigkeiten, wenn man sich der sog. Lagrange-Eulerschen Gleichungen bedient?. 
Wie schon oben gesagt, fiihrt aber die Anwendung des Prinzips von D’ Alembert und des Energic- 
satzes erheblich schneller zum Ziel, ein Weg, der auch bei manchen anderen nichtholonomen 
Systemen mit Vorteil benutzt werden kann. 

Zunachst wird als sogenannte nichtholonome GeschwindigkeitsgréBe die Geschwindigkeit v 
der linken Lenkachse eingefiihrt. Auch diese mu8 sich linear durch die q, ausdriicken lassen: 


9 
2 bys Co (oe S)k (4) 
Ferner wird zu den 8 Gleichungen (3) und (4) als neunte die Identitat 
1% (5) 


hinzugefiigt. Dann hat man ein lineares Gleichungssystem mit 9 voneinander unabhangigen 
Gleichungen, die man sich nach den q, aufgelést denken kann. Die q, miissen demnach lineare 
Funktionen von v und ¥ sein. Infolgedessen muB sich die Geschwindigkeit jedes Fahrzeug- 
punktes linear durch v und y ausdriicken lassen. 

Entsprechendes gilt dann aber auch fiir die virtuellen Verschiebungen, denn die Gleichungen 
(3), (4), (5) bleiben giiltig, wenn die Geschwindigkeiten q, durch die entsprechenden virtuellen 
Anderungen 6q, ersetzt werden und an Stelle von v die virtuelle Verschiebung ds von B ein- 
gefiihrt wird. Demnach gilt: 

Die virtuelle Verschiebung jedes Fahrzeugpunktes laBt sich linear durch die 
Verschiebung 6s der linken Lenkachse und durch die virtuelle Drehung dy des 
linken Vorderrades ausdriicken. 

Daher kann die virtuelle Arbeit 6 A, der eingepragten Krafte stets auf die Form 


5A,, = Kyds + My dy (6) 


gebracht werden, worin der Index 0 andeuten soll, das vorlaufig die Relativbewegung von Motor- 
und Kraftibertragungsteilen nicht beriicksichtigt wird. 

Um die zu den so eingefiihrten allgemeinen Kraftkomponenten K, und M, gehérenden 
Bewegungsgleichungen aufzustellen, fiithrt der folgende Gedankengang ans Ziel: 

Wird bei ruhendem Fahrgestell das Lenkrad gedreht, so werden dadurch nur die Vorderrader 
und Teile des Steuermechanismus bewegt. Demnach kann sich auch eine virtuelle Drehung 67 
des linken Vorderrades nur auf die eben genannten Fahrzeugteile auswirken, fiir die im folgenden 
die zusammenfassende Bezeichnung ,,Vorderradsystem‘: benutzt wird; die virtuellen Verschie- 
bungen aller anderen Fahizeugpunkte, die nicht zum Vorderradsystem gehéren, kénnen in- 
folgedessen nicht von dy sondern nur von ds abhingen. Die an diesen Punkten angreifenden 
eingeprigten Krafte liefern daher nur einen Beitrag zu Ko. My setzt sich somit nur aus 
Momenten zusammen, die von den am Vorderradsystem wirkenden eingepragten Kraften 
herrihren. 

Demnach muB die zur Momentengréke M, gehorende Bewegungsgleichung des Kraftwagens 
gleichzeitig eine Bewegungsgleichung des Vorderradsystems allein sein. Nun wird sich in Ab- 


1 Siehe G. Hamel, Theoretische Mechanik, S. 480, Berlin, Gottingen, Heidelberg 1949. 
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schnitt e) zeigen, daB diese Gleichung verhaltnismafig leicht aufgestellt werden kann, wenn das 
Prinzip von D‘Alembert benutzt wird. Ist diese Gleichung aber einmal bekannt, dann ist das 
Aufstellen der noch fehlenden fiir Ky kein Problem mehr, denn diese 148t sich dann einfach aus 
dem Energiesatz fiir skleronome Systeme ableiten, was im Abschnitt g) geschieht. 

Bevor an die Aufstellung der Bewegungsgleichungen herangegangen wird, miissen die bei 
reinem Rollen geltenden Beziehungen zwischen den Geschwindigkeiten in expliziter Form auf- 
gestellt werden. 

d) Die Bedingungsgleichungen. Besonders leicht lassen sich die Bedingungsgleichungen 
angeben, wenn gleich geeignete nichtholonome GeschwindigkeitsgréBen benutzt werden. Als solche 
werden auBer ¥ noch eingefiihrt: die Geschwindigkeit v,; von B,, die Geschwindigkeit v, von B,. 

Diese Geschwindigkeitsvektoren miissen parallel zu den entsprechenden Radebenen gerichtet 
sein, wenn die Rader nicht seitlich gleiten. Nun besteht die Relativbewegung von B gegeniiber 
B, nur aus einer Drehung um B, mit der Winkelgeschwindigkeit b. Der Vektor der Relativ- 


geschwindigkeit mu demnach senkrecht auf v, stehen und die Grobe 1d, haben. Daraus folgt 

v4==veos7, 19, = vsiny. (7) 
Dies sind die Bedingungen fiir das Nichtstattfinden seitlichen Gleitens der Hinterrader und des 
linken Vorderrades. Aber auch fiir das rechte Vorderrad 1aBt sich die Bedingung leicht angeben. 
Die Bewegung des Fahrgestelles kann namlich in jedem Augenblick als reine Drehung um ein 
Momentanzentrum aufgefaBt werden. Dieses ist aber leicht zu finden; denn wenn die Rader 
nicht seitlich gleiten, muB es der Schnittpunkt der Verlangerungen der Radachsen sein. Daraus 
ergibt sich aber einfach (Abb. 1) 


ctg ¥,= ctg y +A mite A Z 3 (8) 


Diese Bedingung, die von dem Steuermechanismus erfiillt werden mu, ist also holonom. 
AuBerdem wird von jedem Rad vorausgesetzt, da®B es auch in seiner Ebene nicht gleitet. 

Wird der positive Drehsinn der Rader dadurch festgelegt, da beim Vorwartsfahren des Fahr- 

zeuges die Winkelgeschwindigkeiten ~, ~,, Pz, ~3 positiv sein sollen, so lauten die Bedingungen 

offensichtlich 

fiir das linke Vorderrad rg = v — c®, fir das linke Hinterrad 19, = 5 — ad, F i (9) 

fiir das rechte Vorderrad rgy,= v,-+c¢d,, fir das rechte Hinterrad r,~, = vy + bd, : f 


Dann findet man nach einfacher Rechnung, daBb die benutzten Geschwindigkeitsgréfen in fol- 
gender Weise mit v, y und y zusammenhangen: 


v,;=Vvcosy, (10a) 
i= = sin 7, (10b) 
: v yeas 

ps Sank (:os 4— 7 sin i) ; (10c) 
; v bs 
Gaps (cos oon 7 sin x) } (10d) 
vo = g(y)v mit g(y) = +71 +Asin2y +272 sin? y, (10e) 
— = ae = 1 

Wom f(y) xX mit f(%) = 1 +7 sin 24 +72 sine x > (10f) 
B= = Sime nays (10g) 
d= bein x + £04 (10h) 
Sin RaD Cine . 

=F (1—f sing) eg, ay 
P= — (sta | sin x) +e f(x) x. (10k) 


e) Aufstellung einer Bewegungsgleichung mit Hilfe des Prinzips von D‘ Alem- 
bert. Im Abschnitt c) hatte sich gezeigt, daB sich eine Bewegungsgleichung des Kraftwagens 
aus dem Vorderradsystem allein ergeben mu8. An diesem halten sich nach dem D‘Alembertschen 
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Prinzip die um die negativen Massenbeschleunigungen vermehrten eingepragten Krafte das 
 Gleichgewicht. Um also dieses Prinzip anwenden zu kénnen, muB zuerst die entsprechende 
allgemeine Gleichgewichtsbedingung fiir die am Vorderradsystem angreifenden eingeprigten 
Krafte bestimmt werden. Das ist gleichbedeutend mit der Aufgabe, die durch (6) zunichst rein 
formal eingefiihrte MomentengréBe M, in allgemeiner expliziter Form durch die am Vorderrad- 
system wirkenden eingepragten Krafte auszudriicken; denn bei Gleichgewicht muB ja nach dem 
Prinzip der virtuellen Arbeiten 6A,,=0 sein, also auch M)=0. Zur Bestimmung dieser Gleich- 


- gewichtsbedingung kann aber in (6) ohne weiteres ds—0 gewahlt werden. Dann bleibt also 
0A, ,—= M, dy, und dieses 14B8t sich leicht an Hand von Abb. 2 bestimmen. Dort bedeutet My, 
das eingepragte Moment am Lenkrad, M, bzw. M,, das Moment der am linken bzw. rechten 
Rad angreifenden eingepragten Kriafte um die Lenkachse durch B bzw. B, (solche Momente 


Abb. 2. Schema des Vorderradsystems. Abb. 3. 
kénnen z. B. bei Seitenwind auftreten), M: bzw. M; das Moment der eingepragten Krafte um 
die entsprechende Radachse (f-Achse bzw. &,-Achse), z.B. das Bremsmoment und das Moment 
des Rollwiderstandes, y= p(y) den Ausschlagwinkel des Lenkrades; dieser ist eine Funktion 
von x, die sich aus den Einzelheiten des Steuermechanismus ergibt. Hiernach ist zunachst 
6A,, = Mz oy + M, oy + M., 6%. + Mz dy + Mz dogs. (11) 
Bei reinem Rollen und 6s = 0 ist aber 
dp=y' (xox, Sye—=—flyox, rop——edy, réop,.—ef(y) dy, (12) 
so daB aus (11) wird 


64,,=|Miy'(d) + M.— > Mz +S () (Me +> Ms) bx (13) 


worin der Klammerausdruck mit M, identisch sein muf. 
Demnach lautet die gesuchte allgemeine Gleichgewichtsbedingung fiir das Vorder- 
radsystem 


M, = Miy’ (x) + M. : M; + f(x) (M., +5 M,,) =, (13) 


Um daraus die Bewegungsgleichung zu erhalten, hat man nur noch die negativen Massenbe- 
schleunigungen zu den eingepragten Kraften hinzuzufiigen. Nun kénnen die Massenbeschleu- 
nigungen der Steuerung sicherlich vernachlassigt werden, so daB zum ersten Glied in (14) nichts 
hinzukommt. An zweiter Stelle steht dort das am linken Vorderrad um dessen Lenkachse 
wirkende eingepragte Moment M,. Von diesem ist also das Moment der Massenbeschleunigung 
um die gleiche Achse abzuziehen. In Abb. 3 ist der Vektor von B, dem Schnittpunkt von Rad- 
und Lenkachse, nach einem Massenelement dm des linken Vorderrades mit Tg bezeichnet. Ferner 


sei 6 der Vektor der Absolutbeschleunigung des Massenelementes und (rg X b), die z-Komponente 


des Vektorprodukts Tz x 6. mn 

Dann ist vom Moment M, der Wert I (rz xX 6), dm abzuziehen. Mit den iibrigen in (14) 
auftretenden Gliedern ist in entsprechender Weise zu verfahren, um zur Rohform der Bewegungs- 
gleichung des Vorderradsystems zu gelangen. Diese lautet in nun ohne weiteres verstandlicher 


Schreibweise 
My = | (Fax 6), dm — © { xb); dm +f (x) ice xb), dm +© | Fs, x Be an al 
24% 
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Jetzt ist es zweckmabig, Tz X b umzuformen. In Abb. 3 soll der Punkt 0 ein ruhender Bezugs- 
punkt sein und S der Schwerpunkt des Rades, wobei der Einfachheit wegen angenommen ~ 
werden soll, da® dieser ttber dem Berithrungspunkt des Rades mit der Fahrebene liegt. Dann 
gilt zunachst 


[Fox 0) dm —$ f (nx) dm +x JF dm, tera onl 


Dabei ist aber [ (7x7) dm der Drall des Rades beziiglich B. Wird dieser mit Dz bezeichnet, 
die Masse des Rades mit mg und die Geschwindigkeit des Schwerpunktes S mit ¥,, lautet (16) 


a dD ape 
J Fax) dm = + mp Bx %. (17) 
Da jedoch ¥ parallel zu v, gerichtet ist, gilt einfach 
= dDp 
[ (ex 2) dm = "2 (18) 


Nun ist offensichtlich 
Dp = mg7,X0 +AD4+COG, (19) 
mit A bzw. C als Haupttragheitsmoment des Rades um die z- bzw. &-Achse. 


Wird der Einheitsvektor in €-Richtung mit i bezeichnet und der Einheitsvektor in z-Richtung 
mit k, erhalt man aus (19) 


{Pe _ (_mpev + ABE + CHETCHS. (20) 


Hiervon werden aber in (15) nur die z- und die &-Komponente gebraucht. Zu diesen liefert aber 
das letzte Glied in (20) keinen Beitrag, weil di/dt senkrecht auf k und i steht. 
Durch Einsetzen von (18) und (20) in (15) erhalt man dann mit den Abkiirzungen 


’ C ; C2 
schlieBlich die Bewegungsgleichung des Vorderradsystems: 
M, = A(8-+ f(x) Bs) +e m' (f(z) %, —%) . (22) 


Hierin sind nun noch #, 8, v. mit Hilfe der Bedingungsgleichungen (10) durch v und y und 
deren Abteilungen auszudriicken, um zur eigentlichen Bewegungsgleichung zu gelangen. Vor- 
laufig bedient man sich aber zweckmaBigerweise der kurzen Form (22). 
f) Die kinetische Energie Ty. Fiir die kinetische Energie Tp des linken Vorderrades 
ergibt sich leicht aus Abb. 1 
2 Tp = mgv? + AP —2megc Hv + Cg, 


woraus unter Beriicksichtigung der ersten Bedingungsgleichung von (9) und den Abkiirzungen 
(21) wird 


2 Tr = m'v? + A'S? —2em'vd. (23) 
Entsprechend bekommt man fiir das rechte Vorderrad 
2 Tr, = m'v2 + A'G?2 + 2em'v, By. (24) 


Ferner gilt offensichtlich fiir das tibrige Fahrzeug (dieses wird im Folgenden kurz mit F be- 
zeichnet) 

2T;p = mp vj + JOY + 2myz sv, 0, + Cy (g? + 9), (25) 
worin my die Masse von F' bedeutet, J das Tragheitsmoment von F um die Vertikale durch B,, 
s der Abstand des Schwerpunktes Sy von F von BB, (Abb. 1), C, das Haupttrigheitsmoment 


eines Hinterrades um seine Achse. 
SchlieBlich erhalt man mit Hilfe der Bedingungsgleichungen (10) aus (23), (24) und (25) 
nach kurzer Rechnung fiir die gesamte kinetische Energie T, des Fahrzeugs 
2 Ty = h(y) v® + m'(v? + 02) + A’ (9? + 2) + 2em’ (vx, 8, — vd) (26) 
mit 


h(y) = mp, + my, sin 2y -++ mp, sin? y , (27) 
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G, 
Le le ier | 


il Cc 
tn ( mp s + & (b— “) ; | (28) 
wae | 
a a ans my 12 = (2.22 — @ )) io 
g) Anwendung des Energiesatzes. Die noch fehlende Bewegungsgleichung wird mit 
Hilfe des Energiesatzes fiir skleronome Systeme aufgestellt. Fiir die Leistung L, der eingepragten 
Krafte gilt 


L.=4,=T. 
Aus (6) folgt aber einfach 


A, = Kyo My x. 


denn bei skleronomen Systemen gehéren die wirklichen Verschiebungen auch zu den méglichen 
virtuellen. Es muB also 


ae [T>.— Myx 27) 


sein. Setzt man hier (26) und (22) ein, ergibt sich nach einiger Rechnung fiir die gesuchte Be- 
wegungsgleichung zunachst 


Ky= h(y) 0+ ™ Ovi 4 m'(o-+e(in) + ¥ sing b-+5,) + | 
+ <7 sin (in — 8) + em’ (g(x) 6 — 8). | 


Hierin sowie in (22) sind nun noch die Bedingungsgleichungen (10) einzusetzen. Die Rechnung 
ergibt, unter Einfithrung der Abkiirzungen 


(30) 


m == my +2mp (Gesamtmasse) , } 
C G | 
ee 7 
M,—= Mp, 275, | 
ko, (y) = m, + (mp, + Am’) sin 2 y+ (m, + 22 m’) sin? y + 
422% sin y (YIP Asin 2y +72 sin? y — 1), 
/ i t . 
Foo (%) = (mp, + Am’) cos 27 roy (me + 72m’) sin 2 4 + 
cm’fA . 2cos2%+Asin2 7 I ary SERS cas) (32) 
—— | ; cos +Asin 2 ¥ +/? sin? ‘ 
al E Speci isin ya aint “aN 
hos (7) = 4 sin 1 : | om’ , 1). 
03 4) = l x 1+ Asin 2 y + /? sin? y V1 +A sin 2 7 +22 sin? x 
AG , : 5 2cos2y%+Asin2y 
ky, (x) = AT sin y-+em yl +) sin 2 ¥ + sin*7) (ib dsin2y fisin? gy?’ 
Moi (%) = kos (Z) » 
pe ve 1 Acm’ 2cos 27 +Asin24 { 
Mp2 (7%) —— Teo z(1 | eae ze 2 V1 +Asin2y +2? sin? x | 
(33) 


2 1 
Moa(%) = A (1 | I +Asin2y eieapilt 
, 2cos2y+Asin27 
Mog (Z) = —AA (1+ Asin 27 +/2sin? x)? 
fiir die Bewegungsgleichungen des hier behandelten Kraftwagens bei nicht laufendem, aus- 
gekuppeltem Motor: 
. : SOM By 
Ky = koi(y) 0 + Koo (x) ¥ x + kes (4%) % + kioy (x) x% =k, | 
1 ¢ mk 
M, = Mz y’ (x) +M. —— M; 4 (M,, + — Mz.) = | (34) 


& |] +Asin2y7+/* sin? y 


= hg (x) v + M92 (%)¥ % Mog (%) % Moa (X) 1° = Mo - 
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h) Beriicksichtigung der rotierenden Motor- und Kraftiibertragungsteile und 
Ubergang zum Vierrad-Antrieb. Bei den rotierenden Motor- und Kraftibertragungsteilen 
handelt es sich vorwiegend um achsensymmetrische Kérper, deren Achsen im Fahrzeug parallel 
zur Fahrebene gelagert sind (z. B. Schwungrad, Getrieberdder). Dadurch lassen sich diese 


besonders leicht behandeln; denn sowohl die eigene Achse als auch die Vertikale zur Fahrebene | 


durch den Kérperschwerpunkt sind stets Haupttragheitsachsen, was zu einfachen Beziehungen 
fiihrt. Nicht so einfach liegen dagegen die Verhaltnisse bei anderen Teilen, beispielsweise bei 


| 
| 
| 
! 


der Kurbelwelle. Aber dieser kurze Hinweis mége hier geniigen, weil ein Kingehen auf derartige © 


Feinheiten iiber den Rahmen dieser Arbeit hinausfiihren wiirde. Mit anderen Worten: 

Genau genommen beziehen sich die folgenden Untersuchungen auf starre Kérper mit rotations- 
symmetrischen Tragheitsellipsoiden beziiglich der Kérperschwerpunkte, deren Achsen mit den 
parallel zur Fahrebene gelagerten Kérperachsen zusammenfallen. 

Auferdem wird jetzt zum allgemeineren Fall des Vierrad-Antriebes ithergegangen, und zwar 


so, daB sich aus ihm dann ohne weiteres sowohl der reine Hinterradantrieb als auch der reine | 


Frontantrieb als Sonderfille ergeben. Im Gegensatz zum Hinterradantrieb, auf den sich die 


vorhergehenden Untersuchungen bezogen, ist beim Vierrad-Antrieb auch der Motor an das | 
Vorderradsystem angeschlossen. Aber offensichtlich andert dies grundsatzlich nichts an den | 


Uberlegungen des Abschnittes c), die zu den Bewegungsgleichungen (34) fithrten. Demnach 


kénnen der Vierrad-Antrieb und die rotierenden Motor- und Kraftiibertragungsteile lediglich | 


zu gewissen Zusatzgliedern in den Bewegungsgleichungen fithren. 
Bei voll eingekuppeltem Motor wird statt (6) gesetzt 


0A, = K ds + M by = (Ky + K,) 6s + (My + M,) 6x, | 
o4= 0A, 6AZ | 

und entsprechend fiir die kinetische Energie 
= pate (36) 


worin K, und M, aus den am Antrieb wirkenden eingepragten Kraften zu bestimmen sind und 
T, der Zuwachs an kinetischer Energie ist, der sich bei Rotation des Antriebes ergibt (unter 


,,Antrieb** seien alle zu beriicksichtigenden rotierenden Motor- und Kraftibertragungsteile 


gemeint). 
Entsprechend wird an Stelle der Bewegungsgleichungen (34) zunachst rein formal gesetzt 


K=K,+K,=h+kh, M=M,+M=m+m, (37) 


wobei aber keineswegs wegen (34) K, = k,, M,= m, gefolgert werden darf. Denn auch wenn 
keine eingepragten Krafte auf den Antrieb einwirken, also auch bei K, — M,—0, werden 
ia. k, und m,~0 sein. Die Gleichungen (34) und (37) gelten nicht gleichzeitig, sondern in (37) 
bedeuten k, und m, nur eine Zusammenfassung derjenigen Glieder, die bei voll eingekuppeltem 
Motor zu den bereits bekannten, zu ky und my zusammengefaften, hinzutreten. Bei der Be- 
stimmung dieser Zusatzglieder k, und m, wird genau so vorgegangen, wie bei der Ermittlung 
von ky und mp. 


Zur Kinematik des Antriebes ist zu sagen, da die Bedingung des reinen Rollens die Ein- | 


schaltung von 3 Ausgleichgetrieben zwischen Motor und den Radern erforderlich macht. Namlich 


je eins, das fiir den Ausgleich der Winkelgeschwindigkeiten der beiden Vorder- bzw. der Hinter- | 
rader sorgt,.und ein mittleres fiir den Ausgleich zwischen den beiden Antriebswellen, die von 


diesem zu den beiden anderen Ausgleichgetrieben fiihren. 


Nun mégen von der Kurbelwelle bis zum mittleren Ausgleichgetriebe m rotierende Teile zu _ 


beriicksichtigen sein. Das Tragheitsmoment der v.Masse (v=1, 2...m) um ihre Lagerachse sei 
J, und y, ihr Drehwinkel, wobei der positive Drehsinn so festgelegt sei, daB beim Vorwarts- 


fahren y, >0 ist. Ferner seien n bzw. n, Teile zwischen dem mittleren und dem vorderen bzw. | 
hinteren Ausgleichgetriebe zu beriicksichtigen. Die Drehwinkel dieser Teile seien mit pr, baw. | 
wa, und ihre Traégheitsmomente mit Jy, bzw. Jy, bezeichnet. Nun ist bekanntlich in camens | 
Ausgleichgetriebe die Winkelgeschwindigkeit des Treibrades proportional der Summe aus den | 


Winkelgeschwindigkeiten der beiden angetriebenen Wellen. Infolgedessen gilt 
Yr = bry (P+H2) (res aes eet) ae 


Hy = Many (P1 +s) (v=1,2,... 7), (38) | 


Wy = L(Y 2) My, (pi+ $s) (v= 1, 2,...m), 


(35) | 
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mit konstanten positiven Koeffizienten ft. Beim Auto ist gewéhnlich T=Ty5 fy =fay3 €8 iat 
aber zweckmafig, davon noch keinen Gebrauch zu machen, um spater aus den Gleichungen des 
Vierrad-Antriebes ohne weiteres Hinter- und Vorderrad-Antrieb als Sonderfiille zu erhalten. 


Fir die Momentengleichung des Vorderradsystems ergibt sich jetzt folgendes: Zunichst folgt 
aus (35) fiir ds =0 


5 Ang 2 My by . (39) 
Andererseits ist 
0:4. = py. M, oy, =f PS My, Opry, s) (40) 


worin die M, und die My, die an den betreffenden Teilen des Antriebes wirkenden eingepragten 
Momente um deren Achsen sind. Von diesen Momenten ist nur eins ein antreibendes, namlich das 
vom Gasdruck in den Zylindern herriihrende, wahrend alle anderen sich aus den Lagerreibungen 
und sonstigen Widerstanden ergeben. Die kurze Rechnung ergibt 


M, = < (f(z) —1) Mn mit Mi Moe eS a ae (41) 


M,, ist, wie man leicht nachweisen kann, im Gleichgewichtsfall gleich dem an jedem Vorderrad 
wirkenden Antriebsmoment. AuBerdem ersieht man aus (41) noch etwas wichtiges. Beim Kraft- 
wagen ist stets 


nl<z. 5 ea also = f(y) S$ 1, wenn Wee 03 
Im besonderen ist bei y <1: f(y) —1~ —2Ay, nach (41) also 
Meo jhe M,, 7: (42) 


Nun ist bei antreibendem Motor bei der Vorwartsfahrt M,, > 0, bei der Riickwartsfahrt M,, < 0, 
was zur folgenden, nicht unwichtigen Erkenntnis fihrt: 

Bedingt durch die kinematischen Verhaltnisse des Vorderradsystems ist das Antriebsmoment 
beim Vierrad-Antrieb und beim Frontantrieb bestrebt, den Betrag des Winkels y bei der Vor- 
wartsfahrt zu verkleinern, bei der Riickwartsfahrt dagegen zu vergréBern. 

DaB dieser Unterschied gegeniiber dem Hinterradantrieb, bei dem ja M,—0 ist, durchaus 
von Bedeutung sein kann, wird sich in Ziff. 3., c) bei den Stabilitatsuntersuchungen der Gerade- 
ausfahrt zeigen. 

Die Massenbeschleunigungsglieder m, und k, werden nun genau so ermittelt wie my und kp. 
Die hier iitbergangene Rechnung liefert mit den Abkiirzungen 


Dp to incl >, Diy iyo Jy = > So bacte DSi p. bg I ty ive (43) 
folgendes: 


k,, (y) = [44 2Asin2y—(4 == /*) sin? y] re 


a Yi+a sin 2 y+ /? sin? y+ 2V1+A sin 2y4+2 sin? 7] + 


+ 72M (2 cosy Asin x) (1+ J1-+ Asin 2 +7? sin® 7), 
TT, 


Ju i 
has it) = a9 VeAcos 2 — (4-40) sin 2.4] + 
ahs niet Te Ah + | 
+ 5 ah (2 cos 2 I, +A sin x) ( if yl Asin 9 x +22 sin? x ! (44) 


Ju | 7 Beare Oe (2 cos y+ Asin x) + 
rmy| 2 V1+Asin2 y+/2sin? x 


+ (Acos x — 2 sin x) (1+ [T+ Asin 2 747 sin? 7) 


c= Sa a Ju ; sin 2 y +Asin® y 
Vu ea\ | yl | Asin2y+/7? sin? 7) +78 Bos ¢-+ Asin 2)| rasa sia 


ie 


(y) = alr Ly t Asin 27+? sin? 7) + (2 cos y+ Asin x) 


| 2cos24+Asin2 x 
r 


a +Asin24-+/%sin? x)? d 
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1a.) = Wa \Ve 


aie eo olhe 2cos2y¥+Asin 2 7 _ Ju Lass sin 2 ¥ +Asin* ¥ 
malin ers aan FAsin2y+isin'y |" (A008 sin 7) 1+Asin 2y + 2? sin® x’ 
e272 (sin2y+A/sin® x)? (45) 


mM, (x) es Jy Te (1--A sin 2 X, +22 sin? x)? e 


aay e242 (2 cos 24+ Asin 2 y) (sin 2 % + Asin? x) 
ma (%) = Jv re (1 + Asin 2 7 + A® sin? x) ‘ 


Wir fassen zusammen: Bei voll eingekuppeltem Motor lauten die Bewegungsgleichungen 
der in dieser Arbeit behandelten Fahrzeuge 


K = (kg, + fa,) © + (Fog + Bag) UX + (Hog + May) % + (ky, ha.) 07> (46) 
M = My y'Q) + M.—2 M+ SQ) (Ma +> Ma) + (FQ) 1) Mu | yg 
= (1g, $1mg,) B+ (194+ Mas) © F (Ios + Mag) ZF CMay + Me) 2 - 
Ferner gilt fiir die kinetische Energie T: 


2 T = (ko, + hay) 0? + 2 (Ko, + has) 0% + (M05 + Mag) 2" » (48) 


und es ist nach dem Energiesatz 


T=Kv+Mj. (49) 


_In diesen Gleichungen beziehen sich die Koeffizienten mit dem Index a auf die Teile des An- 
triebes, d.h. auf die rotierenden Motor- und Kraftiibertragungsteile bis zu den Ausgleichge- 
trieben, diese mit einbegriffen. Fiir diese Koeffizienten gelten (44) und (45). Fiir die Koeffizienten 
mit dem Index 0, die sich auf das iibrige Fahrzeug beziehen, gelten (32) und (33). Beim Front- 
antrieb sind in den Koeffizienten (44) und (45) Jy und Jy durch 0 zu ersetzen, beim Hinterrad- 
Antrieb dagegen Jy und Jy. Demnach verschwinden beim Hinterradantrieb in den Bewegungs- 
gleichungen (46) und (47) alle Koeffizienten mit dem Index a bis auf k,, und k,,; auBerdem ent- 
fallt dann selbstverstandlich das Antriebsmoment M,,. 

Bei ausgekuppeltem oder nicht voll eingekuppeltem Motor zerfallt die Momentengleichung (47) 
in zwei Gleichungen, weil dann noch ein weiterer Freiheitsgrad hinzutritt. Eine dieser Glei- 
chungen gilt fiir das Vorderradsystem und die damit verbundenen Teile des Antriebes bis zur 
Kupplung, die andere fiir die tbrigen Antriebsteile. Aber nur in dieser kann der dritte Freiheits- 
grad wegen der iiber den Antrieb getroffenen Voraussetzungen in Erscheingung treten, so daB 
sich an der Form der Kraftgleichung (46) und der Momentengleichung (47) des Vorderrad- 
systems nichts andert. Dort sind dann nur in den Koeffizienten mit dem Index a alle Gré8en, 
die von den vor der Kupplung liegenden Motorteilen stammen, durch Null zu ersetzen. In den 
meisten Fallen kénnen dann auch noch die in diesen Koeffizienten verbleibenden Glieder, die 
vom Schaltgetriebe, von den Ausgleichgetrieben und den Antriebswellen herrithren, vernach- 
lassigt werden, so daB dann auf den rechten Seiten der Bewegungsgleichungen nur noch die 
. Koeffizienten mit dem Index 0 auftreten. Somit gelten die bereits im Abschnitt g) aufgestellten 
rechten Seiten der Bewegungsgleichungen (34) auch bei laufendem, ausgekuppeltem Motor. 

i) Uber die Kraftkomponente K, soweit sie aus Antrieb, Rollwiderstand 
und Luftwiderstand resultiert. Wahrend die Momentengré8e M in (47) bereits explizite 
in allgemeiner Form angegeben ist, sei jetzt noch die Kraftkomponente K fiir die wichtigsten 
beim Kraftwagen zu beriicksichtigenden eingepragten Krafte ermittelt. 

1. Der Antrieb. Nach den Ausfiihrungen des Abschnittes h) itber den Antrieb gilt fir 
die Leistung L, der am Antrieb wirkenden eingepragten Krafte (Antriebsmoment, Gleitreibungen 
und sonstige Widerstande) zunichst: 


1B = Dy M, Py =e Ps My, Wry Te > My, Puy . (50) 
Wegen (38) wird daraus 

I, es Was (p Se P2) ae Min (V1 El 3) = K, Vv = M, Y (51) 
und durch Einsetzen der Bedingungsgleichungen (10) 


Mm : 1S ar M m : 
K,= : (1 Ly1- Asin 2447? sin® x) + "(2 cos y+ Asin x). (52) 
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Hieraus wird naherungsweise fiir y <1 
Mm M 
Sole meee eee (53) 


2. Der Rollwiderstand. Wird die Grife des Momentes des Rollwiderstandes um die 
Radachse mit M, bezeichnet, ist die Leistung des Rollwiderstandes aller 4 Rider 


Ee = —(M,|9] + Mrs |G1| + M,, | Gal + M,, | Gal) - (54) 
Durch Einsetzen der Bedingungsgleichungen (10) findet man fiir die Kraftkomponente K, 
K, = — sgn ie ( at i ae a | 8(z) ery 
(55) 
Mr, | a | 
== = (08 % —~> sin 4 capial COs +psin 
1 


wobei beim Kraftwagen die Glieder mit a/l und c/l weet werden kénnen. Ist aufer- 
dem y <1, lautet (55) naherungsweise 


K, = — sen of (M,+ M,)-+ = (Ma + M,, eee athe (56) 
1 1 
3. Der Luftwiderstand. In Abb. 4 sei P, diejenige Komponente des resultierenden Luft- 
widerstandes, die parallel zur Fahrebene gerichtet ist. Dann ist die Leistung L, des Luftwider- 
standes 
DL, = — P, cosa vp 
Nun ist vy —v,-+e#,, so daB sich wegen (10a, b) fiir die Kraft- 
komponente K, des Luftwiderstandes ergibt 


KG == — PB cosa (cos y+ i sin x) ‘ (57) 

und bei y¥< << 1 
Ki ~ — Peosa(1+-7). (58) 
3. Beschrankung auf kleine Winkel y. a) Der allgemeine Fall. 
Die Bewegungsgleichungen werden bedeutend einfacher, wenn y < 1 
ist, was beim Kraftwagen in der normalen Fahrt meistens der Fall 


ist. Entwickelt man die Koeffizienten (32), (33), (44) und (45) der Be- 
wegungsgleichungen (46) und (47) nach Potenzen von y und bricht 


hinter den linearen Gliedern ab, so findet man nach kurzer Zwischen- eee 
rechnung 
ky, = m,-+ 2 (Mpg +A my, = Mrs ! Am’ (. 22 m’ + 2A— m') 1: 
‘ Vic (59) 
hy = (2 — hem’) x. hy, = —2 Alem! + (7 —2hem’|), | 
My, = ko, » My, =27 them! —24(7 “+ aem') 7, | (60) 
Me aA (2 Ay), ny, == = 2A (1 —SAq) 
hee ie Padi 
oy 7 4 in | re | Sd 14%), 
ha = 22 (8 ! Jy 5 =) (7 ! = a 2 Alc ale 
i i ip ry Ta 2 7 : ; 
Ty J ( (61) 
k,, — —4 A — + aA ae ry | 
Jy, Jm 
hy = —2eA(E +R) 544). 
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ie == hk; m= 2B (TE fe {> | 
2 (62) 1 
m,, = 0 m,, = 42 Iv x 


Ferner erhalt man fiir die MomentengréBe in (47) 


(Mn ote M:,) “i M., Hie (63) | 


b) Ein Zahlenbeispiel. Um einen Uberblick iiber die GréSenordnungen der Koeffi- 
zienten k und m zu bekommen, sei nun ein Zahlenbeispiel durchgerechnet, und zwar mit Werten, 
wie sie bei einem mittleren Personenkraftwagen auftreten. 


c 
r 


M= Mz y' (x) + M+ M., +— (M,, — M,) 2A| 


Abmessungen: 

l=26m: t= 1,3m; «1 Bem; s]0,6om;," T= 034m ic —.0 aaa 
Massen: 

m= 130 kgm—sec?;  mg—= 2,5 kgm—'sec?; mp = 125 kgm“1sec?. 
Tragheitsmomente: 


J, = 230 mkgsec? (Tragheitsmoment von my um Vertikale durch Schwerpunkt von mp); 
C = 0,13 mkgsec?; A = 0,065 mkgsec?. 


Antrieb (nach einem Zahlenbeispiel von W. Kamm u. C. Schmid '): 


J Ubersetzungsverhialtnis 
one 1. Gang 
Motortriebss ane 0,447 
Schwungrad ... . | 3,04 Dee Aso 
Kupplung. .. .— 0,090 |f 97 |j = 25,65 
Getriebe und Gelenk- | 
Welle Ameer nae 0,05 Zale osOo 
ba he aay 
Ausgleichgetriebe . . | 0,60 1 1 


Durch Einsetzen dieser Zahlen in die Beziehungen (59), (60), (61) und (62) ergibt sich fiir die 
Massenbeschleunigungsglieder der Bewegungsgleichungen (46) und (47) 


bei Vorderradantrieb in direktem Gang 


ky = (134,5-++ 67,2 7) v + (33,6 — 53,0 7) vy — 0,038 y y — (0,181 — 0,155 x) 72, 
k, = (10,2+ 5,1 7) v+ (2,55 + 0,64 y) v 7 — 0,255 v 7 — (0,255 — 0,638 y) 72, 


bei Vorderradantrieb im 1. Gang 


ky = (205 + 102 x) o + (51,2 + 12,8 y) vy —5,12 7 ¥ — (5,12 — 12,8 x) 72, 


64 
bei Hinterradantrieb im direkten Gang “g 
k, = (10,2 + 5,1 y) » + (2,55 — 9,56 x) vy, 
bei Hinterradantrieb im 1. Gang 
k, = (205 + 102 y) v + (51,2 — 192 y) vx, 
bei Vorderradantrieb im direkten Gang 
m, = — 0,0384 x v + (0,143 — 0,117 x) v % + (0,136 — 0,136 x) ¥ — (0,0678 — 0,170 x) 72, 
Mm, = — 0,255 ¥ v — 0,0638 x v x + 0,006 x 72, (65) 


bei Vorderradantrieb im 1. Gang 


m, = — 5,12 yo — 0,128 y vy + 0,013 y y*. 


ae W. Kamm und C. Schmid, Das Versuchs- und MeBwesen auf dem Gebiet des Kraftfahrzeugs, Berlin 
8. 
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Damit erhalt man bei voll eingekuppeltem Motor die folgenden Bewegungsgleichungen: 
Bei Vorderradantrieb im direkten Gang 


K = (144,7 + 72,3 x) v + (36,2 — 52,4 7) v y — 0,293 y 7 — (0,436 — 0,793 y) 72, l (66) 
M = — 0,293 x v + (0,143 — 0,181 x) v x + (0,136 — 0,136 x) % — (0,0678 — 0,176 x) x2. f 

Bei Vorderradantrieb im 1. Gang 

K = (340 + 169 x) v + (84,8 — 40,2 y) v x — 5,16 ¥ ¥ — (5,30 — 13,0 x) v7, \ (67) 
M= —5,16 xy v + (0,143 — 0,245 y) v x + (0,136 — 0,136 x) ¥ — (0,0678 — 0,183 y) ¥?. J 

Bei Hinterradantrieb im direkten Gang 

K = (144,7 + 72,3 x) v + (36,2 — 62,6 x) vy — 0 0384 y y — (0,181 — 0,155 y) 72, } 

im |. Gang | 

K = (340 + 169 x) v + (84,8 — 245 x) v x — 0,0384 y ¥ — (0,181 — 0,155 x) 7”, , (68) 
M =m, = — 0,0384 x v + (0,143 — 0,177 x) v ¥ + (0,136 — 0,136 x) ¥ na 


— (0,0678 — 0,170 y 
Diese ‘Zahlen zeigen folgendes: In den Kraftgleichungen betragt der Anteil der Antriebsteile 


an den Gliedern mit v und vy beim Vorderradantrieb etwa 1/,,, dagegen an den Gliedern mit 
und y? mehr als die Halfte. Im ersten Gang entfallt auch bei den Gliedern mit v und vy mehr 
als die Halfte auf den Antrieb, an dem allein die Schwungscheibe mit 85 % beteiligt ist; bei den 
Gliedern mit y und y? betragt der Anteil des Antriebes fast 100 %. Die Kraftgleichungen des 
Vorderradantriebes unterscheiden sich von den entsprechenden des MHinterradantriebes im 
wesentlichen in den Gliedern mit 7 und y?. Die geringste Rolle spielen diese Glieder beim Hinter- 
radantrieb, die gré{te beim Frontantrieb im 1. Gang. In den Momentengleichungen machen 
sich die Antriebsteile in den Gliedern mit vy, 7 und y? kaum bemerkbar, dagegen vergréfern 
sie den Betrag des Gliedes mit v im direkten Gang auf das 7,6fache des Leerlauf-Betrages und 
im 1. Gang sogar um das 134fache. 

c) Ist die Geradeausfahrt mit konstanter Geschwindigkeit stabil? Als erste 
Anwendung sei folgender verhaltnismaBig einfacher Fall untersucht: Das Fahrzeug fahre zu- 
nachst mit y = 0 und konstanter Geschwindigkeit Vv =v). Dann miissen die eingepragten Kraft- 
komponenten K und M gleich 0 sein. An au®eren eingepragten Kraften seien lediglich konstante 
Rollwiderstande und der Luftwiderstand zu beriicksichtigen, der bei v= v, auch konstant 
sei, so da®B dann auch die auf die Rader ausgetibten Antriebsmomente konstant sein miissen. 
Ist nun eine solche Bewegung stabil? Das soll hier heiBen: Kommt das Fahrzeug bei kleinen 
Abweichungen von vy = 0 und v = v, von selbst wieder auf y= 0 oder bedarf es dazu standig 
einer Steuerung durch den Fahrer? Selbstverstandlich hat eine solche Frage nur Sinn bei nicht 
selbsthemmenden leichtgangigen Steuerungen. 

‘Wird fiir die Geschwindigkeit v = v, + w gesetzt und wird in itblicher Weise angenommen, 
daB y und w samt ihren Ableitungen so klein sind, da8 in den Bewegungsgleichungen alle Glieder 
fortgelassen werden kénnen, die von hdherer als erster Ordnung klein sind, so ergibt sich aus 


der Momentengleichung (47) unter Beachtung von (60) und (62) 
A ae a 
M= (24 them’) vj +24 7. (69) 


Bei den vorausgesetzten kleinen Stérungen spielen also die Teile des Antriebes noch keine Rolle. 
Far M bekommt man aus (63), wenn dort M;, = M,= M,, — 0 gesetzt wird und fiir die Mo- 
mente um die Radachsen (¢- und &,-Achsen) die beretis in (54) eingefiihrten Momente der Roll- 
reibung eingesetzt werden: 


= < (M,, M,) 2/4— (M,, M,,) x - (70) 


Da aber bei y= y= ¥ = 0 auch M=—0 sein mu, miissen die als konstant angenommenen 
Momente der Rollwiderstande gleich gro8 sein, also M,,— M,. Demnach entsteht aus (69) 
und (70) 


2A E+ (2 + Jem’) vy) % +242 (M,, — M,,)4=0. (71) 
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Das an jedem Vorderrad wirkende Antriebsmoment M,, wird bei der gestérten Bewegung im 
allgemeinen nicht konstant bleiben. Da aber auch konstanter Luftwiderstand bei v = vy voraus- 
gesetzt wurde, werden die Anderungen von M,, um so kleiner sein, je kleiner die Abweichungen 
von der stationaren Bewegung sind. Infolgedessen fiihren diese Anderungen in der Bewegungs- 
gleichung (71) lediglich zu GréBen, die von zweiter Ordnung klein sind, so daf sie fortgelassen 
werden kénnen. Demnach ist in (71) fiir M,, das bei v= v,) und y—0 wirkende konstante 
Antriebsmoment einzusetzen. Dann sind aber samtliche Koeffizienten in der Bewegungs- 
gleichung konstant, so da® nun (71) von der Kraftgleichung getrennt untersucht werden kann. 
Im Laufe der Zeit geht y dann und nur dann gegen 0, wenn alle Koeffizienten gleiches Vorzeichen 
haben. Nun ist stets A’ > 0, so daB also bei Stabilitat auch 


(24 +4em’) pel nim san be sa (72) 


sein mu. Dieses trifft nur zu, wenn 
Y > 0, M,,— M,> 95 (73) 


also nur bei der Vorwartsfahrt und beim Vierrad- oder Vorderrad-Antrieb, da beim Hinterrad- 
antrieb M,, = 0 zu setzen ist, 

Auch hierzu ein Zahlenbeispiel, wobei sich wieder alle Werte auf das m-kg-sec-MaBsstem 
beziehen sollen. Mit den im Abschnitt b) benutzten Werten bekommt man fiir die ersten Glieder 
in (71) 


2 Aly 4 (2 he m’) v7 = 0.136 % +.0,143 v7. (74) 


Im letzten Glied in (71) 1a8t sich das konstante Moment M,, — M,,, das bei v= v) und y= 0 
wirksam ist, aus der Kraftgleichung bestimmen. Bei y= 0 und v= v, mu K= 0 sein. K 
setzt sich hier zusammen aus der Antriebskraft K,, dem Rollwiderstand K, und dem Luftwider- 
stand K;. Nach den fir diese Krafte aufgestellten Beziehungen (53), (56) und (58) erhalt man 


zunachst, wenn dort vy = 0 gesetzt wird 
K— K, -- K+ Ky, | 


Z | M,, + M,,) — Pyeosa=0. | 


(Me) eh (75) 


ro | 


Unter Beschrankung auf den einfachen Fall, daf alle Rollwiderstandsmomente gleich M, sind 
und der Anstrémwinkel « <1 ist, entsteht aus (75) fiir den Vorderradantrieb (M,,, = 0): 


M,, —M,,= M, ++ Py. (76) 


Wird andererseits fiir den Rollwiderstand K, die einfachste GesetzmaBigkeit zugrunde gelegt, 
das er namlich proportional dem Gesamtgewicht G des Fahrzeuges ist, also K,— —f,G, mit 
konstantem ,,Rollwiderstandsbeiwert f,, so findet man durch Vergleich 


Mocitf 


Als Beispiel sei f, = 0,03 (TeerstraBe) gewahlt. Damit wird M,— 3,24 mkg. Ferner wird fir 
den Luftwiderstand P, einfach gesetzt: 


P,= + vec, F, (77) 


mit @ als Luftdichte, F als Fahrzeugquerschnittsflache und c,, als Widerstandsbeiwert. Bei 
mittlerem Luftdruck und mittlerer Temperatur ist 0/2 — 0,0637. Setzt man fiir einen mittleren 
Personenkraftwagen beispielsweise c,— 0,5, F—= 3 m?, ergibt sich fiir den Luftwiderstand 
P, = 0,0955 vg. Mit diesen Werten lautet (76) 


M,, — M,, = 3,24 + 0,0162 v2. 


Wird dieses und (74) in die Bewegungsgleichung (71) eingesetzt, kommt man bei Vorderrad- 
antrieb auf 


xy + 1,053 v, ¥ + (3,50 + 0,0175 oe ga Ue 
bei Hinterradantrieb auf 


x + 1,053 », 7 —3,507—0. 
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Die allgemeine Lésung dieser Gleichungen lautet bekanntlich 


y= Ae? + Be" (wenn J, #A/,). 
Fihrt man die Rechnung durch, ergeben sich folgende Vergleichswerte : 


e i pe 

iid V.-Antrieb | H.-Antrieb V.-Antrieb | H.-Antrieb 
36 —0,53 +0,33 —10,9 
(2, —0,51 +0,18 —21,2 

108 —0,63 +0,12 —31,7 


Demnach verschwindet das Glied mit ¢*' stets sehr schnell, aber auch das Glied mit e”' wird 
beim Vorderradantrieb schnell klein, nach etwa 4sec betrigt es nur noch ungefahr 1/,) seines 


Anfangswertes. Dagegen wachst beim Hinterradantrieb ¢' ziemlich rasch, und zwar bei 
kleineren Geschwindigkeiten starker. 


4, Uber die kraftefreie Bewegung des vierridrigen Wagens. a) Die ersten Integrale der 
Bewegungsgleichungen. Fir den gewodhnlichen vierradrigen Wagen und die anderen in 
Ziff. 2., a) genannten Fahrzeuge mit 2 = b’/l = 0 vereinfachen sich die allgemeinen Bewegungs- 
gleichungen ganz erheblich. (34) lautet dann einfach 


Ky = (m, + mp, sin 2 ¥ + m, sin? y) v + | 


+ (mr, cos2 x +t sin2 x) 0 +25 sin xi | (78) 


M, = 24 sinyb +24 cosyoy +2 4’j. (79) 


Diese Gleichungen lassen sich exakt lésen, wenn die eingepragten Kraftkomponenten K, und Mo 
standig 0 sind, was z.B. der Fall ist, wenn das Fahrzeug auf horizontaler Ebene rollt und die 
Schwerkraft die einzige am Fahrzeug wirkende eingepragte Kraft ist. Zunachst liefert der 
Energiesatz sofort ein Integral, denn bei K, = 0 und M,=0 muB8 nach (29) die kinetische 
Energie T, konstant bleiben. Aus (48) erhalt man 

konst. = 2 T, = (m, + my, sin 2 7 + mg sin? x) v? +44 sin yoy +2A’y?. (80) 


Aber auch die Momentengleichung (79) 148t sich leicht einmal integrieren, denn statt (79) kann 
geschrieben werden 


padeili's : 
My=24' 5 (i, +7 sing) = 0. (81) 


Hierin ist das zweite Glied in der Klammer nach der Bedingungsgleichung (10b) nichts anderes 

als die Winkelgeschwindigkeit #, des Fahrgestelles. Ferner ist die Winkelgeschwindigkeit @ 

des Vorderradgestelles um die Lenkachse y¥ + #1, so daB aus (81) einfach folgt 
M,=2A4'9=0, | 

‘ , : 82 

P= 0 = konst.=% +b, = 4 +7 sing. | (82) 


b) Zwei Sonderfalle. Aus (80) und (82) lassen sich sofort zwei Sonderfalle ablesen: 
Lito 00 Kougt. 
Das Vorderradgestell bewegt sich geradlinig, so daB der Punkt B, die bekannte Huygenssche 
Traktrix beschreibt. 
Poy (ie Ue 4 = y, = konst., v = konst. 
Das Fahrzeug bewegt sich mit einer konstanten Geschwindigkeit v auf einer Kreisbahn, und es 
ist nach (82) 


SIN Yo 


Hierin kann //sin y, nichts anderes als der Halbmesser des von der Lenkachse durchfahrenen 
Kreises sein, was ein Blick auf Abb. 1 sofort bestatigt. 
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c) Ein weiterer interessanter Sonderfall. Von nun ab sei stets w 40, ¥ 40 
vorausgesetzt. Setzt man in (80) 4 = @ — 0, und v= (I/sin y) 0; nach (10b), bekommt man 
nach kurzer Rechnung 


Ate sin = hey (83) 
: l + /1—2xsin? y+ usin 27 x 


mit den Abkiirzungen 


lp PEG salege m 
i peas eee ae (84) 
my, 2m, 2 my, 
Hiermit ergibt sich 
v = (85) 


— + /1—2sin? y+ psin2y ° 

Untersucht man nun die Systemkonstanten x% und yw, so sieht man, da der Fall x = 0, w= 0 
durchaus vorkommen kann. In diesem Sonderfall mu aber nach (85) die Geschwindigkeit v 
der Lenkachse wahrend der ganzen Bewegung konstant bleiben. Das bedeutet, da ja auBerdem 
auch die Winkelgeschwindigkeit § des Vorderradgestelles konstant bleibt, daB die Lenkachse 
mit konstanter Geschwindigkeit v, einen Kreis beschreibt. In diesem Falle la8t 
sich auch die Differentialgleichung (83) fiir y ohne weiteres lésen, so das man dann imstande 
ist, die Schleppkurve irgend eines Punktes des Fahrgestelles zu bestimmen, also beispielsweise 
des Schnittpunktes B, von Lingsachse und Hinterradachse. Die im Abschnitt e) durchgefithrte 
Rechnung fithrt, wie die Abb. 5 und 7 zeigen, zu zwei wesentlich voneinander verschiedenen 
Typen von Schleppkurven, je nachdem ob | kleiner oder gréBer als der Halbmesser @ des von 


der Lenkachse beschriebenen Kreises ist (« = t) ‘ 
d) Zur Integration im allgemeinen Falle (x40, w~ 0). Aus (83) ergibt sich 
GI ra NES 2s ior ETD oe pe See eae, (86) 
+ V1 —2x sin? y+ usin 2 lw 
Durch Trennen der Variablen bekommt man 
or= d= | ee dy (87) 
=e V1—2x sin? y +psin 2 ¥—tTsiny 


und durch einfache Umformung 


oat=?v= Js +7 Jo, 


mit ‘ 
1 L-x +x cos 24% +ysin2y g 

= —_ : 2 
Jt 2 | 1—o+ocos2y+psin2y¥ (2x). { (88) 
Re | pie ea eee ( 

c 1—o+ocos2y+psin2 7 ; 

2 
26 =2x +7 =2x- cee 
Pw? 


Die allgemeine Lésung von J, kann aus jeder Formelsammlung entnommen werden. Auch J, 
ist elementar lésbar; denn die Substitution tg y = =z liefert 


Jom ( zl 2wet 2) 2] : 
2 bea So Vi+2p2+(1—2%) 2 


dz. (89) 


Dieses Integral 148t sich durch Aufspaltung in Partialbriiche auf Integrale von der Form 
| asthe Wasp SiS ble te 
J @—a)’V1+2yz2+0—2o0)2 


(vy >0, ganzzahlig) , 


zuriickfiihren, die elementar lésbar sind. Die allgemeine praktische Durchrechnung ist aber sehr 
umstandlich und diirfte ziemlich umfangreich werden. Deshalb wird sie in dieser Arbeit nicht 
durchgefiihrt. Nur der Fall = 0 1aBt sich ohne grofBe Schwierigkeiten vollstandig durch- 
rechnen, denn bei 4 — 0 fiihrt die Substitution 2? = u auf erheblich einfachere Integrale. 
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e) Die Lésungen des Falles fe = 0. 1. Der Sonderfall x = 0. Nach (87) ist 


a iil dy aoe Tt Ot 

y | 1+tecosy’ eee a) ees ) 
Nun hat sich schon im Abschnitt c) gezeigt, daf§ in diesem Sonderfall die Lenkachse mit kon- 
stanter Geschwindigkeit v) einen Kreis durchfahrt. Fiir den Halbmesser 0 dieses Kreises gilt 
dann offensichtlich @ = v,/w. Demnach ist 7 = o/l. Als Liésung von (90) ergibt sich fiir 


2 1—t 
ibs = ———— ei 
v< CH) i arctg (= te | A (91) 
t=1: b=igt, (92) 
2 Geek y 
se — == 9} — (ipgeies 
a : eal oea arta (ST eZ), | 
wenn die Spitze der Schleppkurve auferhalb des Kreises liegt, sonst (93) 
= 2 ; feb el ee | 
ave Wr Te i tg 9 ’ fA ee a0 se Ie 


Abb. 5. 


Abb. 5 u. 6. Bahnen der Lenkachse und 
Schleppkurven. 


Abb. 6. 


Hiernach sind die Schleppkurven in den Abb. 5, 6 und 7 berechnet worden. Abb. 6 zeigt die 
Schleppkurve im Grenzfalle 9 = / fiir 4 volle Umlaufe des Vorderradgestelles; mit ¢ + oo fiihrt 
die Schleppkurve zum Kreismittelpunkt. 
2. Der Sonderfall 20=2x+7?=0, dh. c=;2=—2x,%<0. Aus (88) wird 
1 
Ju= x [ L—x + xe0s2y) dy), | 
: (94) 


Jy= | sing [12x F 2x cos 7 dy. | 
Das erste Integral 1aBt sich sofort lésen: 


Jy=(l—x) y+ sin2y, (95) 


Ingenieur-Archiy 
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beim zweiten liegt die Substitution cos y= z nahe. Damit wird 


= ie 2% +2ue2dz 


mit der Lésung 


I To Sy ero ape ee ; f = 03 
Jg= —-> cosy Yl — 2x sin BY pa anesin (/ 2% cos 7). 


(96) 


Hiermit bekommt man fiir ? = J,+ 7 J,, wenn wieder zum Winkel y = y + a iibergegangen 


Abb. 7. 


Abb. 8. 


Abb. 9. 
bbe Abb, 10. 


Abb. 7-10. Bahnen der Lenkachse und Schleppkurven. 


-0 
# = 0395 
Uv =18 
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wird, 


or —— Bn ing 
| * sin y yl — 2x cos* y ae arcsin (V/3 = - sin ») 2 (979 


o= (l—x) y +—*sin2y — 


Nun kann die Bahnkurve der Lenkachse errechnet werden; denn aus 


dx dx 0, 
i ae = aa @ = UVcos | 
: , (98) 
pe = MAE 
; ein ek geld sin } | 
bekommt man nach (85) 
x ; cos } dd 
—=—T — ; 
l | /1 —2x sin? y+psin2 x | 99 
eat f sin § dd : 99) 
l | V1—2xsin? y+psin2y 
Hiermit wurde das Beispiel x = — 0,25, also 7 = V0.5 5 = 0,707 zahlenmabig durchgerechnet. 


Die Bahnkurve der Lenkachse (s. Abb. 8) wurde durch numerische Integration von (99) be- 
rechnet. 


3. Die nichtperiodischen Bewegungen mit 2¢>1, d. h. ?>1—2x. Glei- 


chung (88) ergibt 
1—x+x cos 24 
2a Stee 1—o-+ cos 27 d(2 x) 
mit der Lisung 


a eee eee Ar Tq (V2o0—1 tg x) . (100) 
Aus (89) erhalt man mit der Substitution z*? = u 
= | ad teste (101) . 
(i+ u) [1+ (—2o)u] ic + (1—2 x) 
Durch Partialbruchzerlegung entsteht 
—_ rt {OS ee = | Gipge Ua ean (102) 
o J 214+u)/1+0—2x)u Oo } 2f1+(1—2o0)u}/14+ 1 —2x)u 
Diese Integrale werden mit Hilfe der Substitution 
yl+(0—2x)u=w (103) 
gelést. Die Ergebnisse der hier iibergegangenen Rechnung lauten 
15 /2o—1 siny 
o=~ . 
A ae (, V1 —2 x cos? | 
|— op oe or Gin (V; — _ gin ¥) , wenn x> 0, | (104) 
ML 
+ SOY, (tp Shee 2 =i RES) 
ae * are sin (\/ eae ein y) wenn x <0, 
1—2x 
und 
an x /2o—1 sin Y1 
See ae oe za (, et | : 


Be “xe 
| a=: a Wr Gin es sin vs) A wenn x > Us | 

+ J ° Ui = 
| + eee are sin (2 sin vy) » wenn x< 0 | 


Mit (104) wurden die Kurven der Abb. 9 berechnet und mit (105) Abb. 10, und zwar mit 
fa == d20, alone =—i27 


; ee ee ee 


25 


356 Stiickler: Uber die Differentialgleichungen fiir die Bewegung usw. Ingenieur-Archiv 


2 
Vo 


4. Die periodischen Bewegungen mit 2¢ <1, d.h. figs? <1—2zx. Fiir diese bekommt 


man wegen 
—itiAr Tg (uz) = arctg z 


ohne weiteres aus (104) 


2 ert 
yD = a Ww | ae are tg | ieee + 


@ 20/1—26 t /1—2% cos? y + (1—2 x) cosy | 
| = )2% Ar Gin (| pace sin v) » wenn x>0, | ( (106) 
re 20 ar ap 
+7 lS arc sin (ee sin ¥) , wenn x < 0. | | 
2 1— 2x J 
5. Der Grenzfall2o=1, d.h. 7? = jot —1]1-— 2x. In diesem Fall nimmt das zweite 


Glied in (104) die unbestimmte Form ~ an. Die Anwendung der Bernoulli-L’ Hospitalschen 
Regel ergibt dann 


Tsiny ; 


G=2x ———— 
Xi | Vio x cos? y +t cosy 


—V2x Ar Sin be sin ») » wenn x>0, (107) 
J 


y 4 f= 2A ave sin 


| 
S { 

siny) wenn x <0. | 
Fiir + = 0 wird hieraus wegen 7 = I 


2 EN 
eee WEEE 


~ 
eee 


wie es nach (92) sein muf. 


(Eingegangen am 12. April 1952.) 
Anschrift des Verfassers: Dr.-Ing. Bernd Stiickler, Berlin-Steglitz, Mariendorfer Str. 32 b. 
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Zur Berechnung der Randstérungen kreiszylindrischer Tonnenschalen. 
Von W. Zerna. 


1, Einleitung. Die elastizitatstheoretische Berechnung des Spannungs- und Verformungs- 

zustandes kreiszylindrischer Tonnenschalen — einer fiir das Bauwesen bedeutsamen Schalen- 
form — ]a8t sich nach der Membrantheorie in einfacher Weise erledigen, die allerdings im all- 
gemeinen nicht die gewiinschten Randbedingungen lings der begrenzenden Erzeugenden zu er- 
fiillen gestattet. Die von den Randern ausgehenden Stiérungen des Membranspannungszustandes 
kénnen nur mit Hilfe der Biegetheorie erfaBt werden, was nach den bekannten Theorien! auf 
nicht unerhebliche Schwierigkeiten stéBt. 
_ _Inneuerer Zeit durchgefiihrte Untersuchungen iither Schalen? allgemein und Zylinderschalen * 
im besonderen erlauben nun die Entwicklung einer wesentlich einfacheren Theorie zur Berech- 
nung des Randstérungsproblems kreiszylindrischer Tonnenschalen, woriiber nachfolgend be- 
richtet werden soll. 


2. Elastizitatsgesetz und Gleichgewichtsbedingungen. Die Schalenmittelflache einer kreis- 
zylindrischen Tonnenschale sei in iiblicher Weise durch Zylinderkoordinaten x und @ beschrieben, 
der Kriimmungshalbmesser der Mittelflache sei a, W 
die Spannweite zwischen den aussteifenden Binder- = ve 
scheiben | (Abb. 1). SP 


A, 


b 2e| ve Mop 


Se 


YW te f 


Abb. 1. Mittelflache mit Koordinaten und Abb. 2. Schalenelement mit SchnittgroBen, 
Abmessungen. a) Langskrafte, b) Momente und Querkrafte. 


Die SchnittgréBen, deren Bedeutung und positiver Wirkungssinn aus Abb.2 ohne weiteres 
ersichtlich ist, werden wie folgt bezeichnet : 


Nos Ne = Langskrafte, 
Nyo> Nex = Schubkzifte, 
ORE NED = Querkrafte, 
M,, M, = Biegemomente, 


M,.,, M,x = Drillungsmomente. 
Die Ableitungen nach ~ und x/a seien durch Punkte und Striche gekennzeichnet. 
Bedeuten u, v, w die Verschiebungen eines Punktes der Schalenmittelflache in Richtung der 
Koordinaten x, g und der Normalenrichtung (positiv nach aufen), so lassen sich die Dehnungen 


Se und die Gleitung y,, der Mittelflache durch 


vy 
« 


Ww a: u’ v 
Sie ose Cor 9 Yh ee (1) 


ausdriicken. 


1 Vgl. U. Finsterwalder, Ing.-Arch. 4 (1933), S. 435 F. Dischinger, Beton u. Eisen 34 (1935), S. 257; 
: - bsen, Bauing. 20 (1939), S. 394. 
* Ava AE. Green ed W. Zerna, Quart. Journ. Mech. and Applied Math. 3 (1950), S. 9. 
3 R. S. Jenkins, Modern Building Techniques, Bulletin No. 1, London 1947. 


i) 
on 
* 
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Fiir das Elastizitatsgesetz, das den Zusammenhang zwischen den SchnittgréBen und den Ver- 
schiebungen angibt, ergibt sich als erste Naherung die Form? 
No = Deel e2) | 


N, =D(e, +pe,) » (2) 


Nos = Negi= DZ" eg, | 

M, == (w" +p"), 

M, = = (w yw’), (3) 
My = Myx = (1 — 1) 0 


Darin ist 4 die Querkontraktionsziffer und D und B bedeuten die Dehnungs- und Biegungs- 
steifigkeiten, die sich unter Verwendung des Elastizitatsmodul E, sowie der Schalendicke t gemahs 


Et Ee 
dk npetatis pete we 4 
y Lae? le 12(1—,) (4) 
berechnen. 
Die exakten Gleichgewichtsbedingungen der unbelasteten Schale lauten ? 
Nz +Nox=9, (5) 
No +Nig—Q,=9, (6) 
Q, +Q. +N, =0, (7) 
1 © / 
Q, == a (M, Sik: May) 2 | 
es (8) 
a= Ler tag. | 


Wird nun fiir die Gleichgewichtsbedingungen die gleiche Art der Naherung verlangt, wie sie 
beim Elastizitatsgesetz nach (2) und (3) vorhanden ist, so laBt sich nachweisen, daB Q, in der 
Gleichung (6) vernachlassigbar ist *, so daB (6) durch die angenaherte Gleichgewichtsbedingung 

Neh era) (9) 
ersetzt werden kann, 

Es sei betont, da die Theorie, welche die angenaherte Gleichgewichtsbedingung (9) ver- 
wendet, durchaus nicht schlechter ist als die itblichen Theorien 4, die (6) verwenden; denn die 
grundlegenden Annahmen der Schalentheorie bedingen ohne dies nur eine erste Naherung, die 
aber in den iiblichen Theorien nicht konsequent durchgefiihrt worden ist. 


3. Differentialgleichung. Die beiden Gleichgewichtshedingungen (5) und (9) werden durch den 
Ansatz 


Ne OS IN 3 Ds 3 Nag = — &” (10) 
identisch erfillt. Darin ist ® eine Spannungsfunktion. Aus (2) folgt 
N_ —wN 2 
anu? ~ 
cables IT bay 
No ING 
oe == DO —,) ’ (11) 
DS 


und aus (1) die Identitat 


Ep TEx Veo 0 : (12) 
_ + Vgl. R.S. Jenkins, a. a. O. und A. E. Green und W. Zerna, a. a. 0. Die in die Schalentheorie einge- 
fiihrten grundlegenden Annahmen lassen alle Verfeinerungen des Elastizitatsgeseizes durch Hinzunabme 
weiterer Glieder wenig sinnvoll erscheinen. 
2 Vel. W. Fliigge, Statik und Dynamik der Schalen, §. 110. Berlin 1934. 
3 Vel. R. S. Jenkins, a. a. O. 
4 Vgl. FuBnote 1 S. 357, 
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Kinfihrung von (10) in (11) und dies dann in (12) liefert 
4Ap—2E—#) ag 9, (13) 


2 2 
wenn J = a ue ++ i den Laplaceschen Operator bezeichnet. 


Mit (13) ist eine Differentialgleichung fiir die Spannungsfunktion ® und die Normalverschie- 
bung w gefunden. Um eine weitere Gleichung zu erhalten, wird (8) in (7) eingesetzt. Dies ergibt 
wenn noch (3) und die zweite Gleichung (10) beriicksichtigt werden 


3 
AAw+ > 0" 0, (14) 
Die beiden Differentialgleichungen (13) und (14) bilden ein System homogener partieller Diffe- 


rentialgleichungen vierter Ordnung mit konstanten Koeffizienten zur Bestimmung der beiden 
unbekannten Funktionen @ und w. Die beiden Differentialgleichungen werden zweckmafig durch 
eine einzige ersetzt. Um sie zu erhalten, kann wie folgt vorgegangen werden. Die Gleichung (13) 
wird mit einem konstanten Faktor k multipliziert, dann werden beide Gleichungen addiert. Dies 
ergibt 


DUA—p?*)k a4 eae 
AA(w +k ®) - (w pad ane 0) = 0 (15) 
Der Faktor k soll nun so bestimmt werden, daB 
at 
~~ BD —p*)k 
ist, woraus folgt 
k= 7K 
mit 
ee fiedea 
K=(*) = ie yt (16) 
Mit Einfiihrung der komplexen Funktion 
PY=w+tik® (17) 
und 
one t 
aloe (18) 
wird aus (15) 
AAY —{P" = 0. (19) 


Dies ist eine komplexe Differentialgleichung vierter O1dnung, die das hier betrachtete Problem 
der Kreiszylinderschale beschreibt. Ist die Funktion Y gefunden, so liefert gemaf (17) ihr Real- 
teil die Normalverschiebung w und ihr Imaginarteil die Spannungsfunktion ®. Damit sind dann 
aber nach (3), (8), (10), (11) auch alle tibrigen Schnitt- und Verschiebungsgréfen bekannt. 


4, Lésung der Differentialgleichung. Fiir die praktische Rechnung ist es zweckmabig, die 
Koordinaten € = x/l und 9 = ¢/# einzufiihren, wobei | die Tragerspannweite und # den Offnungs- 
winkel der Schale angibt (Abb.1). Die Differentialgleichung (19) ]a8t sich durch den Ansatz 


W—C,,,¢ sinnaé (Gp 13 eos) (20) 


befriedigen, worin C,,, komplexe Konstanten sind. Wird (20) in (19) eingesetzt, so ergibt sich 
die Bestimmungsgleichung fiir m zu 


me —2 (BP) mt 4 (22) (nytt + ENAE— Ho. (21) 


at 


Dies ist eine biquadratische Gleichung, deren Lisung sich leicht angeben JaBt. Es ist 
m= a, + 11, m; = —o1 + ips (22) 
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mit 
o1 = Pi Br=ya- (23) 
Darin bedeutet 
pi a (24) 
| n Vary 21? _eive sal 1 
Pv =p EASE IR BAe ey tee 
; eae ee (25) 
[Fite ee 
Yes | va |/x co ae Re le 
1 ah eee ta - 
ny erl tae (26) 
l /12(1 —p?) 


[4 


& 
BS) 


yd 


~ 
& 


% 


A 


eee 1 all L = 1 L \ 
0 G02 =COK 006 G06 Ow Gn GM G6 Wes 020 O22 G24 


Abb. 3. Diagramme zur Lésung der Differentialgleichung. 


Werden die von den Randern zahlenden Winkel m = “ ie =7- = und wm = 5 (Abb. 1) 
eingefithrt und wird (22) in (20) eingesetzt und schlieBlich Realteil und Imaginarteil getrennt, 
so ergibt sich schlieBlich fiir w und K @ die Darstellung 

w = {e ™° (B, cos fh, @ + By sin B, @) +e *” (Bs cos 8, © + Bysin Bp @) + 
+e “® (B;cos fp, + BysinB,w) +e “° (B,cosB,w + Bysin B,w)}sinnaé, 
K@= {e “” (B, sin B, © — B, cos B, @) +e “* (— By sin fy @ + B,cos By o) + 
+e ™° (B;sin Bo — B,cosB,w) +e “° (—B,sin B,w + By cos fy o)}sinnaé. 


Daraus folgen dann die Schnitt- und Verschiebungsgré8en fiir uw = 0 zu 


N,=— (ya)[e “° (B,sin B, © — By cos B, @) +e “° (— Bysin B, @ + By cos B, a) 


—= HO) . — wo +| (27 
+e “" (B;sin Bj w — Becosf,w) +e “” (— B,sin By o + B,cos B,@)|sinnazé, = ) 


Neg = —yale ™° (Bisin B, @ — Bicos B, @) +e “” (— Blsin B, © + Bicos f, @) 
+e %” (Bisin Bo — BicosB,w) +e ™° (— Bisin B,@ + BicosB,a)|cosnz€, ie 
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N, = ter (Bi sin B, wo — Be cos By @) +e 0 (— B3sin 8, w + Bi cos By @) + (29) 
eee (B? sin B, o — Bzcosf,w) +e “*° (— Bz sin fb, @ + B? cos By w)| sinnwé, | 

2, = ray ie Minas ( B} cos By @ wo + Bs sin B, @ w) )+e ™ e (B3 cos By c w + Bisin By ¢ @ ict | 
er? o (B3 cos B, w + B3sin B, «) vote AX, (BF cos By @ + Bi sin p, w)| — (30) 
mote *° (Bl cos f, @ + Bh sin B, @) +e ™”° (Bleos By @ + Blsin B, @) + | 
+e " (Bicos B, @ + Bisin Bw) +e ™° (Blcos By w + B} sin By «)]| sinnzé, | 

= ja ye ™° (Becos B, @ + B3sin B, @) + e~%9 (B3cos By @ + Besin fy @) + 
a 
+e "®(B?cosf, @ + Bsin Bw) +e ™ > (B? cos By m + B?sin By w) ) — (31) 
( 
—(y x)? le ~ °(B, cos B, @ + Basin B, w) +e %° (Bs cos B, @ + B, sin fp, w) 
+e ™°(B, cos B,@ + B,sinB,w) +e ™“°(B,cosf,w +B, sin By «)|) cosnmé, 

M, = 5 le" (Bheos fy @ + Bisin f, 0) +e (Boos f, 0 + Bjsin fy @ (22 
= aga rn) (B? cos B, o a0 B2sin p, «) eo ® o (B2 cos Bo a) + B? sin B, w) oe nazé 7 : 
M,o= ae? e ® (Bi cos B, @ + Bisin B, 0) +e (Bi cos By @ + Bysinf,@) + — | (33) 
+e ™° (Bcos B,o + Bisinb,w) +e ~*° (Bheos B.w + Bi sin B, w)| cosnmé , | 
—=—_ + aE Cae  (B, cos B, @ + By sin B, w) +e “*°(B,cos By @ + B,sin By @) + (34) 
+e" (B, cos By» + By sin fi, 0) +e (B,cosByo + Bysinfyo)|sinnxé, | 
one = — 28 /e—™.% (BYsin Bi, @ — Bycos fi, 0) +e ™° (— B2sin B, @ + B2cos B, @) + (35) 

+ e-% © (B2sin B, o — Be cos B,) +e “” (— BisinB,o + B3cos fy )| cos navé ; 

Et = we fom: BCE? sin B, © — Bz cos fh, @) + gfe > B2 sin B, w + Bi cos By @) + | (36) 

pe (Pie Bi = Recodo) lee. ° (— B3sin B, w + B2 cos f, w)} sin eae 

iG 5 Sn Py 6 1 
Eww =—#a (Se ™” (Bisin B, © — Bycos By @) + & (— Bysin B,@ + Bycos B, 0) + 
= oo a, @ (BE? sin py oo Bi cos py @) + (es A, w (— B} sin 1035 (00) + B} cos Bo o)| 5 (37) 
— Gaile“ (Bisin B, © — By cos f, @) +e" (— BysinB,@ + Bi cos fy ©) + | 
yn 

6" °(B8sin Bo — Bgcos B, a) +e“ °(— Basin Bw + Bgcos Bw)|j sin nag, | 
a /12 ps AD ; Bast a Sig) B wo + B sin [ w) ++ | 

Etw= (B, cos B, @ + BysinB, w) +e °(B,cos Py o + Bysin Bp (38) 

ae “° (B, cos Bo + B,sin fh, w) +e “° (B, cos By + Bg sin By o)) sin NIE , | 

ie oe oa 12 |e~ "> (Bi cos B, @ + Bisin B, @) te” (Bicos B. @ + Bisin By w) + | (39) 

t 


+6” (Bicos Bi, w + Bisin B, @) +e ~° (Bicos By o + Bysin By) sin nav . | 
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Die darin auftretenden Konstanten B (uw = 0,1, 2,3; t= 1, 2,5. .8) lassen sich aus folgenden 


Rekursionsformeln bestimmen: 


[i 4— 1 u— 1 | =I 
Bi=o,B +8, B, By = 0, B, +f, Bi ; 
; p—l u—l Lu p—l (ei (40) 


Darin ist fiir B} = B; zu setzen. 


5. Zusammenfassung. Die Ermittlung des Spannungs- und Verformungszustandes an den 
Randern belasteter kreiszylindrischer Tonnenschalen wird auf eine komplexe Differential- 
gleichung vierter Ordnung zuriickgefithrt. Der Realteil der Lésungsfunktion liefert die Normal- 
verschiebung, der Imaginarteil eine Spannungsfunktion. Durch diese beiden GréBen lassen sich 
die iibrigen Schnitt- und VerschiebungsgréBen ausdriicken. Die Lésung der Differentialgleichung 
fiihrt auf eine biquadratische Gleichung, deren zahlenmaBige Ergebnisse sich aus Diagrammen 
ablesen lassen. Es werden dann fiir die Schnitt- und VerschiebungsgréBen explizite Formeln 


angegeben. 


(Eingegangen am 21. Mai 1952.) 


Anschrift des Verfassers: Privatdozent Dr.-Ing. W. Zerna, Hannover, EdenstraBe 30. 


‘ 


Einféhrung in die Technische Thermodynamik und in die Grund- 
lagen der Chemischen Thermodynamik, vonDr.-tng. Ernst Schmidt, 0. Pro- 


fessor und Direktor des Instituts fiir Warmetechnik an der Technischen Hochschule 
Braunschweig. Vierte iiberarbeitete und erweiterte Auflage. Mit 244 Abbildungen 
und 69 Tabellen sowie 3 Dampftafeln als Anlage. XVI, 520 Seiten. 1950. 


Ganzleinen DM 30.— 


Die schnelle Entwicklung der Technik in den letzten Jahrzehnten hat eine erweiterte Anwendung d i 
. 5 er Thermodynamik 
ego ener Gebieten zur Folge gehabt. Bei dieser Sachlage entstand der Wunsch nach einer eae toeadics Dar- 
: e ung es-unter dem Namen »angewandte Thermodynamik*‘' stehenden Wissensgebietes. Das vorliegende Buch 
Sore liesem ‘Wunsche Rechnung. Die Tatsache, daB es nach kurzer Zeit bereits in vierter Auflage erscheint, zeigt, 
a8 sein Aufbau dem Bediirfnis der Ingenieure entspricht. In 21 Abschnitten, die reich mit Abbildungen und Tafeln 
ausgestattet sind, ist der Stoff gegliedert. In dem reichen Inhalt sind die Grundlagen besonders entwickelt. Viele 
Beispiele vertiefen die Beherrschung des Stoffes und erleichtern die Lésung der im einzelnen behandelten Probleme, 
Zahlreiche Zahlenangaben erleichtern die Durchfihrung von Rechnungen, so daB das Buch auch in dieser Bezichung 
einen geschlossenen Aufbau zeigt. Von den Sonderabschnitten sei hier nur auf die Behandlung der Strémungs- 
maschinen, der Warmeiibertragung und die Anwendung der Thermodynamik auf chemische Vorgange hingewiesen ... 

» Zeitschrift des Vereins deutscher Ingenicure‘‘ 


Die Grundlagen der angewandten Thermodynamik. Von Dr.-Ing. habil. 


Kurt Nesselmann, Wiesbaden. Mit 311 Abbildungen und 5 Diagrammen im Text. XI, 
320 Seiten. 1950. Z Ganzleinen DM 18.— 


Das Buch fihrt zunichst durch Behandlung der wichtigsten thermodynamischen Grundbegriffe und Grundgesetze 
iiber vollkommene und unvollkommene Gase und Dampfe, iiber die strémende Bewegung von Gasen-und Dampfen, 
iiber den Entropiebegriff und IJ. Hauptsatz der Warmelehre den Leser in das Wesen der Materie ein, um ihn zu 
thermodynamischem Denken zu erziehen. Angesichts der beabsichtigten Beschrankung des Werkes auf einen maSigen 
Umfang bespricht der Verfasser dabei nur das Grundsitzliche der verschiedenen Beziehungen und Prozesse und ver- 
zichtet bewuBt auf die Darstellung technischer Einzelheiten... : 


Die Darstellung des umfangreichen Stoffes ist-trotz mancher Schwierigkeit der Materie und trotz der riumlichen 
Beschrankung, die sich der Verfasser auferlegt, leicht verstandlich, wobei zahlreiche gute Abbildungen, Diagramme 
und praktische Rechenbeispiele wertvolle Dienste leisten. Fir ein eingehenderes Studium sind im Text gentigend 
viele Hinweise auf die neueste Fachliteratur enthalten, so daB das Buch nicht nur Studierenden, sondern auch Inge- 
nieuren, der Praxis zu empfehlen ist, Ausstattung des Werkes ist mustergiiltig. 

»,Heisung — Liftung — Haustechnik‘ 


Leitfaden der Technischen Warmelehre ncbst Anwendungsbeispiclen. 
Von Dr.-Ing. habil. Hugo Richter, Gummersbach. Mit 384 Abbildungen, 1 Diagramm 
und 104 Zahlentafelnim Text und Anhang. XII, 617 Seiten. 1950. Ganzleinen DM 34.50, 


Das gut ausgestattete Werk bringt in seinem Hauptteil die Grundlagen der Technischen Wirmelehre: Die Eigen- 
schaften und GesetzmaBigkeiten der Stoffe und Dampfgemische, die Warmesitze, Zustandsdnderungen und Kreis- 
prozesse, Warmeinhalt und Entropie, ferner Stromungslehre, Warmetibertragung und Verbrennung. Ohne zu hoch 
gegriffene Voraussetzungen wird der Leser mit der mathematischen Behandlung der Probleme vertraut gemacht und 
bis an den gegenwartigen Stand herangefihrt. Das Buch eignet sich daher sowohl zum Selbststudium, dem zahlreich 
eingestreute Anwendungsbeispiele entgegenkommen, wie auch zur Erginzung der Kenntnisse fiir den im Berufsleben 
stehenden Ingenieur, und kann als zuverlassiger Fihrer hierzu bestens empfohlen werden. Durch die zahlreichen 
Tabellen und Diagramme stellt.es zugleich ein wertvolles Handbuch und Nachschlagewerk dar... 
»»Kaltetechnik** 


Warmeibertragung im Gegenstrom, Gleichstrom und Kreuzstrom. 


Von Dr.-Ing. Helmuth Hausen, o. Professor an der Technischen Hochschule Hannover. 
(Bildet Band VIII von ,,Technische Physik in Einzeldarstellungen™.) Mit 230 Text- 
abbildungen. XII, 464 Seiten. 1950. DM 69.—; Ganzleinen DM 73.— 


... Das auBerordentlich Mar und fliissig gesehriebene Buch wird jedem etwas geben: Der Wissenschaftler, der sich 
vor allem fiir die Theorie der Warmetibertragung und der Apparate interessiert, wird nicht nur viele Anregungen 
finden, sondern auch die Eleganz bewundern, mit der hier die Mathematik dienstbar gemacht worden ist. Der Prak- 
tiker wird dagegen angenehm iberrascht sein, wie einfach sich auch die schwiexigsten Aufgaben lésen lassen, ohne 
daB die. Genauigkeit darunter leidet, Ex wird insbesondere eine Reihe von rechnerischen und zeichnerischen Niahe- 


rungsyerfahren finden, die abgewanidelt auch auf andere Probleme anwendbar sind... 
»» Chemie-Ingenieur- Technik‘ 
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Vor kurzem erschien: 


Physikalisches Woérterbuch 


Herausgegeben von Wilhelm H. Westphal, Berlin. Zwei Teile in einem Band. Mit 
etwa 10500 Stichwértern und 1595 Textfiguren. VII, 1628 Seiten. 1952. 
Halbfranz DM 148,— 


Inhaltsiibersicht: Erster Teil: Stichwérter A—L. Zweiter Teil: Stichworter M—Z. 
— Geschichte der Physik. — Lebensdaten der Physiker, — Tabellen. — Nachtrage. 


Das bewahrte Physikalische Handwérterbuch von Berliner und Scheel, das zuletzt 1932 in 
2. Auflage erschien, ist seit langem vergriffen, In den seitdem verflossenen 20 Jahren hat 
die Physik ihren seit der Jahrhundertwende eingeleiteten Fortschritt in eher noch ver- 
starktem Tempo fortgesetzt. Insbesondere ist die Physik der Atomkerne mit ihren welt- 
bewegenden Auswirkungen fast ganz ein Kind dieser jiingsten Zeit. Aber auch auf den 
meisten anderen Gebieten sind wesentliche, teilweise umwalzende neue Erkenntnisse zu 


verzeichnen. Als Beispiel sei nur die Physik des Weltalls genannt, der es mehr und mehr ~ 


gelingt, auf dem Boden der allgemeinen Relativitatstheorie einerseits, der Atom- und Kern- 
physik andererseits, das Werden und den heutigen Zustand des Weltalls und der Himmels- 
korper zu entritseln. \ 


Es war deshalb notig, ein véllig neues Werk zu schaffen, das dieser Entwicklung Rechmung 
tragt. Um Raum fiir den gewaltig angeschwollenen Stoff zu gewinnen, war eine gewisse 
Beschrankung auf das eigentliche Erkenntnisgut und das fiir die reine Forschung Wesent- 
liche unter Verzicht auf die Behandlung rein technischer Anwendungen unumgianglich. 
Trotzdem hat sich die Zah] der Stichwérter mit etwa 10500 gegeniiber dem Physikalischen 
Handwérterbuch nahezu verdoppelt. 


Die Zuverlassigkeit der Darstellung wird durch die Verteilung des Stoffes auf etwa 80 Mit- 
arbeiter gewahrleistet, die als Sachverstandige auf ihren Spezialgebieten bekannt sind. Hin 
Nachtrag bringt Erganzungen, die noch wahrend der Drucklegung eingingen. Ein Anhang 
enthalt unter anderem einen kurzen AbriB der Geschichte der Physik und eine Tabelle mit 
den Lebensdaten von etwa 800 Physikern. 
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Grundlagen der Atomphysik 


Eine Einfihrung in das Studium der Wellenmechanik und Quantenstatistik 


Von Dr. phil. Hans Adolf Bauer, Professor.an der Technischen Hochschule‘und der 
Universitat in Wien. Vierte, umgearbeitete und bedeutend erweiterte Auflage. Mit 


244 Textabbildungen. XX, 631 Seiten. 195). ‘ ' Ganzleinen DM 45.— — 


Inhaltstibersicht: Die Teilchenstruktur der Materie. — Die Wellenstruktur der Materie. 
— Die Vereinigung des Teilchen- und Wellenbildes in der Wellen-( Quanten-)Mechanik. — 
Schrédingers Stérungstheorie. — Relativistische Verallgemeinerung der Wellenmechanik 
(Diracsche Theorie). — Quantenstatistik. — Namen- und Sachverzeichnis. 


. .. Dieses sowohl vom experimentellen, wie vom theoretischen Standpunkt geschrieb=ne, 
leicht faBliche Werk, das auch die neuesten kernphysikalischen Erfahrungen beriicksichtigt, 
ist fur alie geeignet und empfehlenswert, die die fiir eine ernste Beschiftigung mit der Atom- 
physik unerlaSliche Vorbildung besitzen. Das Buch geht so weit, daB es nicht nur dem 
Anfanger auf dem Gebiet der Atomphysik zum Studium dieses Gebietes bestens dienen kann, 
sondern da es auch dem Fortgeschrittenen noch manches zu bieten vermag, was er bisher 
nur in spezielleren Darstellungen finden konnie. 
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